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part  i 

INTRODUCTION    TO   TKEmiODYNAllIC  ANALYSIS 


CHAPTER  1 
INTRODUCTION 

The  United  States  depends  extensively  on  fossil  fuel  imports  for  energy 
and  chemicals.   To  reduce  this  foreign  dependency  and  become  energy  self- 
sufficient  we  must  convert  alternate  energy  resources  to  liquid  fuels  and 
chemical  feedstocks.   Many  forecasters  predict  that  coal,  oil  shale,  and 
tar  sands  among  other  nonrenewable  energy  resources  can  supply  this  country's 
energy  needs  for  several  hundred  years  once  the  technology  is  realized  and 
the  economics  become  attractive.   Ultimately  future  energy  resources  will 
have  to  be  renewable.   A  major  renewable  energy  resource,  biomass,  results 
from  the  utilization  of  plant  material,  an  accumulation  of  solar  energy. 
In  addition  to  these  alternate  energy  resources  the  United  States  must  pro- 
ceed with  a  program  of  technological  improvements  in  energy  conservation 
since  it  is  quite  clear  that  energy  will  never  again  be  -is  cheap  and  abun- 
dant as  in  the  past. 

With  the  rapidly  rising  cost  of  energy  one  area  of  conservation  that 
has  been  receiving  increased  attention  is  improved  waste  heat  recovery. 
The  recovery  of  heat  destined  to  be  dissipated  into  the  environment  yields 
many  benefits  including  a  reduction  in  fuel  consumption  and  reductions  in 
air  and  water  pollution  because  of  the  reduced  heater  and  cooling  tower 
loads.   Thermodynamic  analysis  should  play  an  important  role  in  the  develop- 
ment of  energy  conservation  programs  such  as  this  one  to  correct  current 
wasteful  energy  practices  and  in  the  development  of  alternate  energy  re- 
sources to  help  fend  off  foreign  dependency. 

Currently  the  Chemical  Manufacturing  Association  as  reported  by  Faliwell 
and  Greek  (1)  is  claiming  to  have  reached  their  first  goal  of  a  15%  reduction 
in  energy  use  through  rather  simple  straightforward  conservation  techniques. 


They  have  new  set  an  energy  conservation  goal  of  yet  another  15%  by  1985. 
To  achieve  this  additional  savings  more  sophisticated  conservation  techniques 
must  be  employed  that  require  full-scale  planr-by-plant  energy  audits. 
These  audits  would  determine  which  segments  of  the  processes  are  the  most 
inefficient,  and  they  would  show  where  capital  can  best  be  spent  in  order 
to  realize  a  large  energy  savings.   Second  law  analysis  can  be  used  most 
effectively  to  locate  the  process  segment  inefficiencies  and  determine  which 
are  the  most  inefficient.   It  can  also  be  used  to  compare  different  design 
alternatives  that  correct  these  inefficiencies. 

A  thermodynamic  analysis  of  a  process  can  be  based  upon  either  the  first 
or  the  second  law  of  thermodynamics.   The  first  law  is  an  energy  conserva- 
tion principle  whereas  the  second  law  assigns  a  quality  or  worth  to  energy. 
Each  energy  source  has  a  certain  potential  to  perform  work  related  to  the 
second  law,  i.e.  a  work  equivalent.   This  is  a  function  of  temperature, 
pressure,  chemical  composition,  and  the  environment  in  which  it  exists . 
For  example,  ^  million  kcai  of  thermal  energy  in  a  waste  flue  gas  stream 
at  5Q0°F  is  most  certainly  worth  more  than  the  same  amount  of  thermal  energy 
in  a  wastewater  stream  at  150"F.   In  comparing  the  relative  worths  of  energy 
sources  the  work  equivalent  is  the  best  basis  for  comparison. 

This  idea  may  be  extended  to  form  a  second  law  analysis  thereby  pro- 
viding a  common  basis  in  which  different  process  schemes  can  be  compared 
and  their  inefficiencies  identified.   Although  second  law  analyses  have  not 
been  widely  used,  recently,  several  existing  processes  or  Industries  have 
been  examined  in  this  manner  (2,  3,  4,  5,  6).   The  second  law  analysis  can 
be  used  throughout  the  development  of  a  new  process  to  determine  which 
approaches  utilize  energy  most  efficiently  and  to  identify  areas  of  improve- 
ment.  For  example,  currently  the  subject  of  Gasohol  is  attracting  much 


attention  and  so  are  the  two  processes  fox-  i.-vanufacturir.g  the  required 
ethancl:   the  farming- fermentation  process  and  the  ethylene  hydration  pro- 
cess.  Second  law  analyses  (7)  have  shown  that  the  farming  fermentation 
route  is  roughly  twice  as  energy  efficient  as  the  other  route  and  that  it 
has  a  greater  potential  for  improvement  than  the  other  route  to  ethanol 
production.   In  this  thesis  the  Char  Oil  Energy  Development  (COED)  process 
(8)  which  has  been  researched  extensively  was  chosen  to  illustrate  the  pro- 
cedures and  necessary  considerations  in  order  to  apply  the  thermodynamic 
principles  of  this  second  law  analysis  to  one  alternate  energy  technology. 

Second  law  analysis  may  also  be  used  to  carry  out  an  optimization  of 
a  process  design  from  the  standpoint  of  an  overall  efficiency.   The  concepts 
of  work  equivalents  or  availabilities  which  are  similar  have  recently  been 
incorporated  into  heat  exchange  synthesis  (9)  and  distillation  system  syn- 
thesis (10).  These  system  syntheses,  however,  examine  the  heat  exchange 
network  in  a  general  manner  from  an  ultimately  economic  point  of  view.   Here 
the  optimization  of  a  process  design  will  be  considered  in  detail  from  a 
totally  thermodynamic  point  of  view  wherein  the  process  will  be  designed 
for  a  maximum  in  the  overall  recond  law  efficiency.   That  efficiency  as 
will  be  discussed  shortly  is  distinguished  from  the  second  law  efficiency 
by  the  additional  energy  inputs  and  inefficiencies  included  therein.   In 
this  thesis  a  flue  gas  heat  recovery  process  was  chosen  to  illustrate  this 
thermodynamic  optimization  of  one  energy  conservation  technology. 

Specifically,  the  investigation  herein  will  deal  with  the  thermodynamic 
principles  involved  with  the  analysis  and  optimization  of  processes  with 
the  appropriate  considerations  needed  to  proceed  with  both.   The  fundamental 
principles  involved  with  earlier  and  more  recent  thermodynamic  analysis  will 


be  presented  in  the  remainder  or  Part  I  while  the  second  law  analysis  and 
thermodynamic  optimization  applications  will  be  presented  in  Parts  II  and 
III  respectively.   This  investigation  will  show  that  two  different  effi- 
ciencies may  be  used  in  a  second  law  analysis  and  that  each  cne  has  dis- 
tinct advantages  and  disadvantages.   Also  the  concept  of  an  overall  second 
law  efficiency  will  be  dealt  with.   This  involves  consideration  of  energy 
inputs  and  inefficiencies  not  normally  included  in  an  efficiency  analysis 
or  thermodynamic  optimization.   One  such  additional  energy  input  is  the 
energy  associated  with  the  fabrication  of  plant  equipment  and  with  the 
construction  of  a  plant.   This  factor  has  not  previously  been  included  in 
efficiency  analyses  or  thermodynamic  optimizations .   It  will  be  discussed 
extensively  in  Chapter  S  with  specific  reference  to  its  application  to  the 
COED  plant.   Additional  inefficiencies  that  must  be  considered  in  an  overall 
second  law  efficiency  are  those  associated  with  raw  materials  processing  and 
electrical  generation.   Finally  an  overall  second  law  efficiency  will  be 
presented  for  the  COED  process  and  two  thermodynamic  optimums  will  be 
given  for  two  corresponding  definitions  of  the  flue  gas  heat  recovery 
system.   The  COE'D  second  law  efficiency  will  be  compared  with  the  first 
law  efficiency  of  that  process,  and  one  of  the  thermodynamic  optimums  of 
the  flue  gas  heat  recovery  system  will  be  compared  with  a  corresponding 
economic  optimum  for  that  system. 


CHAPTER  2 
FIRST  LAW  ANALYSIS  OF  CHEMICAL  PROCESSES 

2.1  Derivation  of  the  first  law  efficiency 

This  chapter  deals  with  the  first  law  approach  in  finding  an  effi- 
ciency for  a  chemical  process.   Until  recently  this  has  been  the  primary 
approach  given  in  the  literature  for  evaluating  the  energy  efficiency  of 
any  process.   Because  most  previous  applications  have  lacked  rigor  it 
will  be  advantageous  to  examine  the  thermodynamic  principles  and  to  provide 
the  basis  for  a  general  first  law  efficiency. 

Specifically,  the  first  law  efficiency  may  be  derived  from  the  state- 
ment of  the  first  law  for  an  open  steady-state  system  in  which  potential 
and  kinetic  energies  have  been  ignored: 

&H  =  Q  +  W  (2.1) 

where 

AH  =  total  change  in  enthalpy 

0  ■  net  heat  flow  across  system  boundaries  with  heat 
flow  in  being  a  positive  quantity  and  heat  flow 
out  being  a  negative  quantity 

W  •  net  shaft  work  between  system  and  surroundings  with 
the  convention  that  work  done  on  the  system  is  posi- 
tive and  work  done  by  the  system  is  negative 

The  total  enthalpy  change  may  then  be  expressed  in  terms  of  the  standard 

heat  of  reaction  and  the  appropriate  reactant  and  product  heat  effects: 

iH  "  ~H298  *  AHP  -  iHR  C2.2) 

where 

AH    ■  the  standard  heat  of  reaction  at  298K  and  1  atal 

"H    =  enthalpy  change  as  the  products  in  their  standard 
states  are  taken  from  298K  and  1  atm  to  their 
actual  temperatures 

t&  ■  analogous  enthalpy  change  for  reaccants 
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These   terms   can  be  expanded   to   represent  the   total   enthalpy  change  of  a 

general  chemical   process   involving  any   number   of   chemical  reactions   as 

shown  in  Figure  1.      This  process  has   n  reactant  streams  with  molar  flow 

rates    r    ,    r_,    ...    r.,    ...    r      and  n  product   streams  with  molar   flow  rates 
12  i  a 

p_ ,    p„,    ...    p p     with    the  sum  of   all  heat   flows    equal   to  Q   and   the 

sum  of   all  work  exchanges   equal   to  W.      The  components   of  AK  may  now  be 
expressed  in   the  following  manner: 

f  (2.3.1 


(2.4) 


(2.5) 

where 

AH    =  standard  heat  of  formation  for  component  a  (a  =  i  or  j) 

Cp   =  constant  pressure  heat  capacity  for  component  a  (a  =  i 
or  j)  with  units  consistent  with  those  cf  AH0 

T    =  reference  temoerature  of  298K 
o 

Next,  the  total  standard  heat  of  combustion  is  formed  for  the  reactants 

and  products  in  terms  of  standard  heats  of  formation.   The  chemical  species 

are  taken  to  be  composed  of  the  elements  C,  H,  0,  N,  and  5  and  the  combus- 
tion equation  for  an  arbitrary  species  i  is 

C  H,  0  8.  S,   +  (a,   +  3i   Yi  ,  ,  .  . 

H  si  Yi  *i  Ai    x     T  "  T  +  Ai)02 

8.  5. 

=  Bl«J2  +  -f  H20  +  -f  N2  +  \i  S02  (2.6) 

In  this    equation  H„  was    taken  as    the  oxidation  product   instead  cf  NO 
I  x 

because  of   the  unfavorable  equilibrium  of    the   latter   for  all   Dut  very  high 
temperatures.      The  H„0   is   either   liquid  or  vapor  depending  upon   the  method 


of  calculating  AH  which  will  he  considered  in  section  2,2.   The  standard 
heat  of  combustion  for  co-npour.d  i  is 

iHi  "  aiiHco2  +  lF  42o  +  xiAHso2  "  iHi  (2-7) 

Using    this   equation,    the  total  standard  heat  of   combustion  for  reactants, 
AH    ,    and  for  products,    AK_,    is   found: 

n  ,  n      ■      3  n 

AH*?  =  AH*  I     r.ct.    +  Ah£  E      r .  -f  +  AE„_        I      r   A . 

R  CO.,    .    ,      i   i  Ti.,0    .    .      1    2  SO.    .  i  1 

2   i=l  2      i=l  2   1=1 

-  \     r.AHf  ^3> 
1=1      X      * 

iHP  "  iHC02  ^  *fi   +  &HH20    .^j  if  +  AHS02  jx  Pj\i 

m  f 

-  I      p. AH.  (2.9) 

4.1      J       J 

Elemental  balances   show   that  AH„Q_   can  be   expressed  as    a  combination 

c  c 

of  AH     and  AH    .      The   elemental  balances   are 

n  rn 

I      r.ct.    =      :      p. a.  (2.10) 

1-1     X   X        J-1     J    J 

n  g.         «  S. 

I      r    -f  m     I      p     -J-  (2.11) 

1-1     ^    2         j-1     J     Z 

n  in 

E   r  A     =     Z     p  A  (2.12) 

1-1  J-1     J   J 

Conceptually  these  are  seen  to  be  true  because  the  amount  of  a  combustion 

product  such  as  CO.  produced  from  the  combustion  of  the  reactants  equals 

the  amount  likewise  obtained  from  the  products.   Using  these  equalities 

Eq.  (2.9)  is  subtracted  from  Eq.  (2.8)  giving 

AH=  -  48?  -  I     p  AH.;  -  I     r  AHf  (2.13) 

S     P   J-1  *      J    1-1   i  ~ 


Comparison  with  Eq.  (2.3)  reveals  that 

AH°9g  =  AHp_  -  AHj  (2.14) 

Eq.  (2.2)  then  becomes 

AH  =  AK£  -  AKp  +  AH°  -  AH°  (2.15) 

Now,  this  result  is  combined  with  the  first  law  statement  to  obtain 

a  basic  equality  for  the  energy  out  terms  with  the  energy  in  terms.   First, 

Eq.  (2.15)  is  substituted  into  Eq.  (2.1).   3econd,  Q  and  W  are  replaced  by 

Q,-„>  Q„„,>  w,-„>  and  W    where  for  convenience  these  quantities  are  taken 
in   out   m      out 

to  be  positive.   These  actions  result  in 

W  (2.16) 

out 


AH£  -  AH=  +  AH°  - 

AH°  =  Q.  -  Q  +  W. 
R    m    out    in 

Rearranging  yields 

AH^  -  AH°  -  W  _ 
P     P    out 

-  q    „  -  ah£  -  AH°  -  W. 

xout     R     R    i] 

Qln  (2.17) 

The  equation  is  now  multiplied  through  by  -1  to  give 

-AHf  +  AH°  +  u         +  Q         -  -&K?  +  &H?  +  V.     +0.  (2.18) 

t  P  out  out  R  R  in         'in 

Because  all  heats   of   combustion  are  always   negative  -AHC  will  be  written 

as  an  absolute  value,    | AH    [.      Also,   because  most  processes   occur  with 

T_.    >    298  and  T.    >   298   it  should  be  noted   that  AE°   and  AH°  will   almost 
*  J  PR 

always  be  positive.      This   results    in   the  following  general   first   law 
statement: 

|  AH|  |    +  AH°  +  WQut  +  QQut   =    [AH<|    4-  AH°  +  W^  +  Q^  (2.19) 

Again,    it  should  be  remembered   that  when  using   this   equation  all  Q   and  W 
terms   are   taken  as   positive.      Thus,    each   term  in  Eq.    (2.19)    is    expected 
tn  be  positive  making   each  side  of    the   equation   equal    to    the  sum  of  all 
types   of   energy  associated  with   inputs    or  outputs. 


10 


This  equality  should  serje  well  as  a  basis  for  defining  the  first  law 

efficiency.   This  efficiency  can  be  expressed  as  all  of  the  utilised  energy 

out  divided  by  all  of  the  energy  in.   The  situation  when  all  forms  of 

energy  are  realised  to  their  full  potential  will  be  termed  the  unity  case 

because  for  that  situation  the  first  law  efficiency,  n.   ,  equals  unity. 

It  should  be  noted,  however,  that  if  all  possible  inputs  are  considered  it 

is  not  possible  to  have  a  first  law  efficiency  greater  than  unity.   If  any 

form  of  energy  resulting  from  the  process  is  not  realized  to  its  full 

potential  n,    is  less  than  unitv.   For  example,  some  conditions  that 
1st  -  r 

would  cause  n-    to  be  less  than  unitv  are 
1st 

1.  incomplete  combustion  resulting  in  a  product  such  as  CO  in 
the  stack  gas 

2.  nonrecovery  of  the  sensible  heat,  £H_  ,  of  a  product 

J 

3.  mechanical  inefficiencies  resulting  in  a  heat  flow  to  the 
surroundings 

4.  heat  loss  to  the  surroundings 

2.2  Calcualtion  of  AH 

Bailie  and  Doner  (11)  have  defined  a  first  law  efficiency  whose  numeri- 
cal value  depexids  upon  the  type  of  thermcchemical  data  used  to  evaluate  AH. 
Even  though  their  formalization  was  not  precise  their  definition  is  essen- 
tially the  same  as  the  one  given  in  the  previous  section  wherein  the  first 
law  efficiency  equals  utilized  energy  out  divided  by  total  energy  in.   They 
considered  heat  of  formation,  higher  heating  value,  and  lower  heating  value 
data  in  the  calculation  of  a  first  law  efficiency.   They  pointed  out  that 
the  higher  heating  value  data  are  primarily  used  in  combustion  processes 
and  that  for  most  other  types  of  chemical  processes  the  heat  of  formation 
data  are  used.   They  analyzed  a  simple  conversion  system  utilizing  the 
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reaction  CH,  ■  C  +  2H  .   In  this  example  the  carbon  product  was  considered 
waste.   The  first  law  ef f iciericies  they  calculated  were  0,  0.64,  and  0.60 
corresponding  to  the  use  of  heat  of  formation,  higher  heating  value,  and 
lower  heating  value  data  respectively. 

Using  these  methods  of  calculating  AH  three  different  efficiencies  are 
obtained,  however,  only  two  of  these  are  realistic.   First,  consider  cal- 
culation of  AH  from  heat  of  formation  data.   It  is  possible  to  show  through 
manipulation  of  Eqs .  (2.1),  (2.2),  and  (2.3)  that 

AlC  +  AH°  +  Q    +  W   „  =  AH*;  +  AH°  +  Q.   +  If  (2.20) 

P     P    out    out     R     8    in    m 

Hence,  it  is  thermodynamically  legitimate  to  use  the  above  as  the  basis  of 
an  efficiency,  however,  it  is  not  physically  significant  to  refer  the  en- 
thalpy of  species  to  their  elements  when  the  energy  possessed  by  any 
species  is  realized  through  an  oxidation  process.   Thus  it  seems  more 
realistic  to  refer  specie  enthalpies  to  oxidation  products  by  using  heats 
of  combustion.   In  this  way  ridiculous  and  meaningless  results  such  as 
those  obtained  by  3ailie  and  Doner  will  be  avoided.   Next,  consider  the 
calculation  of  AH  from  heat  of  combustion  data  for  two  general  chemical 
process  cases.   In  the  first  case  the  chemical  process  does  not  involve 
water  in  the  reactant  or  product  streams.   As  pointed  out  by  Bailie  and 
Doner  in  their  example  the  higher  and  lower  heating  value  data  yielded  two 

different  first  law  efficiencies.   This  is  seen  to  be  justifiable  because 

f         f 
as  long  as  the  AH,    and  AH.  terms  in  Eqs.  (2.7),  (2.8),  and  (2.9)  are 

handled  consistently  (both  referring  to  same  H„0  standard  state)  either 

heat  of  combustion  method  can  correctly  be  used  to  calculate  AH.   In  the 

second  case  the  chemical  process  does  involve  water  in  the  reactant  and/or 

product  streams.  Here  in  order  to  correctly  apply  both  methods  of  calcu- 

f        f 
lating  AH,  the  higher  and  lower  heating  value  methods,  the  AHi  „  and  AH. 
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terms   in  Eqs.    (2.7),    (2.8),   and   (2.9)   must  agai^.  be  handled  consistently 
with  respect   to  water's   reference  state,   but   also   the  latent  heat  of  vapor- 
ization of  water  must  be   either  added   to   or  subtracted   from  Eqs.    (2.4)    and 
(2.5)   when  appropriate.      This   again  will   lead   to   two   different  first   law 
efficiencies,   however,    in  both  of   the  cases   considered   the  higher  heat  of 
combustion  method  is   perhaps   preferred  because  it  uses    the  maximum  energy 
obtainable   from  any   stream  in   the  efficiency   calculations.      It  should  be 
noted,   however,    that   the  difference  between   these   two  first  law   efficiencies 
is   expected   to  be  small. 
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CHAPTER  3 
THEORETICAL  APPROACH  TO  SECOND  LAW  ANALYSIS 

3 . 1  Previous  work 

This  chapter  deals  with  the  second  law  approach  in  determining  an 
efficiency  for  a  chemical  process.  Two  approaches  to  this  thermodynamic 
efficiency  analysis  based  on  the  second  law  have  been  given  in  the 
literature. 

Denbigh  (  5 )  delineated  the  first  approach  and  applied  it  in  his 
analysis  of  an  ammonia  oxidation  process.   He  used  the  concept  of  availa- 
bility to  calculate  the  actual  and  the  ideal  work  associated  with  the 
process.   With  this  information  Denbigh  found  what  he  termed  the  thermo- 
dynamic efficiency  for  the  ammonia  oxidation  process.   Specifically,  he 
obtained  an  efficiency  of  0.06  for  that  process.   Denbigh  recognized  that 
for  all  practical  purposes  large  scale  industrial  reactions  could  not  be 
carried  out  reversibly  so  under  that  limitation  an  efficiency  of  unity  was 
not  possible.   He  therefore  proposed  the  concept  of  a  "prescribed  degree 
of  irreversibility"  in  which  chemical  reactions  and  the  mixing  of  reagents 
are  allowed  to  proceed  irreversibly.  The  reactions'  heat  effects,  however, 
are  realized  in  a  reversible  manner  and  all  other  physical  or  mechanical 
operations  are  assumed  reversible.   Using  this  concept  Denbigh  calculated 
a  new  efficiency,  the  practical  efficiency.   He  found  that  efficiency  to 
be  0.11  for  the  ammonia  oxidation  process. 

The  other  approach  was  given  by  Gyftopolous,  Lazaridis,  and  Wldmer 
(4)  who  used  the  concept  of  available  useful  work  to  evaluate  the  poten- 
tial for  more  effective  use  of  fuel.   They  evaluated  the  second  law  effi- 
ciency of  several  basic  industries.   Riekert  (6)  formalized  this  approach 
and  its  environmental  datum.   He  employed  the  available  useful  work,  which 
he  called  the  work  equivalent,  to  calculate  the  "efficiency  of  energy 
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utilization."  Riekert  calculated  that  efficiency  for  an  ammonia  oxidation 
process  similar  to  Denbigh's  and  found  an  efficiency  of  0.16.   He  explained 
the  discrepancy  with  Denbigh  in  the  differences  in  their  approaches  even 
though  he  assumed  slightly  different  operating  conditions.   Specifically, 
Riekert  counted  the  work  equivalent  of  the  product  as  an  output  whereas 
Denbigh  gave  the  product  no  value.   Riekert's  approach  will  be  considered 
in  this  chapter. 

He  defined  the  work  equivalent  as  the  maximum  work  obtainable  when 
a  substance  is  moved  reversibly  into  equilibrium  with  a  specified  environ- 
ment.  The  work  equivalent  may  be  developed  from  the  statement  of  the 
first  law  for  an  open  steady-state  system  in  which  potential  and  kinetic 
energies  have  been  ignored,  Eq.  (2.1).   Since  the  processes  taking  place 
are  reversible 

Q  -  /  TdS  (3.1) 

where 

T  =  heat  transfer  temperature 

dS  -  differential  entropy 

Substitution  into  tne  first  law  statement,  Eq.  (2.1),  yields  an  equation  for 

the  maximum  shaft  work  done  by  the  system  which  after  multiplication  by  -1  is 

W    »  AH  -  /  TdS  (3.2) 

max 

For  the  case  where   the  only   transfer  of  heat   takes  place  reversibly  at 

constant   temperature,    T  ,    Eq .    (3.2)    becomes 

W  =  AH  -    T  AS  (3.3) 

max  o 

If  a  reference  state  with  properties  (T  ,  H  ,  SQ)  exists  from  which  no 
work  can  be  obtained  through  interaction  with  the  environment,  then  the 
work  equivalent,  e,  may  be  defined  as 
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e=H-H  -T   (S-S)  (3.4) 

o    o       o 

For  a  homogeneous  mixture  containing  n.  moles  of  species  i  a  partial  molar 

work  equivalent  can  be  defined 

g  =  I   a.*.  =  H  -  H  -T   (S-S)  (3.5) 

.  i  i       o    o  %     o 

where 

t,    -  H.  -  H  .  -  T   (3.  -  S  .)  (3.6) 

1      1      Ol      O     1      01 

is  the  partial  molar  work  equivalent  of  species  i  in  the  mixture  with  respect 
to  the  environment. 

Riekert  considered  the  environment  to  be  fixed  with  respect  to  temper- 
ature, T  ,  pressure,  p  ,  and  datum  level  comoonents .   These  datum  level 
o  "  o 

components  are  assumed  to  be  available  from  the  environment  in  "unlimited" 
supply  without  any  expenditure  of  work.   This  is  analogous  to  an  infinite 
heat  reservoir  of  temperature  T  .   The  properties  of  the  datum  level  compo- 
nents imply  that  no  work  can  be  obtained  from  them  at  T  and  p  ,  i.e.  -;  -  0. 

o      o 

This  would  prohibit  spontaneous  chemical  reactions  between  these  components 
at  T  and  p  .   In  addition  the  work  equivalent  e.  of  any  compound  i  in 
equilibrium  at  T  and  p  with  the  datum  level  components  is  zero.   It  may 
be  noted  that  the  datum  level  of  a  specific  element  is  the  state  of  that 
element  as  it  exists  in  "unlimited"  supply  in  the  environment.   Another 
property  of  the  datum  level  components  is  that  they  provide  the  set  of 
reference  properties  (H  . ,  S  .)  in  Eq.  (3.6).   These  reference  properties 
are  determined  from  the  manner  in  which  substance  i  is  in  equilibrium  with 
the  environment,  i.e.  by  itself  or  as  part  of  a  mixture  . 
3.2  Two  efficiency  definitions 

Work  equivalents  may  be  calculated  for  all  the  process  streams  and 
utilities  associated  with  a  given  process.   In  order  to  do  this  and  then 
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determine  an  efficiency  a  large  amount  of  information  must  be  gathered  - 
temperature,  pressure,  composition,  and  flow  rates  of  all  streams  entering 
or  leaving  the  process,  the  work  inputs  and  outputs,  and  the  heat  inputs 
and  outputs  with  their  temperature  levels.   Once  this  is  accomplished  either 
one  of  the  two  following  efficiencies  may  be  found. 

The  first  of  these,  the  second  law  efficiency  or  efficiency  of  energy 
utilization  according  to  Riekert,  is  defined  on  a  unit  time  basis  as  the 
total  work  equivalent  output  in  all  the  useful  outgoing  streams  divided 
by  the  total  work  equivalent  input  or 
(  ;.  +  H 

out     out  /,  -,s 

n=7T~Tw—  C3-7) 

in    In 

Hie  W   and  W    refer  to  the  total  in  and  out  work  equivalents  of  utilities 
in      out 

such  as  shaft  work,  steam,  condensate,  and  fuel  gas.   In  this  efficiency 
definition  the  work  equivalent  of  any  stream  is  zero  when  it  is  unuti- 
lized and  it  disappears  into  the  environment.   If  all  possible  inputs  are 
considered  it  is  not  possible  for  the  second  law  efficiency  to  be  greater 
than  unity.   This  efficiency  is  seen  to  be  analogous  to  the  first  law 
efficiency  defined  in  Chapter  2  where  energy  has  been  replaced  by  the  work 
equivalent. 

The  other  efficiency,  Denbigh's  thermodynamic  efficiency,  is  referred 
to  here  as  the  incremental  efficiency.   It  is  defined  on  a  unit  time  basis 
as  a  ratio  of  the  ideal  work  involved  with  the  process  streams'  transfor- 
mations to  the  actual  work  equivalent  that  brought  about  these  changes. 
For  a  work-requiring  process  the  incremental  efficiency  is 

-  E°Ut  "  "ia  (3.8) 


'I      0 

net 


For  a  work  producing  process  the  inverse  of  Eq.  (3.8)  equals  n-,.      As  with 
the  first  and  second  law  efficiencies  the  incremental  efficiency  cannot  be 


17 


greater  than  unity  if  all  possible  inputs  are  considered.   However,  in 
contrast  with  these  two  prior  efficiencies  the  incremental  efficiency  can 
be  less  than  zero.   Negative  values  can  occur  when  many  of  the  product 
streams  are  unutilized  in  a  work  requiring  process. 

While  the  second  law  or  absolute  efficiency  measures  the  effective- 
ness of  the  entire  process,  often  the  large  work  equivalents  associated 
with  high  energy  chemical  streams  can  overshadow  the  physical  work  terms 
which  are  more  often  identified  with  fuel  utilization.   In  this  regard 
Gyftopolous  et  al.  (  4)  have  shown  that  the  absolute  efficiency  for  petro- 
leum refining  is  0.91,  but  the  incremental  efficiency,  which  measures  fuel 
utilization,  is  only  0.09.   This  difference  between  the  two  efficiencies 
will  further  be  examined  in  Chapter  &  in  light  of  a  specific  process 
evaluation. 
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CHAPTER  4 
WORK  EQUIVALENT  CALCULATIONS 

4.1  Datum  level  components 

In  this  chapter  the  method  of  calculating  work  equivalents  in  a  specific 
environment  will  be  detailed.   As  discussed  in  the  previous  chapter  the 
environment  must  be  specified  with  respect  to  temperature,  pressure,  and 
datum  level  components.   The  environment  considered  here  is  fixed  with 
respect  to  these  properties: 

(a)  T  =  298K 

o 

(b)  p  =1.0  atm 

(c)  Datum  level  components: 
N2(g):   x^-0.79 
02(g):   x^-0.21 

C02(g):   XC0,  "  °-000314 

H2°U)=  \0  =   1-° 
CaS04-2H20(s"):   x^.,^  -  1.0 

CaC0Q(s):  x„   „  =1.0 

Riekert    (  b )    considered  stack  gases    to  be  part   of  his    environment  so  he 
used  0.17   as    the   C0„  mole  fraction.      So   as    to  be  more  general   the   environ- 
ment considered  here  does   not  include  stack  gases    thus    the  value   used  above 
is    the  mole  fraction  of   CO,-,    in  air.      The   last   two   components    listed  above, 
gypsum   (CaS0,'2H-0)    and   limestone   (CaCO-J  ,    provide  for   the  calculation  of   a 
work  equivalent   for  sulfur.      The  CaS0^*2H70  was   used  as   a  sulfur  source 
because  of   its   reactive  stability,    i.e.,    a  very   low  Gibbs    free  energy  of 
formation,  -^G^,    and  because   it   exists    in  plentiful  supply.      Likewise  CaCO., 
was    chosen  as   a   calcium  source   for  similar  reasons.      In  general  a  datum 
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level  component  should  be  abundant  in  the  environment  and  should  not  react 

with  other  datum  level  components;  this  implies  a  low  AG..   The  state  of 

aggregation  of  datum  level  components  should  also  be  considered.   If  it  is 

a  liquid  or  gas  the  mole  fraction  is  important.   If  it  is  a  solid  then  its 

interaction  with  water  can  be  important.   For  example,  a  solid  can  become 

hydrated  or  it  can  dissolve.   If  the  latter  occurs  the  component's  datum 

level  state  must  be  aqueous  so  that  solution  affects  can  be  accounted  for. 

4.2  Standard  state  vork  equivalents 

For  the  environmental  properties  described  in  the  previous  section 

the  work  equivalents  of  the  elements  8«>  0„,  H_,  C,  and  S  in  their  pure 

forms  at  T  and  o  ,  the  standard  state,  can  be  calculated.   These  standard 
o     "  o 

state  work  equivalents,  e.'s  of  pure  N_,  0-,  and  CO,  are  equal  to  their 
standard  free  energies  of  unmixing  from  the  environment.   For  example,  sn 
is  calculated  from 

E°  ■  -RT  in  xr  (4.1) 

°2      °     °2 

Using  these  air  derived  values,  H  and  C  can  be  calculated  from  the  H,0  and 

CO,  formation  reactions.   The  change  in  the  standard  state  work  equivalent 

for  these  reactions  equals  the  standard  free  energy  of  formation. 

o  71O    o    1  o      a  ,,  ,, 

6   '  %  "  2  =0,  *  iGf  H,0  (4-2) 


^ 
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£co2  -  £c  -  ;o2  -  iGf  co2  (4-3) 

The.  work  equivalent  for  water  in  Eq.  (4.2)  is  set  equal  to  zero  because 

liquid  water  is  a  datum  level  component. 

o 

To  calculate  E  the  comoound  CaSO,  -  2H„0  is  taken  to  be  the  stable  con- 

figuration  of  sulfur.  Sulfur  is  put  into  this  compound  from  CaCO^,  the 
stable  source  of  Ca.  The  reaction  to  bring  sulfur  or  a  sulfur  compound 
into  equilibrium  with  the  environment  is 
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S03  +  CaCO,  +  2H,0  =  CaS04  •  2H20  +  C0»  (4.4) 

It  is  noted  that  ^j's  are  available  for  all  of  these  components  except  for 
)3,  thus  e° 


SO,,  thus  e<,n     ^^y   ^e  calculated  from  Eq.  (.4.5)  after  the  standard  free 


3 

energy  of  reaction,  AG  ,  is  determined. 

0  0        ,0 

O        O     71  O  0  ?  _   OS1  ,0  ,,  ,, 

£C02  +   eC:a>e<2H20   "    ES03  V$*i  "    2^0      =   ^  (^5) 

Next,  e  is  calculated  from  the  formation  reaction  of  SO- 

5  J 

z°       _  ,«  _  1  c°  =  aG°  (4.6) 

£S03   £s   2  E0,     f  SC3  k   ' 

'Die  standard  state  work  equivalent  of  any  substance  containing  C,  H, 
0,  N,  and  S  can  be  calculated  using  the  previous  results  if  standard  free 
energy  of  formation  data  are  known.  This  is  achieved  through 

EC  HO  HJS,  "  nec  +  2  3eH,  +  2  V%   +   2  S% 

(4.7) 

The  standard  state  work  equivalents  for  Nn,  0,,  H_,  C,  and  S  are  given  in 
Table  1.   It  may  be  noted  that  with  the  proper  selection  of  datum  level  com- 
ponents and  with  the  above  procedure  a  consistent  set  of  e.'s  for  all  the 
elements  should  be  obtainable. 
4.3  Heat  of  combustion  approximations 

For  complex  hydrocarbons  such  as  coal  and  coal-derived  oils  the  method 
for  finding  standard  state  work  equivalents  cannot  be  applied  because  of 
the  absence  cf  AG,  data.   Therefore  the  approximation  that  e.  is  equal  to 
the  standard  heat  of  combustion,  iH293  (H2°f))'  nust  be  considered-   Table  2 
compares  these  two  values  calculated  for  several  groups  of  hydrocarbons, 
including  coal  tar  constituents  (12),  for  which  AG,  data  could  be  found. 
This  table  shows  that  for  hydrocarbons  the  ratio  -e_.  /aH^q,  increases  as  the 
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TABLE  1.   Standard  State  Work  Equivalents 
for  Basic  Substances 


e   (298K,  1  atm) 

Substance 

State 
8 

kcal/knol 

N2 

143 

°2 

8 

932 

H2 

8 

65224 

C 

c, 

graphite 

98112 

S 

c, 

rhombic 

139660 

Table   2.      A  Camparis 0 

a  or   c ,   and 
i 

AH29S   for  '/a 

* 
rious   Compounds 

22 

Compound 
Normal  Alkana  - 

o 

e  . 
"1 

kcal/kmol 

"AH298 
kcal/kmol 

-£°i/iH298   X  10° 

^4 

193420. 

212820. 

93.23 

C2H6 

357040. 

372810. 

95.77 

n"C4H10 

669820. 

687640. 

97.41 

n-C8H18 

1294830. 

1317440. 

98.28 

n-C20H42 

3170980. 

3206750. 

98.88 

Cyclic  Alkane  - 

C6H12 

933600. 

944780. 

S8.82 

trans  -  C^ 

1504690. 

1511790. 

99.53 

** 

Benzoid  Aromatic- 

Monocylic: 

C6H6 

788340. 

789080. 

99.91 

C7H8 

939840. 

943580. 

93.60 

p-xylene   CgH10 

1094970. 

1098280. 

99.70 

Bicvclic: 

tetralin  C     H, 

1358360. 

1357010. 

100.10 

nanththaiene   C.-H- 
10   8 

1259460. 

1249850. 

100.77 

a-methylnaphthalene  C -,  -,  H ,  „ 

1412390. 

14040S0. 

100.59 

biphenyl  C.JL. 

1525400. 

1513720. 

100.77 

Tr  ley  lie: 

anthracene   C, ,H,„ 
14   10 

1733930. 

1713270. 

101.21 

phenanthrene   C,,H,_ 
14   i.0 

1728290. 

1707730. 

101.20 

Oxygenated 

phenol  C,H,0 

0    0 

749490. 

746230. 

100.44 

p-cresol  C-H„0 

/     0 

904300. 

901660. 

100.29 

** 

Heterocyclic  Aromatic-- 

pyridine 

676650. 

674550. 

100.31 
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Table  2  -    Continued 

Compound 
Biochemicals  - 

0 

£i 

kcal/  Kinoi 

-AH298 
kcal/kmol 

"Ei>H298   X  10° 

Amino  Acids : 

L- Alanine   CJMLH 

403340 

387210 

104.17 

L-Aspartic  Acid   C.H-0 

4N 

415730 

372720 

108.62 

L- Cysteine  CoH-O^NS 

549300 

532200 

103.21(118.3)* 

Glycine  C2H502N 

248870 

232600 

106.99 

L-Glutamic  Acid  CgH„0 

4N 

540700 

537450 

106. 19 

L-Leucine  Ci,!)^ 

871380 

856090 

101.79 

L- Tyrosine  CgH     0 ,N 

1096070 

1058450 

103.55 

L- Tryptophan  C.Ji   „07 

N 

'2 

1515340 

1472870 

102.88 

Monosaccharides : 

D-Glucose  C,H     0 

0     12    6 

711210 

669950 

105.16 

Galactose   C,H. -0, 
6   12  6 

711210 

669950 

106.16 

Disaccharides : 

B-Lactose  C'fl-.O-ii 

1426410 

1346000 

105.97 

Maltose  CjJSmOj, 

1388330 

1349300 

102.89 

Sucrose  C^H^ 

1431750 

1348200 

106.20 

*      298K,    1  atm,    ideal  gas 

**      aromatic  constituents 

of   coal   tar    (12) 

***     &H;L.   calculated  from 

iHf    = 

-127880  kcal/kmol  with   pre 

iducts   of  N_, 

H2°(£)'    C02'    and  S0-> 

Thermodynamic   data  - 

(13), 

(14),    (15), 

and    (16) 
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carbon  to  hydrogen  ratio  increases  and  in  most  cases  as  the  molecular  weight 
increases.   These  results  will  be  taken  as  empirical  and  no  explanation 

will  be  offered.   Therefore,  in  the  light  of  these  results  the  assumption 

o  c 

that  Ej  equals  -AH^nq  will  be  used  in  the  case  of  highly  aromatic  coal  pro- 
ducts.  This  application  to  coal  seems  reasonable  since  Whitehurst  (17)  has 
suggested  that  bituminous  coals  have  an  aromaticity  of  approximately  60% 
and  contain  significant  amounts  of  polycylic  aliphatic  rings. 

Another  class  of  compounds  of  interest  are  the  biochemicals .   They,  in 
the  form  of  biomass,  may  be  important  in  alternate  energy  conversion  tech- 
nologies therefore  second  law  analysis  is  appropriate  and  consequently 
standard  state  work  equivalents  will  be  needed.   Since  W3-  data,  and  for 
that  matter  AH^qq  data,  are  seldom  available  several  approximations  are 
needed.   First,  the  heat  of  combustion  approximation  will  be  considered. 
The  latter  part  of  Table  2  deals  with  applying  this  approximation  to  some 
biomass  related  compounds  for  which  thermodynamic  data  are  available.   The 
work  equivalent  -  heat  of  combustion  comparison  ratios,  -e  /AH „Q  's,  for  the 
biochemical  compounds  are  seen  to  be  slightly  larger  than  those  for  zhe 
hydrocarbon  compounds.   Also,  these  data  suggest  that  in  actual  biochemical 
calculations  it  mav  be  desirable  to  let  £.  equal  1.05  AK^".-  rather  than 
equating  them  one  to  one.   It  should  be  pointed  out  that  this  work  equivalent  - 
heat  of  combustion  approximation  has  been  established  for  compounds  containing 
C,  H,  0,  and  N.   The  assumption  may  not  be  valid  for  other  types  of  compounds. 
For  instance,  look  at  this  assumption  as  it  applies  for  carbon  and  sulfur. 
The  combustion  reactions  and  heats  of  combustion  for  these  elements  are 

C  +  02  =  C02       "AH298  =  94°51  kcal/'k3 
S  +  0,  =  S02       -AH^98  =  70960  kcal/kg 
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For  comparison  the  standard  state  work  equivalenr.3   for  carbon  and  sulfur  are 
98112  and  139660  kcal/kg  respectively.      Thus    for   carbon   the   comparison  ratio 
is   about  one  but   for  sulfur   it   is   about   two.      Clearly,    the  work  equivalent  - 
heat  of   combustion  approximation  is    invalid   for  sulfur.      This   may  be  ex- 
plained  in   that  sulfur's   combustion  product,    SO,,    is   a  datum  heat   of   com- 
bustion product,    zero   energy  value,    but  it   is   not  a  work  equivalent  datum 
component.      Therefore,    the  addition  of  sulfur   increases    the  heat  of   combustion 
a  specific  amount,   but   it  increases   the  work  equivalent  by  more  than  that 
amount  since  a  significant  amount   of  work  can  be  obtained   from  SO-   as    it   is 
moved  reversibly   into   equilibrium  with   the   environment.      It  may  be  noted   that 
for  this    reason   the  comparison  ratio   for  sulfur-containing  L-cysteine   is 
somewhat  higher   than  that  for  other  amino  acids  when  iH298   calculated  from 
AH     data  with  standard  combustion  products    is   used  instead  of    the  literature 
value  of   the  heat  of   combustion. 

The  second  approximation  concerns    estimating   the  heat   of   combustion. 
To   discuss    this   it   is   instructive   to   consider   the  problem  of   estimating    the 
heat  of   combustion   for  wheat.      Minkevich,    Eroshin,    and  coworkers    (18)    found 

that   based  on  the  following   reaction 

Y,  j 

CH  ON     +  ~  0.   =   CO,  +  -   (p  -    3q)H_0  +  qNH,  (4.8) 

p  n  q         4       2  11  2  j 

the  heat  of  combustion  of  a  wide  variety  of  organic  material  is 

AH!j„0  =  27v  kcal/g  mole  C  (4-9) 

29  o      b 

where 

Y,  ■  4  +  p  -  3q  -  2n 
b 

Thus,  if  the  elemental  composition  of  wheat  is  known  AH„  „  can  be  estimated. 

However,  in  the  likely  event  that  the  composition  of  wheat  is  unknown  further 
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approximations   must  be  made.      Minksvich  et  al.    also   found   that   for  most 
biomass   Yh   is   about  4.29   and   that    the  weight   fraction  of   carbon  is   about 
0.462.      This  weight   fraction  corresponds    to   0.0385  g  mole   C  per  g  biomass 
and  results    in  AH  «   2025  kcal/lb    (4460   cal/g)    for  biomass.      Therefore, 

this   AH„Qq  value  estimates    the  heat   of   combustion  for   dry  wheat.      For 
comparison  Merrill  and  Watt    (19)    report   a  value  of   4.03  kcal/g  or  1828 
kcal/lb   for  wheat  containing   10.8%  moisture.      If    the  heat   of    combustion 
estimate  for  dry  wheat  is   corrected  for  10.8%  moisture   then   this    estimate 
becomes   1806  kcal/lb  which  is  very  close  to   the  reported  value.      It  would 
appear   that   the  approximation  of  AH,gg  =   2025   kcal/lb    (4460   cal/g)    for 
biomass    is   reasonable;    in  the   absence  of   other   data   this   combined  with   the 
empirically   established  equality   of   e,    and  AH9c.q   at  best  allows    e.    to  be 
estimated   for  biomass. 
4.4      Total  work  equivalent 

To  arrive  at  the   total  work  equivalent   for  a  substance  both   the  standard 
state  work  equivalent,    £ . ,    and   the   change  in  work  equivalent,    As    ,    to   another 
state    (T,    p,    x  )    must  be  calculated.      Using   the  molar  work  equivalent  as    an 
approximation   for   the  partial  molar  work  equivalent  in  Eq.    (3.6),    Ac's 
follow  from  that   equation's   partial   derivatives.      In   the   following   deriva- 
tions   the  ideal  gas   and   the   ideal  solution  assumptions    are  used.      First, 
take    the  partial  derivative  of  Eq.    (3.6)    with  respect    to   temperature   to 

obtain 

3e.         3H.  ?S. 

~rzr  -•  -^r   -  T  — —      (constant  o  and  x.)  (4.10) 

jT    3T     o  ol  i 

Next,  apply  the  constant  pressure  heat  capacity  relations  to  arrive  at 

3E1        To 
-TT  =  (1  -  -S)    Cp.      (constant  p  and  x.)  (4.11) 

oL  J.      X  1 
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This  equation  may  be  integrated  with  a  two  parameter  heat  capacity  equation, 

CPj  ■  a.  +  b.T,  from  T  to  T  to  give 
ri    i    1  •       o 

Ae.      =    (a.    -  b.T  )(T  -   T  )    +  •§  b,    (T2  -   T   2)    -   a.T  2n  T/T  (4.12) 

iiio  o2i  o  10  o 

T 

For  constant  heat   capacity   Eq.    (4.11)   becomes 

Ae,      -  Cp.    (T  -   T     -   T      In  T/T  )  (4.13) 

iT  *t  o  o  o 

Now,  take  the  partial  derivative  of  Eq.  (3.6)  with  respect  to  pressure  to 
find 

3e.    3H.      3S. 

— i  ■  — i  -  T  -S-*  (constant  T  and  x.)  (4.14) 

3p    3p     c  ap  1 

Then  apply  the  proper  thermodynamic  identities  to  obtain  these  results 

U  3V.  3V. 

<~W\,*A  m  vi  -  t(Tt%,x.  +  To  hff*»A  C*.l» 

»S  sv. 

C-rf)  ■  v.  -  (T  -  T  )   frrffl  W.16J 

3p  i  o         3T    p,x. 

This  equation  may  be  integrated  from  p   to  p  at  T  to  give  for  an  ideal  gas 

A£i  =  R  T  in  o/p  (4.17) 

p      o    "   o 

Finally,  consider  the  change  of  work  equivalent  at  constant  T  and  p  due  to 

a  change  of  mole  fraction.   For  simplification  let  this  change  only  take 

place  at  T  and  p   then  At.  equals  the  change  in  free  energy  or  in  the  case 

of  an  ideal  solution 

Ae.   =  R  T  In   x.  (4.18) 

1         0       1 

x 

Since  the  work  equivalent  is  a  state  function,  its  calculation  based  on  the 
datum  can  be  performed  in  the  following  manner: 
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State  3 
T   P 


adjust 


o 


-> 


constant 
P 
x. 


•Ae._ 


State   2 

adjust 
PH-P0 

State   1 

unmix 

x.->l 

i 

State  C 

X.      E  . 
1        * 

To    Po 

x.      e1 

X            1 

To    Po 

-  ^        ° 

constant 

To 

X. 

l 

cons  tant 

To 
Po 

i         i  | 

Ae.  — > 
l 
P 


-Ae. 


Ae.  =  Ae.  +  Ae.  +  Ae. 


(4.19) 


In  evaluating  a  process  stream  of  materials  the  total  work  equivalent 

can  be  divided  into  a  chemical  work  equivalent.  E  ,   ,  and  a  physical  work 

chem 

equivalent,  e    .   The  chemical  work  equivalent  can  be  defined  as 
phys 


"chem 


(4.20) 


where  n.  equals  the  moles  of  species  i  in  a  stream  on  a  unit  time  basis.   In 

other  words,  e  ,    represents  that  part  of  the  total  work  equivalent,  e. 
chem  i        ,   , 

which  is  due  to  chemical  energy  at  T  and  p  .   The  physical  work  equivalent 
can  be  defined  as 

(4.21) 


'phys 


n.  (Ae,   +  As.   +  Ae.  ) 
i    i_     i      1 

T      p      x 


i.e.  s  ,    is  the  part  of  z   resulting  from  a  state  different  than  the  datum, 
pays 

In  evaluating  the  process  utilities  the  work  equivalents  for  steam, 
condensate,  and  fuel  gas  can  be  calculated  in  the  same  way. 


PART  II 

THE  ENERGETICS  OF  THE  COED  PROCESS 
FROM  A  SECOND  LAW  PERSPECTIVE 
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CHAPTER  5 
COED  PROCESS  DESCRIPTION 

Part  II  is  an  application  of  the  thermodynamic  analysis  presented 
in  Part  I.   Specifically,  the  Char  Oil  Energy  Development  process  of  PMC 
Corp.  (8)  developed  under  the  sponsorship  of  the  United  States  Energy 
Research  and  Development  Administration  will  be  analyzed  from  a  second 
law  perspective  with  the  first  law  approach  used  as  a  comparison.   Both 
the  second  law  and  the  incremental  efficiency  will  be  used  to  show  that 
each  has  distinct  advantages  and  disadvantages.   Also  the  concept  of  an 
overall  second  law  efficiency  which  includes  an  energy  debit  for  plant 
equipment  will  be  dealt  with.   This  factor  has  not  previously  been  included 
in  efficiency  analyses.   The  energy  associated  with  manufacturing  plant 
equipment  and  with  plant  construction  will  be  discussed  and  specifically 
calculated  for  the  COED  plant.   Also  the  inclusion  of  outside  inefficiencies 
involved  with  raw  materials  processing  and  electricity  generation  will  be 
discussed. 

However,  the  process  must  first  be  described  before  the  analysis  may 
proceed.   The  COED  process  fits  into  the  general  category  of  liquefaction 
via  pyrolysis.   Figure  2  is  a  schematic  diagram  illustrating  the  commercial 
plant  design  for  processing  25,000  tons  (22,680  metric  tons)  of  Illinois 
No.  6  -  seam  coal  per  day.   Following  coal  preparation,  the  coal  is  dried 
in  a  fluidi zed-bed  dryer  at  350°F  (450K)  using  the  first  stage  pyrolysis 
off-gases.   The  dried  coal  is  then  fed  through  four  stages  of  pyrolysis, 
each  at  succeedingly  higher  temperatures;  550"F  (561K) ,  850°F  (72SK) , 
1050°F  (839K),  and  1550°F  (111SK) .   The  fourth  stage  pyrolysis  is  affected 
by  injection  of  steam  and  oxygen.   The  char  leaving  the  fourth  stage  passes 
directly  to  the  air  blown,  Winkler  type  fluidized-bed  gasifiers.   Gas 
produced  there  consists  mainly  of  N  ,  CO,  and  H_.   This  gasifier  gas  then 
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goes  through  the  acid  gas  removal  scrubber  and  becomes  a  high  energy  pro- 
duct.  The  pyrolysis  gas  consists  mainly  of  H_Q,  CO-,  CO,  H_,  and  hydro- 
carbons.  It  goes  through  a  recovery  section  in  which  the  pyrolysis  oil  is 
condensed  and  separated  from  the  product  gas.   Following  oil  filtration, 
the  oil  is  upgraded  to  a  synthetic  crude  oil,  syncrude,  by  hydro  treating. 
The  hydrogen  to  do  this  is  produced  by  gas  reforming  a  portion  of  the 
deacidized  pyrolysis  product  gas  followed  by  shift  conversion  and  methana- 
tion.   In  the  gas  clean-up  sections,  sulfur  is  recovered,  via  a  Claus 
process,  as  a  byproduct. 

The  COED  commercial  plant  design  also  contains  various  systeins  to 
satisfy  the  utility  requirements.   Included  are  the  cooling,  boiler,  and 
process  waste  water  systems.  Waste  heat  boilers  produce  a  net  excess  of 
steam  primarily  from  the  compression,  hydro  treating,  and  gasification  units. 
The  fuel  gas  needed  mainly  by  the  pyrolysis  and  hydrogen  units  is  taken 
from  the  product  gas.   The  net  inputs  and  outputs  of  quantities  used  in 
this  COED  process  are  given  later  in  Figure  5. 

This  COED  commercial  design  was  based  upon  pilot  plant  studies  for 
only  the  pyrolysis,  filtration,  product  recovery,  and  hydrotreating  sections. 
The  designs  for  other  sections,  most  of  which  are  well-known,  were  obtained 
from  various  companies.   A  very  simplified  approach  was  taken  in  process 
water  treatment  that  may  not  represent  an  actual  design. 
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CHAPTER  6 
SECOND  LAW  ANALYSIS  OF  THE  COED  PROCESS 

6.1   Computational  scheme 

In  this  chapter  the  COED  process  will  be  analyzed  with  respect  tc  the 
second  law  of  thermodynamics •   A  specific  unit  of  this  process,  the  char 
gasification  unit,  will  serve  to  illustrate  the  detailed  computational 
approach.   However,  first  the  general  computational  scheme  will  be  outlined 
along  with  several  specific  assumptions. 

The  bulk  of  this  analysis  of  the  COED  process  was  accomplished  through 
the  use  of  a  FORTRAN  computer  program  developed  to  analyze  a  general  coal 
conversion  process.   With  only  minor  modifications  any  chemical  process  may 
be  analyzed  througn  its  use  as  long  as  thermodynamic  data  are  available  for 
calculating  standard  state  work  equivalents  for  all  components  and  as  long 
as  the  process  is  fully  described  with  respect  to  stream  analysis,  tempera- 
ture levels,  pressures,  and  utilities.   The  computer  program  checks  the 
material  balance  for  each  unit,  calculates  enthalpies  and  work  equivalents 
for  each  stream  and  utility,  checks  the  energy  balance  for  each  unit  and 
calculates  first  and  second  law  efficiencies  for  each  unit  and  the  overall 
process.   The  first  law  efficiencies  are  calculated  using  higher  heating  value 
data  since,  as  mentioned  in  section  2.2,  this  represents  the  maximum  energy  ob- 
tainable froma stream.  The  second  law  efficiencies  are  calculated  after 
Eq.  (4.7)  and  Table  1  are  used  to  construct  a  table  of  standard  state  work 
equivalents  for  compounds  present  In  the  COED  process.   These  are  presented 
in  Table  3.   A  flowchart  of  this  computer  program  is  presented  in  Figure  3. 
Appendix  A  lists  the  program  and  a  final  output,  and  also  contains  a  detailed 
block  diagram  of  the  COED  commercial  process  with  all  interconnecting  streams. 
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Table  3.   c   for  Compour 

i  . 

ids  in  the  COED  Frocess 

0 

Compound 

(koal/kmol) 

H2 

56224 

°2 

932 

N2 

143 

co2 

4776 

CO 

65770 

CH4 

198420 

C2H4 

308672 

C2H6 

357036 

C3H6 

448018 

SH8 

513613 

-c4Hio 

669  814 

H2S 

187992 

NH3 

59657 

S 

139660 

Coal 

6667 

* 
Char 

6167 

* 

Syncrude 

10611 

* 
estimated  as  equal  to  the  standard  heat  of  combustion, 

value  in  kcal/k.g 

c 
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Figure  3.  Flowchart  of  the  COED 
Thermodynamic  Analysis 
Program 
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Several  assumptions  were  made  enroute  to  the  final  analysis.   First, 
to  be  consistent  all  of  the  stream  and  utility  data  were  used  as  given. 
However,  one  additional  stream  was  created  as  an  output  from  the  hydrogen 
plant  unit  to  account  for  a  discrepency  in  the  material  balance.   An  ele- 
mental balance  dictated  that  this  stream  contain  primarily  carbon  dioxide 
and  water,  thus  it  was  considered  a  waste  gas  stream.   Second,  an  average 
saturated  condition  was  assumed  for  all  utility  steam  used  and  produced. 
Third,  the  physical  properties  of  the  derived  coal  liquids  were  approximated 
by  those  of  anthracene  oil  (13).   Finally,  the  heating  values  for  char  and 
unhydrotreated  oil  which  were  not  given  in  the  commercial  process  description 
were  determined  from  the  range  of  values  listed  in  a  COED  pilot  plant  report 
(20)  and  from  energy  balance  results  of  the  units  involved.   These  values, 
of  course,  do  not  enter  into  the  analysis  of  the  overall  COED  process.   These 
and  other  les3  important  assumptions  are  detailed  in  the  computer  program 
comments  in  Appendix  A. 
6.2  Char  gasification  unit  example 

The  char  gasification  unit  is  one  of  the  major  units  of  the  COED  process. 
Figure  4  gives  a  diagram  of  the  input  and  output  process  streams  with  utili- 
ties for  this  process  unit.   The  essential  data  are  presented  in  Tables  4 
and  5.   Each  stream's  chemical  and  physical  work  equivalent  is  calculated  from 
Eq.  (4.20)  and  Eq.  (4.21)  respectively.   The  results  of  these  calculations 

are  given  in  Table  6.   The  negative  value  of  £  ,    for  stream  6  occurs  because 

phys 

the  stream  is  a  mixture  and  therefore  the  ability  to  perform  work  according 
to  Eq.  (4.18)  was  lost  when  pure  gases  were  mixed  irreversibly.   Application 
of  the  second  law  efficiency  definition,  Eq.  (3.7),  produces 
-  -  6 steam =Q_81 


El  +  E2  +  S3  +  H  +   E5  +  £ 


J. Coal  Fines 


2   Char 


3   Water 


4   Air 


5   Steam 
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Figure  4.   Char  Gasification  Diagram 
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Table    5.      Char   Gasification  Utilities 


Utility 

Value 
(kcal) 

Work 

2.533  x  107 

**Cooling  Water 

2.838  x  107 

***Steam 

4.758  x  108 

**  AT  =  16.7°C 

***  avg.  conditions  138°C,  3.4  atm 
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Table  6.   Work  Equivalents  for  Char  Gasification 


Stream 
No. 

Component 
Coal 

1 

chem 
(kcal) 

.875   x  108 

"phys 
(kcal] 

(kcal) 

1 

2.160 

X 

io3 

1.375 

x  10® 

2 

Char 

3 

.020   x   109 

4.670 

X 

107 

3.067 

Q 
X    10 

3 

Water 

— 

7.021 

X 

io4 

7.021 

4 
x  10 

4 

Air 

— 

1.727 

X 

io7 

1.727 

x   10 

5 

Steam 

— 

3.032 

X 

10  7 

3.032 

X   10 7 

6 

Product  Gas 

2. 

.622   x  109 

-3.690 

X 

IO8 

2.585 

x  IO9 

7 

Waste  Water 

— 

1.395 

X 

105 

1.395 

x   IO3 

8 

Ash 
Work 

: 

8.081 

X 

IO5 

8.081 
3.167 

x  IO3 
x   IO7 

Cooling  Water 

— 

7.651 

X 

io3 

7.651 

x   IO3 

Steam 

— 

1.213 

X 

io8 

1.213 

x  IO8 
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In  this  efficiency  streams  3,  7,  and  8  plus  the  cooling  water  do  not  repre- 
sent useful  work  equivalents.  Stream  4,  air,  is  included  because  expansion 
work  is  possible.  The  incremental  efficiency,  Eq .  (3.8),  is  calculated  for 
a  work  producing  process  (remember  the  convention  that  work  out  is  negative)  as 


cw    steam  _  . , 


'6 


+  -7   +  E8    (£1  +  £2  +  E4  +  E5): 


These  two  efficiencies  are  comparable  to  values  given  by  Gaggioli  (21) 
(22)  for  a  Koppers-Totzek  gasifier  with  coal  preparation.   He  found  similar 
second  law  and  incremental  efficiencies  to  be  0.78  and  0.15  respectively. 
For  comparison  a  first  law  efficiency  is  calculated  as  0.94  for  the  gasifi- 
cation unit. 

This  char  gasification  example  may  also  be  used  to  illustrate  Denbigh's 
concept  of  a  "prescribed  degree  of  irreversibility"  discussed  in  section  3.1. 
In  utilizing  this  concept  the  criteria  for  comparison  are:   chemical  reactions 
are  allowed  to  proceed  irreversibly  with  the  heat  of  reaction  recovered  as 
mechanical  work  in  a  Carnot  type  heat  engine,  and  all  other  parts  of  the 
process,  except  mixing,  are  carried  out  reversibly.   For  the  char  gasifica- 
tion example  the  gasification  reactions  proceed  irreversibly  at  approximately 
1500°F  (1090K)  with  the  heat  of  reaction  partially  recovered  as  low  tempera- 
ture steam  and  the  rest  of  it  lost  to  the  atmosphere  and  cooling  water.   The 
pumps  and  compressors  are  80%  efficient  in  converting  electricity  to  shaft 
work.   To  apply  Denbigh's  concept  the  heat  of  reaction  becomes  totally  recov- 
erable as  saturated  steam  at  1500°F  and  the  mechanical  equipment  is  considered 

reversible,  100%  efficient.   Thus  a  new  work  equivalent  is  obtained  for  the 

8  7 

steam,  2.15  x  10  kcal,  and  a  new  shaft  work  value  is  found,  2.533  x  10  kcal. 

Using  these  work  equivalents  the  two  maximum  efficiencies  obtainable  under 
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the  prescribed  conditions  are  0.342  and  0.26.5  for  the  second  law  and  incre- 
mental efficiencies  respectively.   It  should  be  noted  that  in  calculating 
the  latter  efficiency  the  work  equivalent  of  the  cooling  water  is  zero  because 
all  of  the  heat  of  reaction  was  recovered  as  high  temperature  steam.   Since 
these  are  maximum  efficiencies  it  may  be  helpful  to  compare  these  with  the 
actual  efficiencies  by  defining  a  practical  efficiency  whereby  it  equals 
the  actual  efficiency  divided  by  the  maximum  efficiency  under  the  prescribed 
conditions.   Thus  the  practical  second  law  and  incremental  efficiencies  are 
0.96  and  0.49  respectively.   The  large  practical  second  law  efficiency  is 
the  result  of  the  large  chemical  work  equivalents  of  the  char  and  gas  over- 
shadowing the  smaller  utility  work  equivalents. 
6.3  Overall  COED  Results 

An  overall  flow  diagram  of  the  COED  process  is  given  in  Figure  5,  and 
the  corresponding  stream  and  utility  work  equivalents  are  presented  in 
Table  7.   The  second  law  efficiency  of  the  overall  process  is 

5    6    7    steam         ~ 


(«-  +  e,  +  e,  +  e  +  s   ,    ~  ) 

1    2    3    w    condensate 


This  value  is  slightly  higher  than  the  0.68  given  by  C-aggloli  for  Koppers- 
Totzek  coal  gasification.  The  first  law  efficiency  of  the  overall  process 
calculated  here  is  0.80  while  the  COED  designers'  calculated  value  was 
0.75.   However,  in  their  calculation  they  considered  electrical  work  to  be 
a  negative  output  term  rather  than  a  positive  input  term  which  was  used 
here  in  order  to  remain  consistent  with  the  first  law  efficiency  development 
in  Chapter  2.   Table  8  gives  the  first  and  second  law  efficiencies  plus  the 
energy  balance  closure  for  each  unit.   It  is  observed  that  for  the  overall 
process  and  for  the  units  second  law  and  first  law  efficiencies  are  quite 


42 


1 

Feed   Coal 

2 

Oxygen 

3 

Air 

v 

4 

Water 

P 

CHAR  OIL  ENERGY 
DEVELOPMENT  PROCESS 


"T"~   X 


Syncruae 2_ 


Sulfur      6 


Product  Gas  ' 


t      t 


Work  Condensate  Cooling    Steam 
Water 


Figure  5.  Overall  COED  Process  Diagram 
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Table  7.   Work  Equivalents  for  Overall  COED  Process 


Scream 
No. 

Component 

(k, 

e 
:al) 

1 

Coal 

6.24 

x  109 

2 

Oxygen 

6.06 

x  106 

3 

Mr 

2.27 

x  107 

4 

Water 

5 

Syncrude 

1.69 

9 
x  10 

6 

Sulfur 

1.34 

x  108 

7 

Product  Gas 

Cooling  Water 

3.03 

9 
x  10 

Steam 

3.22 

x  107 

Work 

2.57 

x  108 

Condensate 

3.95 

x  105 
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close.   The  reason  for  this  is  that  the  dominant  terms  in  the  efficiency  cal- 
culation are  those  which  represent  chemical  energy  and,  as  has  been  shown, 
the  heat  of  combustion  and  the  work  equivalent  are  practically  equal. 

la  Table  8  the  second  law  efficiency  is  less  than  the  first  law  effi- 
ciency for  most  units  and  overall,  as  would  normally  be  expected.  However, 
for  four  units,  three  of  which  primarily  involve  separations,  the  reverse 
is  true.  This  occurs  for  the  three  units  involved  with  separations  because 
the  second  law  analysis  takes  into  account  the  work  of  separation.   For 
example,  the  removal  of  CO2  and  HoS  in  the  gasifier  acid  gas  removal  unit 
increases  the  mole  fractions  of  the  other  components  thereby  increasing 
the  work  equivalent  of  the  output  product  gas  stream.   The  large  difference 
between  the  first  and  second  law  efficiencies  in  the  pyrolysis  acid  gas 
removal  unit  is  in  part  explained  by  the  previous  discussion,  but  is  perhaps 
primarily  caused  by  a  non  closing  energy  balance.   The  balance  accounted  for 
only  0.91  of  the  energy  input  to  that  unit.   The  reason  for  this  speculation 
is  that  if  the  lost  energy  was  recovered  for  example  in  the  form  of  steam 
the  first  law  efficiency  would  gain  as  additional  output  one  kcal  for  every 
kcal  of  recovered  energy  whereas  the  second  law  efficiency  would  gain  as 
additional  output  only  about  0.2  kcal  of  work  equivalent  per  kcal  of  recovered 
energy.   Therefore  the  first  law  efficiency  would  increase  more  noticeably 
than  would  the  second  law  efficiency.   The  other  unit  in  which  the  second 
law  efficiency  is  greater  than  the  first  law  efficiency  is  gasifier  gas  com- 
pression.  This  occurs  in  part  because  of  the  condensation  of  water  from  the 
gasifier  gas,  a  separation  effect,  but  it  occurs  primarily  because  the  second 
law  analysis  places  value  on  the  expansion  work  that  may  be  obtained  from 
a  stream  not  at  the  standard  pressure  and  the  defined  first  law  efficiency 
does  not.   Without  the  separation  and  pressure  affects  the  second  law  effi- 
ciency becomes  0.95  in  contrast  with  the  first  law  efficiency  of  0.98. 
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The  lower  efficiencies  generally  associated  in  the  literature  with 
physical  processes,  such  as  separations  and  compressions,  result  from  using 
an  incremental  efficiency.   For  example,  incremental  efficiencies  for  purge 
acid  gas  removal  and  oil  filtration  are  0.025  and  0.31  respectively.   In 
addition,  these  efficiencies  for  gasifier  gas  compression  and  expansion  are 
0.68  and  0.81  respectively.   The  corresponding  second  law  efficiencies  are 
much  higher.   In  a  process  such  as  fuel  gas  compression  the  work  equivalent 
of  the  gas  is  much  larger  than  the  amount  of  work  lost.  In  other  words,  the 
difference  between  the  ideal  and  the  actual  work  input  is  very  small  relative 
to  the  total  work  equivalent  of  the  fuel  gas .   For  a  gas  with  a  very  low 
work  equivalent,  such  as  air,  the  second  law  efficiency  for  compression 
approaches  the  corresponding  incremental  efficiency.   Consider  for  instance 
that  the  chemical  work  equivalents  of  the  fuel  gas  being  compressed  in  the 
gasifier  gas  compression  unit  are  equal  to  zero.   The  second  law  efficiency 
for  this  compression  process  then  becomes  0.69  which  is  very  close  to  the 
incremental  value  of  0.68  given  above.   This  example  tends  to  show  that  for 
physical  processes  involving  lew  chemical  energy  streams  the  second  law 
efficiency  approaches  the  incremental  efficiency. 

The  second  law  efficiency  as  applied  up  to  this  point  charges  the  overall 
COED  process  with  only  those  inefficiencies  directly  created  by  the  process. 
However,  the  process  uses  items  such  as  oxygen  and  electricity  with  which 
there  is  an  associated  inefficiency.   Electricity,  for  example,  is  generated 
with  an  efficiency  of  approximately  0.38.   This  value  may  be  included  in  the 
work  tern  to  give  an  overall  work  equivalent  for  shaft  work.   An  oxygen  plant 
may  be  added  to  the  overall  process  to  supply  the  necessary  oxygen  at  an 
efficiency  of  about  0.15  (4).   With  these  changes  the  second  law  effi- 
ciency of  the  coal  conversion  endeavor  decreases  from  0.75  to  0.68. 
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In  the  commercial  COED  design  several  items  were  disregarded  that  may 
cause  the  second  law  efficiency  found  here  to  he  too  high.   First,  no  energy 
charge  was  made  for  new  catalyst,  regeneration  of  old  catalyst,  or  general 
chemicals.   The  catalyst  use  amounted  to  approximately  73  kg/hr  (160  lh/hr) . 
An  approximate  energy  charge  for  the  catalyst  and  chemicals  will  be  given 
in  the  next  chapter.   Second,  the  energy  charge  for  waste  treatment  was 
much  less  than  expected  because  of  a  simplified  approach  to  the  difficult 
water  treatment  problem. 
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CHAPTER  7 
EQUIPMENT  AND  CONSTRUCTION  ENERGY  COST 

7.1  Introduction 

When  considering  the  total  energy  input  to  a  process,  the  energy  needed 
to  produce  the  equipment  and  to  construct  the  plant  must  be  included.   The 
energy  input-output  matrix  as  developed  by  Herendeen  (23)  allows  this  addi- 
tional energy  to  be  estimated.   This  type  of  energy  evaluation  uses  Input- 
Output  analysis,  a  technique  borrowed  from  economics,  and  takes  advantage 
of  the  large  data  base  of  intersector  sales  available  from  the  U.S.  Department 
of  Commerce.   The  data  were  last  compiled  in  1963  for  367  sectors,  each 
sector  being  a  segment  of  the  economy  which  produces  goods  or  services, 
e.g.  dairy  farm  products,  machine  shop  products,  and  research  and  development. 
Also  this  energy  evaluation  accounts  for  all  the  steps  in  the  complex 
manufacturer-sales  chain..  For  example,  Herendeen  gives  results  which  indi- 
cate that  the  auto  manufacturer  alone  uses  only  6%  of  the  total  energy  neces- 
sary to  produce  and  to  market  an  automobile. 

The  primary  results  of  the  energy  input-output  matrix  of  interest  here 
are  the  sector  energy  coefficients.   The  two  most  applicable  are  the  energy 
supplied  directly  to  a  sector  per  dollar  of  that  sector's  output,  D  ,  and 
the  total  energy  required  per  dollar  of  market  value,  I  ,  of  which  Dc  is  a 
part.   The  I  and  D  values  for  pertinent  sectors  reproduced  from  Herendeen 
are  listed  in  Table  9.   To  correctly  use  these  sector  energy  coefficients 
the  problem  of  tine  variance  must  be  considered.   The  sector  energy  coef- 
ficients were  developed  in  1973  using  1963  sales  data.   Therefore,  to  use 
these  coefficients  they  must  be  taken  as  time  invariant  in  the  sen3e  that 
only  a  cost  index  is  needed  to  change  their  time  basis.   However,  new  manu- 
facturing techniques  serve  to  reduce  the  sector  energy  coefficients  by 
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TABLE  9.   T  and  D  Sector  Energy  Coefficients 

e     c       * 

for  Various  Sectors 


Sector 


(1)  fabricated  plate  work 

(2)  fabricated  metal  products 

(3)  pumps  and  compressors 

(4)  fabricated  structural  steel 

(5)  ready-mixed  concrete 

(6)  asbestos  products 

(7)  instruments  and  controls 

(8)  pipe,  valves,  and  pipe 
fittings 

(9)  new  construction,  16.7  1.74 
nonresidential  buildings 

(10)  new  construction,  24.8  6.74 
highways 

(11)  electrical  equipment  17.9  0.81 

(12)  clay  refractories  40.3  26.94 

(13)  industrial  inorganic  and  81.4  45.30 
organic  chemicals 


103kcal/3( 

1963) 

103kcal/$(1963) 

29.1 

2.12 

23.2 

1.00 

14.7 

1.90 

31.1 

1.27 

36.1 

2.49 

29.4 

3.46 

10.3 

0.94 

19.5 

2.30 

*Taken  from  reference  (23) 
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reducing  energy  consumption  per  product  unit.   These  new  manufacturing  tech- 
niques, however,  can  be  shown  to  have  little  affect  on  T  with  a  larger 

affect  on  D  .   Remember  that  D  is  a  part  of  T  and  from  Table  9  it  is  seen 
c  c  c 

that  D  is  5-10%  of  T  for  most  sectors.   Also  note  that  D  can  be  thought 
c  c  c  ■ 

of  as  the  energy  it  takes  per  dollar  of  product  to  fabricate  a  product  from 

raw  materials  whereas  the  difference  T  -  D  can  be  considered  the  energy 

consumption  per  dollar  of  product  involved  with  producing  these  raw  materials, 

e.g.  mining  iron  ore,  transporting  iron  ore,  and  processing  iron  ore  into 

steel.   Intuitively  the  term  T  -  D  will  most  likely  be.  time  invariant  with 

increases  in  processing  efficiencies,  less  energy  per  dollar,  being  negated 

by  an  increasing  difficulty  in  mining  ore,  more  energy  per  dollar.   Therefore, 

T  will  probably  be  affected  little.   It  should  be  noted  that  T  obscures  D 
c      r      '  c  c 

in  the  combined  energy  coefficient  calculations  which  lead  directly  to  the 
equipment  and  construction  energy  cost  while  D  will  only  become  important 
when  an  energy  credit  for  scrap  is  determined  for  the  steel  sector. 
7.2  Computational  approach  and  application  to  the  COED  plant 

Knowing  the  fixed  capital  investment  for  the  COED  project,  an  appro- 
priate breakdown  of  its  components  as  prescribed  by  Peters  and  Tinmerhaus 
(24),  and  the  energy  coefficients  D  and  T  ,  an  energy  cost  may  be  assigned 
to  all  of  the  equipment  and  materials  needed  for  the  project's  construction. 
An  equation  will  now  be  briefly  derived  which  will  calcualte  the  total  plant 
energy  cost.   First,  consider  an  expression  for  part  of  the  fixed  capical 
investment. 

S  "  Cl  +  C2  +  C3  +  C4  +  C5  (7'« 


where 


C  =  total  direct  plant  cost  less  land 

C-  =  purchased  equipment  cost 

C9  =  purchased  equipment  installation  cost 
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C,  =  instrumentation  and  controls  cost 
C4  =  piping  cost 
C-  =  offsite  facilities  cost 
The  individual  C  terms  can  be  expressed  as 

C.  =  f±CD  1  -  1,  2,  3,  4,  5  '7.2) 

where 

f .  =  a  fractional  multiplying  factor 

The  energy  associated  with  each  direct  cost  is  determined  from 

E±  =  gtC.  1-1,  2,  3,  4,  5  (7.3) 

where 

E.  =  energy  associated  with  direct  cost  i 

g.  =  combined  energy  coefficient  in  kcal/$  for  direct  cost  i 

determined  from  T  !s  for  appropriate  sectors 

Combining  Eqs .  (7.2)  and  (7.3)  yields 

ET  =  CD  I  g.f,  (7.4) 

1=1 

where 

E_,  =  total  equipment  and  construction  energy  cost 

Eq.  (7.4)  is  used  to  calculate  the  total  plant  energy  cost.   The  multi- 
plying factors,  f.'s,  are  determined  from  percentages  listed  for  a  solid- 
liquid-processing  plant  (24).   These  percentages  were  adjusted  to  fulfill 
the  requirement  given  in  the  COED  economic  analysis  ( 8 )  that  the  offsite 
facilities'  cost  is  30%  of  the  other  direct  costs.   Table  10  presents  these 
multiplying  factors.   The  combined  energy  coefficients,  g.'s,  are  determined 
from  a  combination  of  appropriate  T  sector  energy  coefficients.   As  an 
example  consider  g, ,  the  combined  energy  coefficient  for  purchased  equip- 
ment.  The  major  components  of  this  include  vessels  (fluidized  beds,  tanks, 
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CABLE  10.  Multiplying  Factors  and  Combined  Energy 
Coefficients  for  the  COED  Plant 


Cost  Component 

f  . 
1 

10 

3k 

** 

8i 
cal/$(1963) 

purchased  equipment 

0.42 

22.3 

equipment  installation 

0.16 

24.1 

instrumentation  and 

controls 
piping 

0.06 
0.13 

5.4 
11.9 

offsite  facilities 

0.23 

19.2 

*  as  given  by  Peters  and  Timmerhaus 

**  calculated  from  total  sector  energy  coefficients,  T  's, 
for  appropriate  sectors  c 


53 

and  towers),  heat  exchangers,  cyclones,  and  pumps  and  compressors.   The 
combined  energy  coefficient  is  determined  from  an  average  of  Tc's  in  sectors 
(1),  (2),  and  (3)  as  listed  in  Table  9.   Hie  other  g.'s  are  similarly  found, 
however,  for  these  the  labor  cost  is  significant  and  must  be  included  as 
a  zero  term  in  the  weighting  of  the  appropriate  sector  energy  coefficients. 
To  illustrate  this  additional  point  look  at  the  numerical  calculation  of  g^, 
the  combined  energy  coefficient  for  piping.   The  major  nonlabor  components 
in  this  area  with  corresponding  sectors  are  pipe,  valves  and  fittings, 
structural  supports,  and  insulation  corresponding  to  sectors  (8),  (4), 
and  (6)  (asbestos  is  used  as  a  close  approximation  to  actual  insulation). 
The  combined  energy  coefficient  is  determined  by  a  weighting  of  each  com- 
ponent:  labor  -  50%,  pipe,  valves  and  fittings  -  30%,  and  structural  sup- 
ports and  insulation  -  10%  each.   Therefore  g4  is  calculated  as 

g4  =  0.5  (0)  +  0.3  (19.5)  +  0.1  (31.1)  +  0.1  (29.4) 

g4  -   11.9   x   103   kcal/S(1963) 
The  combined  energy   coefficients   are   listed   in  Table   10.      Finally    the   total 
COED  plant  energy   cost   can  be  determined.      The    total  direct   COED  plant   cost, 
C    ,   was   given  as   505.7  MM$(1975).      With    this    information  plus    the  cost 
index  ratio,    1. 75$(1975) /$(1963) ,    Eq .    (7.4)    may  be  used   to   calculate  a 
total  plant   energy   cost  of   5.85  x  10       kcal.      Thus   a  combined  energy   coet- 
ficient   for   this   plant  is    1.14   x  10      kcal/$(1975) .      This  value  can  be   con- 
sidered a  rough   first   approximation   for  a  solid-liquid-processing  plant. 
With  a  service   life  of  9  years,    as   assumed  by    the  COED  designers    ( 8 ) ,    and 
an  operation  of    330   days/year   the   energy   input  assigned   to    the  equipment   is 
found   to  be  8.21  x   10'    kcal/hr.      This   quantity   is   approximately   1.2%  of   the 
total  energy   input    to    the  COED  plant. 

At   this   point    tne  catalyst  and  chemical   energy   charge  discussed  in 
section  6.3  may  be   estimated  using   total  sector   energy   coefficients,    I 'a. 
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The  COED  economic  analysis  gives  a  combined  catalyst  and  chemical  cost  of 
40.5  MM$(1975).   A  total  sector  energy  coefficient  must  now  be  determined. 
The  catalyst  support  material  is  assumed  tc  be  the  primary  source  of  the 
catalyst's  total  energy  charge.   Sector  (12),  clay  refractories,  seems  to 
be  the  most  representative  sector  corresponding  to  the  support  material. 
The  chemicals  are  assumed  to  correspond  to  sector  (13) ,  industrial  inorganics 
and  organic  chemicals.   The  total  sector  energy  coefficient  for  catalyst 
and  chemicals,  taken  to  be  an  average  of  the  two  separate  coefficients,  is 
60.8  x  10  kcal/S(1963) .   After  the  cost  index  ratio  is  applied  a  final 
energy  charge  for  catalyst  and  chemicals  is  calculated  to  be  1.40  x  10""  kcal. 

Applying  the  operating  time  conversion  the  energy  input  assigned  to  the 

7 
catalyst  and  chemicals  is  found  to  be  1.96  x  10  kcal/hr.   This  non- 
recoverable  energy  will  be  included  in  the  equipment  and  construction 
energy  input  to  give  a  total  energy  input  of  1.02  x  10  kcal/hr  or  approxi- 
mately 1.5%  of  the  total  energy  input  to  the  COED  plant. 
7.3  Salvage  value  considerations 

At  the  end  of  the  plant  equipment's  service  life  an  energy  credit  from 
the  scrap  steel  is  obtained.   Since  the  actual  amount  of  steel  is  unknown 
an  approximation  method  is  used  to  determine  this  credit.   First  look  at 
the  steel  sector  and  consider  what  energy  is  saved  by  the  use  cf  scrap 
steel  over  raw  ores.   This  savings  may  approximately  be  expressed  per  dollar 
of  finished  steel  as 

E   =  T       -  (1-h)  D  (7.5) 

sc     C     -  C     - 

steel  steel 

where 

E   =  energy  credit  coefficient  of  scrap  steel 
sc  r 

h   -   fraction  of  energy  saved  in  the  steel  sector  alone  by 
using  scrap  instead  of  ores 
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In  other  words,  on  a  dollar  of  finished  steel  basis  Eq.  (7.5)  means  that 
the  energy  value  of  scrap  equals  the  total  energy  required  to  produce  steel 
(energy  consumption  from  mining,  transportation,  and  processing)  minus  the 

net  energy  used  only  by  the  steel  sector.  Gyftopolous  (4  )  gives  h  as 

4 
about  0.33  and  Herendeen  gives  Tc      and  Dc      as  6.61  x  10   and  4.13  x 

steel       steel  , 

104  kcal/$(1963)  respectively.   Thus  by  Eq.  (7.5)  Esc  becomes  3.84  x  10 

kcal/$(1963) .  Now  it  is  convenient  to  express  the  energy  savings  when 

scrap  is  employed  as  a  fraction. 

f   =  E  /T  (7-6) 

sc    sc   c  „   , 
steel 

where 

f   =  energy  credit  fraction  for  scrap  steel 
sc       aJ 

This  fraction,  f   ,  is  the  fraction  of  the  energy  charged  to  steel  which  is 
recovered  when  salvage  is  used.   Numerically  fgc  is  approximately  0.58  which 
is  larger  than  h  (0.33),  the  fraction  of  energy  recovered  in  the  steel  sec- 
tor alone  when  salvage  is  used,  because  the  use  of  scrap  steel  saves  all 
the  energy  expended  to  mine  and  process  the  ores.   Next,  this  energy  credit 
fraction  for  scrap  steel  is  used  to  determine  the  combined  energy  credit 

coefficients,  g   's  analogous  to  g.'s.   To  accomplish  this,  first  a  new 

i 
sector  energy  coefficient,  the  total  sector  energy  credit  coefficient,  I  , 

is  determined.   Here  it  may  be  helpful  to  review  the  analogy  of  calculating 

combined  energy  credit  coefficients  so  that  the  equipment  energy  cost  could 

be  found.   In  that  instance  the  total  sector  energy  coefficients,  T  '»i  were 

weighted  to  find  combined  energy  coefficients,  g.'s,  which  were  then  used 

with  the  proper  multiplying  factor,  f.,  and  direct  cost,  C  ,  to  ultimately 

obtain  E  ,  the  total  plant  energy  cost.   Returning  to  the  total  sector  energy 
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credit  coefficient,  T   is  approximated  by  assuming  that  steel  is  the  only 

raw  material  and  it  is  calculated  from 

T   =  (T  -  D  )f  (7.7) 

cc     c    c  sc 

In  other  words  the  total  sector  energy  credit  coefficient  equals  the  part 

of  the  energy  consumption  per  dollar  of  product  involved  with  producing  raw 

materials  that  is  recovered  when  salvage  is  used.   Of  course  this  analysis 

applies  only  to  those  sectors  having  steel  to  salvage,  otherwise  T   equals 

zero.   The  combined  energv  credit  coefficients,  g   's,  are  then  determined 
°'  c, 

1 

analogous  to  the  combined  energy  coefficients  in  the  previous  section.   These 
are  listed  in  Table  11.   Since  the  f .  terms  do  not  change  Eq.  (7.4)  may  be 

applied  after  using  the  cost  index  ratio  to  obtain  an  energy  credit  for  the 

12 
COED  plant  of  2.3  x  10   kcal  corresponding  to  an  energy  credit  flow  of 

3.2  x  10  kcal/hr  or  39%  of  the  total  input  equipment  and  construction  energy 

flow.   This  value  can  be  considered  a  maximum  because  factors  such  as  the 

energy  needed  to  dismantle  the  plant  for  scrap  and  transport  the  scrap  have 

not  been  taken  into  account.   Also,  not  all  of  the  steel  may  be  salvageable. 

Using  this  energy  credit  value  the  net  energy  flow  into  the  plant  is  about 

7.0  x  10  kcal/hr  or  about  1.1%  of  the  total  energy  input,  less  equipment, 

to  the  plant. 


57 


TABLE  11.   Combined  Energy  Credit  Coefficients 
for  the  COED  Plcnt 


3     ■■• 

Cost  Component  10  kcal/$(1963) 

purchased  equipment  12.2 

equipment  installation  4,4 

instrumentation  and 

controls  0 

piping  4.8 

offsite  facilities  6.9 
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CHAPTER  8 
CONCLUSION 

The  second  law  thermodynamic  analysis  of  the  Char  Oil  Energy  Develop- 
ment process  yielded  an  efficiency  of  0.75  when  applied  to  the  immediate 
process.   The  scope  of  the  analysis  was  then  broadened  to  encompass  indirect 
inefficiencies  in  electrical  generation  and  oxygen  production  and  to  include 
the  equipment  and  construction  energy  costs.   This  led  to  a  final  overall 
second  law  efficiency  of  0.67.  This  shows  that  the  second  law  efficiency 
decreases  as  the  scope  of  the  analysis  broadens.   The  analysis  can  be  ex- 
tended to  the  point  where  other  energy  inputs  such  as  coal  mining  and  ash 
disposal  are  included  but  that  was  considered  to  be  remote  and  too  detailed 
for  this  second  law  analysis. 

The  COED  units  involving  physical  processes  had  second  law  efficiencies 
that  were  higher  than  initially  expected.   The  incremental  efficiency  was 
found  to  give  lower,  more  representative  values  for  the  physical  process 
units.   However,  for  units  involving  chemical  transformations  the  second 
law  efficiency  was  deemed  to  be  the  more  useful  of  the  two. 

An  analysis  was  developed  for  determining  energy  charges  for  equipment 
and  construction.   This  charge  was  found  to  be  only  a  small  part,  approxi- 
mately 1%,  of  the  total  energy  involved  with  the  COED  process.   However,  this 
energy  may  be  significant  when  calculated  for  a  process  with  few  high  energy 
chemical  streams.   Also  the  equipment  and  construction  energy  charge  analy- 
sis worked  well  in  the  estimation  of  an  energy  charge  for  indirect  materials 
consumed,  specifically  for  catalyst  and  chemicals  for  which  only  an  economic 

cost  was  available.   The  other  result  of  this  analysis  was  the  approximate 

4 
1.14  x  10  kcal/$(1975)  value  of  the  combined  energy  coefficient  which  may 

prove  useful  for  estimating  the  energy  cost  for  a  solid-liquid-processing 

plant. 


PART  III 


OPTIMIZATION  OF  FLUE  GAS  HEAT  RECOVERY 
FROM  A  SECOND  LAW  PERSPECTIVE 
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CHAPTER  9 
THERMODYNAMIC  AND  ECONOMIC  OPTIMIZATIONS  APPROACHES 

9.1  Heat  recovery  process  description 

Part  III  is  a  further  application  of  the  thermodynamic  analysis  pre- 
sented in  Part  I.   However,  in  this  case  it  will  be  used  to  optimize  the 
design  of  a  process  from  the  standpoint  of  an  overall  efficiency.   Specifi- 
cally, a  flue  gas  heat  recovery  process  design  will  be  optimized  from  a 
second  law  perspective.   Included  in  this  thermodynamic  optimization  will 
be  the  concept  of  an  energy  debit  for  equipment  and  construction  as  dis- 
cussed in  Chapter  7.   The  economic  optimum  design  will  also  be  determined 
for  use  as  a  comparison.   For  both  the  thermodynamic  and  the  economic  opti- 
mizations there  exist  two  system  definitions  that  may  be  employed  when  the 
respective  objective  functions  are  determined.   This  will  be  discussed  in 
section  9.2. 

The  flue  gas  heat  recovery  system  which  is  used  here  Is  similar  to  one 
given  by  James  and  Stokes  (25)  for  the  heat  recovery  from  reformer  furnace 
stack  gases  in  ammonia  plants.   Figure  6  gives  a  block  flow  diagram  of  the 
flue  gas  heat  recover)'  system  plus  furnace.   The  heat  recovery  section 
serves  as  a  combustion  air  preheater  in  which  heat  exchange  is  accomplished 
through  a  number  of  parallel  counter-current  shell  and  tube  heat  exchangers. 
The  flue  gas  and  air  streams  were  given  by  James  and  Stokes  to  enter  the 
heat  recovery  section  at  500°F  (533K)  and  100°F  (311K)  respectively.   The 
pressure  drops  associated  with  the  flow  of  these  streams  through  the  flue 
gas  heat  recovery  system  are  compensated  for  by  forced  and  induced  draft  fans. 

9.2  Determination  of  three  objective  functions 

For  the  previously  described  flue  gas  heat  recovery  system  two  system 
definitions  may  be  employed  that  yield  two  thermodynamic  and  two  economic 
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objective  functions  each  of  which  may  be  optimized.   The  two  system  defi- 
nitions differ  in  that  the  reformer  furnace  unit  is  included  in  one  system 
definition,  the  combined  heat  recovery  system,  while  it  is  not  included  in 
the  other  system  definition,  the  heat  recovery  addition  system.   The  thermo- 
dynamic and  economic  objective  functions  will  now  be  determined  for  each  cf 
these  system  definitions. 

First,  consider  the  combined  heat  recovery  system  as  shown  in  Figure  7 
with  only  inputs  and  outputs.   A  thermodynamic  objective  function  for  this 
system  is  ultimately  obtained  from  the  application  of  Eq .  (3.7)  which  yields 
an  overall  second  law  efficiency,  n  . 

£4  +  V7° 

^1  = (9.1) 

e-   +   £„  +  ex?  +  e    ,   +  W. 
1    2    "6     equip    m 

Next,  let  each  work  equivalent  term,  e  ,  equal  the  sum  of  the  work  equivalent 

i 
of  the  furnace  unit  before  flue  gas  heat  recovery,  e  ,  and  the  work  equiva- 

lent  of  the  flue  gas  heat  recovery  addition,  Ae  .   Therefore  Eq.  (10.1) 

n 

becomes 

i+° 

n,-— i ; i 1 i —  (9.2) 

e,    +  0  +  £,  +  Ae,  +  £        .      +  Ae         .      +  W.      +  AW. 
1  2  2  equip  equip  xn  m 

Since  all  but  the  delta  quantities   are  constant   this    equation  may  be  written 

as 


Cl 


C„  +  Ae„  +  Ae        .      +  AW. 
Z  2  equip  in 


(9.3) 


where 

T 

C-    ■   constant   1      (e.) 
1  4 

Til  T 

C0  =   constant   2      (e,    +  cn   +  s        .      +  W.    ) 
I  12  equip  m 

Therefore,    to  maximize   the  second  law  efficiency   the  group   of    terms 

(Ae„  +  AeeauiD  +  AHj    )    must  be  minimized.      Thus    that  group   of   terms 
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Figure  7.   Overall  Block  Diagram  of  the 
Combined  Heat  Recovery  System 
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comprise  the  thermodynamic  objective  function  for  the  combined  heat  recovery 

system.   The  thermodynamic  optimum  of  this  system  is  given  in  terms  of  the 

i      i         t 

variable  values  that  minimize  (&€-  +  Le  +  AW,  )  instead  of  the  actual 

2     equip     in 

overall  second  law  efficiency,  r\.,    because  the  constant  terms,  C.  and  C~, 
are  unknown.   This  group  of  terms  may  have  a  negative  value  because  as 
more  heat  is  exchanged  less  fuel  is  used  and  Ae  becomes  negative.   It  may 
be  noted  that  the  work  equivalent  of  the  fuel  oil  in  this  objective  function, 
Aejj  is  approximated  by  its  higher  heat  of  combustion  as  has  been  shown  in 
section  4.3.   The  result  of  this  approximation  is  a  thermodynamic  analysis 
that  in  one  respect  closely  resembles  a  first  law  analysis.   Every  thermal 
kcal  recovered  from  the  stack  gas  essentially  reduces  the  fuel  demand  by  a 
corresponding  kcal.   In  other  words  each  kcal  of  thermal  energy  recovered 
from  the  stack  gas  stream  by  the  air  stream  is  roughly  translated  as  in- 
creasing the  work  equivalent  of  the  air  stream  by  that  same  kcal.   However, 
first  law  efficiencies  seldom  take  into  account  the  electrical  generation 
inefficiency  and  the  energy  debit  associated  with  the  equipment.   Corre- 
sponding economic  objective  functions  for  this  system  definition  are  well 
known  and  will  not  be  detailed.   For  example,  James  and  Stokes  (25)  present 
an  economic  analysis  of  a  similar  flue  gas  heat  recovery  system  in  which  they 
give  payout  periods  of  their  system  for  two  purchased  fuel  oil  costs.   The 
economic  objective  function  used  here  is  the  net  present  value  because  it 
gives  a  means  of  direct  comparison  of  alternatives  while  taking  into  account 
the  time  value  of  money.   The  net  present  value  is  optimized  with  respect 
to  fuel  savings,  equipment  costs,  power  costs,  and  various  indirect  costs. 
Second,  consider  the  other  system  definition,  the  heat  recovery  addi- 
tion system  as  given  in  Figure.  3.   A  thermodynamic  objective  function  for 

this  system  is  obtained  directly  from  application  of  the  second  law  effi- 
ciency to  give 
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n,  = *= ' (9.4) 

^     E.  +  Z?°    +    Z  .        +   W. 

5    x6     equip    in 
For  this  heat  recovery  addition  system  the  thermodynamic  optimum  is  deter- 
mined by  the  actual  maximum  overall  second  law  efficiency  value  corresponding 
to  the  optimal  variable  values.   In  contrast  with  the  other  second  law 
efficiency,  n, ,  the  work  equivalent  of  the  preheated  air  stream,  E  ,  does 
not  approximately  equal  the  thermal  energy  recovered  from  the  stack  gas. 
Thus  in  this  system  less  emphasis  is  placed  upon  the  amount  of  heat  re- 
covered since  one  kcal  of  thermal  energy  recovered  equals  only  about  one 
quarter  of  one  kcal  of  work  equivalent.   The  corresponding  economic  objec- 
tive function  for  this  system  is  similar  to  the  previous  one.  However,  in 
actual  practice  this  is  not  as  useful  since  it  is  not  clear  as  to  what 
credit  is  given  to  the  recovered  heat.   Thus  the  economic  objective  function 
for  the  heat  recovery  addition  system  will  not  be  considered  further. 
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CHAPTER  10 
THERMODYNAMIC  AND  ECONOMIC  OPTIMAL  HEAT  RECOVERY  DESIGNS 

10.1  Computational  Scheme 

To  optimize  the  thermodynamic  and  economic  objective  functions  first 
the  physical  conditions,  constraints,  and  assumptions  must  be  established. 
The  inlet  flue  gas  and  air  stream  temperatures  have  been  given  in  section 
9.1.   The  flow  rates  of  these  streams  are  determined  from  the  specific 
heat  recovery  situation  given  by  James  and  Stokes  (25)  where  exit  tempera- 
tures of  300°F  (422K)  and  350°F  (450K)  corresponding  to  exit  flue  gas  and 
air  streams  respectively  produce  a  heat  recovery  of  40  MM  BTU/hr.   Using 
this  information  plus  Perry's  (16)  empirical  formula  for  sulfurless  No.  2 
fuel  oil,  C,   H,  a  combustion  reaction  for  that  fuel  with  10%  excess  air, 
and  constant  pressure  heat  capacities  (13)  an  energy  balance  is  made  to 
determine  the  flow  rates.   They  are  found  to  be  1.124  x  10  kmoi/hr  and 
1.00  x  10  kmol/hr  for  the  flue  gas  and  air  streams  respectively.   It 
should  be  noted  that  these  flow  rates  are  considered  constant  for  this 
analysis  even  though  more  efficient  heat  recovery  will  result  in  the  use  of 
less  fuel  thereby  decreasing  the  air  demand.   However,  calculations  show 
this  to  be  a  good  assumption:   if  the  heat  recovery  doubles  then  the  flow 
rates  are  only  reduced  by  about  0.5%. 

Other  physical  conditions  worth  noting  involve  the  heat  exchanger  place- 
ment, the  heat  exchange  fluid  flow,  and  the  routing  of  specific  fluids  inside 
the  heat  exchangers.  The  heat  exchanger  system  as  shown  previously  in  Figure 
5  consists  of  n  exchangers  connected  in  parallel.   The  parallel  placement 


is  used  instead  of  a  series  placement  because  the  large  gas  flow  rates  of 
about  100  m  /s  would  create  an  unreasonable  overall  pressure  drop  in  the 
series  placement  that  would  cause  the  power  terms  to  completely  overshadow 
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the  other  terms  in  the  thermodynamic  and  economic  optimizations.   Another 
important  physical  condition  is  the  heat  exchange  fluid  flow.   The  heat 
exchange  may  take  place  with  either  parallel  or  counter  flow  of  fluids. 
For  a  given  heat  exchange  area  and  set  of  stream  temperatures  the  loss  of 
work  equivalent  in  a  counter-current  heat  exchanger  is  smaller  than  that 
in  a  co-current  heat  exchanger.   A  co-current  heat  exchanger  is  an  inher- 
ently irreversible  process  while  for  counter-current  heat  exchange  the  loss 
of  work  equivalent  can  be  reduced  by  increasing  heat  exchange  area.   For 
instance,  if  a  counter-current  heat  exchanger  is  balanced,  one  in  which 
total  stream  heat  capacities  are  equal,  and  it  has  infinite  area  the  result 
is  a  reversible  heat  exchange  process.   Because  of  these  facts  the  counter 
flow  design  was  chosen.   It  should  be  noted  that  work  equivalent  is  lost 
whenever  a  temperature  difference  exists  between  two  streams  because  a 
Carnot  engine  could  be  operated  across  that  difference.   The  other  important 
physical  conditions  is  the  routing  of  specific  fluids  inside  the  heat  ex- 
changers.  A  choice  must  be  made  as  to  which  fluid  flows  through  the  shell 
and  which  fluid  flows  through  the  tubes.   Peters  and  Timmerhaus  (24)  list 
several  major  factors  involved  in  determining  the  best  fluid  routing.   These 
factors  include  fouling,  corrosion,  pressure  drop,  and  fluid  velocities. 
However,  fouling  and  corrosion  are  assumed  negligible  and  because  of  the 
closeness  of  the  flow  rates  the  pressure  drops  and  fluid  velocities  will 
be  similar.   Since  the  flue  gas  contains  a  condensible  vapor,  water,  it  is 
possibly  easier  to  separate  liquid  from  the  gas  on  the  shell  side  so  the 
flue  gas  is  routed  through  the  shell  and  the  air  through  the  tubes.   This 
fluid  routing,  though,  is  still  somewhat  arbitrary.   Further  physical  details 
and  assumptions  concerning  the  heat  exchanger  design  are  found  in  Appendix  B. 
It  should  be  noted  that  the  most  important  two  remaining  physical  design 
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limitations  are  that  the  tube  length  cannot  exceed  20  ft  and  that  the  tube 
length  may  also  not  be  less  than  the  shell  inside  diameter.   These  will 
be  further  discussed  in  the  next  section. 

An  important  physical  constraint  to  be  considered  when  heat  recovery 
increases  is  the  existence  of  a  pinchpoint.   A  pinchpoint  occurs  when  the 
temperature  of  the  two  streams  approach  each  other.   When  the  heating  capaci- 
ties of  the  fluids  are  constant  the  pinchpoint  will  occur  at  one  end  of 
the  heat  exchanger,  but  when  the  heat  capacity  of  a  fluid  stream  changes  it 
is  possible  for  the  pinchpoint  to  occur  inside  the  heat  exchanger.   If 
condensation  or  vaporization  of  a  component  is  possible  the  effective  heat 
capacity  of  that  stream  will  increase  and  thus  lead  to  the  possibility  of  an 
interior  pinchpoint.   This  concept  is  well  illustrated  by  Tucker  and  Chen  (26) 
who  cite  the  example  of  flue  gas  heating  of  boiler  feed  water.   In  the  pre- 
sent study  because  of  the  existence  of  a  condensible  specie,  i.e.  water 
vapor,  in  the  flue  gas  the  possibility  of  this  occurence  was  examined  but 
was  found  to  be  nonexistent. 

At  the  same  time  that  the  physical  conditions,  constraints,  and  assump- 
tions are  being  specified  the  independent  design  variables  must  be  determined. 
For  this  flue  gas  heat  recovery  design  these  variables  are  the  following: 
number  of  heat  exchangers  connected  in  parallel,  tube  outside  diameter 
(5/8",  3/4",  or  1"),  exit  temperature  of  the  preheated  air,  and  shell  in- 
side diameter.   Further  details  on  the  ranges  of  these  independent  design 
variables  will  be  given  shortly. 

Once  the  physical  conditions,  constraints,  assumptions,  and  independent 
design  variables  are  determined  a  computer  program  can  be  written  to  find 
the  optimum  of  each  objective  function  given  in  section  9.2.   The  bulk  of 
the  optimization  procedure  for  the  flue  gas  heat  recovery  design  remains 
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Che  same  irrespective  of  whether  the  objective  function  is  economic  or 
thermodynamic.   For  example,  the  heat  exchanger  design  section  remains  un- 
changed for  any  of  the  objective  functions.   Briefly  the  computer  program 
written  for  this  optimization  problem  conducts  an  iterative  search  whereby 
the  independent  design  variables  are  incremented  one  at  a  time  over  a  pre- 
determined search  interval.   Specifically,  the  computer  program  is  designed 
to  incrementally  search  over  a  range  of  shell  inside  diameters  and  exit 
preheated  air  temperatures  before  the  tube  outside  diameter  and  the  number 
of  heat  exchangers  is  incremented.   This  results  in  a  table  of  exit  pre- 
heated air  temperatures  versus  shell  inside  diameters  for  a  particular 
value  of  tube  outside  diameter  and  number  of  heat  exchangers.   The  basic 
search  ranges  for  this  table  are  the  following:   shell  inside  diameter  — 
50-300  cm  (19.7  -  118  in)  and  exit  preheated  air  temperature  (stream  3  in 
Figure  5)  —  330-530K  (134.3  -  494. 3°F) .   The  program  computes  values  of 
the  objective  function  for  all  the  design  variable  combinations  and  finds 
the  maximum  objective  function  value  along  with  corresponding  values  of 
the  design  variables.  A  flowchart  of  this  computer  program  is  presented  in 
Figure  9.   Appendix  B  lists  the  program  and  sample  search  results. 


(         START    J 


Declare  and 

Initialize  Variables 

and  Arrays 


Input  data:      shell-   and 
1  tube-side  fluid  proper- 
ties, independent   design 
variable  values,    plus 
indicators   of  which 
optimization  results 
to  be  outputted 


Initialize 

optimization 

loop  variables 


loop  over  range  of  values 
for  number  of  heat  exchangers 
and  tube  outside  diameters 


range  of  temperature  and 
shell  inside  diameters 


Calculate  shell- 
side  fluid  exit 
temperature  from 
an  energy  balance 


Calculate  shell- 
side  heat  trans- 
fer parameters 
Re,  Pr,  and  h 


Calculate  tube- 
side  heat  trans- 
fer parameters 
Se,  Pr,  and  h 
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Determine   tube 

length   and  number 

of  baffles    from 


AT 


In' 


U   ,    and   A 


Calculate  shell- 
side  pressure 
drop  by  iteration 


Calculate  tube- 
side  pressure 
drop 


Calculate  blower 
and  heat  ex- 
changer total 
costs 


Determine  the  net 

present  value  of 

the  combined  heat 

recovery  system 


Find  the  net  work 
equivalent  for 

the  combined  heat 
recovery  system 


Find  the  overall 
second  law  effi- 
ciency of  the 
heat  recovery 

addition  system 
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'  Output  desired  search 
results:  all  values 
/for  either  net  present 
/value,  net  heat  equiv- 
alent, and/or  the 
overall  efficiency 


I  Output  the  desired  con- 
ditions and  parameters  j 
for  the  maximum  net 
/present  value,  minimum 
net  work  equivalent, 
and/or  maximum  overall/ 
second  law  efficiencv 


STOP 


Figure  9 .   Flowchart  of  the  Flue  Gas  Heat 
Recovery  Optimization  Program 
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10.2  Optimization  results 

The  results  of  the  optimizations  of  the  economic  objective  function 
and  the  two  thermodynamic  objective  functions  determined  in  section  9.2  are 
given  here.   At  the  outset  it  will  be  helpful  to  explain  the  basis  of  com- 
parison for  these  three  objective  functions  and  also  to  define  some  of  the 
terminology  used  in  dealing  with  the  optimums  of  these  objective  functions. 
Since  four  independent  design  variables  are  used  the  number  of  comparison 
methods  and  the  detail  of  each  can  easily  be  seen  to  be  large  and  complex. 
However,  it  may  be  revealed  here  that  only  three  of  these  variables  can  be 
used  for  comparison  since  the  tube  outside  diameter  variable  is  found  for 
all  optimums  to  be  one  inch,  the  largest  value  tested.   For  simplicity  only 
one  independent  design  variable,  the  number  of  heat  exchangers,  is  chosen 
as  the  basis  of  comparison.   That  variable  is  chosen  over  the  others  because 
of  its  inherent  unconstrained  range  of  values  and  its  major  role  in  each 
objective  function.   For  instance,  the  shell  inside  diameter  is  limited  to 
about  118  in.  by  design  constraints  and  the  exit  preheated  air  temperature 
is  limited  to  the  temperature  of  the  inlet  flue  gas  by  an  end  pinchpoint 
condition.   The  number  of  heat  exchangers  is  a  major  cost  factor  in  the 
economic  objective  function  and  because  of  its  determination  of  individual 
exchanger  gas  flow  rates  it  is  a  major  factor  in  determining  pressure  drops. 
It  should  be  remembered  that  with  the  very  large  flow  rates,  about  100  a  /s , 
associated  with  this  system  only  a  small  pressure  drop  produces  both  stag- 
gering economic  and  thermodynamic  energy  costs. 

Terminology  pertaining  to  the  optimums  of  the  three  objective  functions 
is  defined  next  in  light  of  the  previous  discussion.   An  objective  function 
optimum  at  a  specified  number  of  heat  exchangers  will  be  referred  to  as  a 
local  optimum  whereas  the  objective  function  optimum  found  with  respect  to 
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all  the  independent  design  variables  including  the  number  of  exchangers  will 
be  referred  to  as  the  global  optimum.   Economic  objective  function  optimums 
will  be  considered  as  net  present  value  optimums  which  correspond  to  the 
combined  heat  recovery  system  (see  Figure  6) .   Optimums  of  the  two  thermo- 
dynamic objective  functions  will  be  considered  as  exchanger  optimums  and  as 
exchanger-plus-furnace  optimums  corresponding  to  the  heat  exchanger  addi- 
tion system  (see  Figure  6)  and  the  combined  heat  recovery  system  (Figure  6) 
respectively. 

The  simplistic  approach  mentioned  in  the  early  part  of  this  section  is 
well  in  tune  with  the  major  objectives  of  this  heat  recovery  optimizational 
study  which  are  to  find  the  conditions  at  each  global  optimum  and  to  compare 
the  economic  and  thermodynamic  optimums.   As  mentioned,  because  of  this  sim- 
plistic approach  it  is  not  considered  pertinent  to  consider  the  complex 
functionality  that  exists  between  the  heat  exchanger  design  criteria  and 
the  design  variables.   The  major  heat  exchanger  design  criteria  of  each 
global  optimum  are  given  in  Table  12  while  Table  13  gives  the  corresponding 
objective  function  value.   Table  12  gives  outlet  flue  gas  temperatures  that 
indicate  the  possibility  of  water  vapor  condensation.   However,  no  arrange- 
ments exist  in  the  computational  procedure  for  condensation  because  water 
vapor  is  only  1.3  mole  percent  of  the  total  gas  corresponding  to  a  very  low 
dew  point.   In  Table  13,  the  negative  value  of  the  thermodynamic  objective 

function  for  the  combined  heat  recovery  system,  (Ae„  +  Ae      +  AW.  ),  occurs 

2     equip     in 

because  as  pointed  out  in  section  9.2  fuel  is  saved  which  results  in  a  de- 

t 
crease  of  total  fuel  work  equivalent,  i.e.  a  negative  Ae_,  that  in  this  case 

is  larger  in  magnitude  than  the  increase  in  equipment  and  electrical  work 

equivalents.   It  is  interesting  in  the  light  of  Part  II  of  this  thesis  that 
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Table  13.      Optimal  Values    of    the  Economic  and 
Thermodynamic  Objective  Functions 


Objective  Function  Value 

Net  Present  Value*  4.429    x  10      $(1979) 

i  i  7 

(Ae-  +  Ae        .     +  AW.    )*  -1.376  x  10     kcal/hr 

2.  equip  in 

Overall  Second  Law  Efficiency  ru**  0.4704 


*     for  the   combined  heat  recovery   system 
**     for   the  heat   exchanger  addition  system 
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a  first  law  efficiency  corresponding  to  the  heat  recovery  addition  system 
is  calculated  without  the  inclusion  of  equipment  energy  or  electrical  gener- 
ating inefficiencies  to  be  approximately  0.70.   A  corresponding  incremental 
efficiency  is  found  to  be  negative  since  the  work  equivalent  change  between 
output  and  input  streams  is  negative  and  since  work  is  required  which  by 
convention  is  positive.   In  other  words  work  is  required  in  a  process  cap- 
able of  performing  work,  i.e.  a  Carnot  heat  engine  could  be  operated  between 
the  stream  temperatures. 

The  other  purpose  of  this  study,  economic  and  thermodynamic  optimum 
comparisons,  is  in  part  fulfilled  by  Figure  10.   In  this  figure  local  opti- 
mums normalized  with  respect  to  corresponding  global  optimums  are  plotted 
versus  the  number  of  heat  exchangers.   The  normalized  net  present  value 
curve  is  observed  to  peak  rather  sharply  at  8  heat  exchangers  in  contrast  to 
the  two  thermodynamic  curves  peaking  gradually  at  27  and  35  heat  exchangers 
corresponding  to  the  normalized  exchanger  curve  and  the  normalized  exchanger- 
plus-  furnace  curve  respectively.   However,  these  latter  two  curves  reach 
0.95,  essentially  the  maximum,  at  about  12  and  15  heat  exchangers  respectively. 
It  is  observed  that  the  normalized  exchanger  curve  for  small  numbers  of  heat 
exchangers  has  a  lower  slope  than  does  the  normalized  exchanger-plus-fumace 
curve  yet  it  reaches  its  global  optimum  before  the  normalized  exchanger- 
plus-furnace  curve  does.   The  normalized  exchanger  and  exchanger-plus- 
furnace  curves  as  it  should  here  be  remembered  correspond  to  the  heat  ex- 
changer addition  system  and  the  combined  heat  recovery  system  respectively. 
As  discussed  in  section  9.2  in  the  heat  exchanger  addition  system  less 
emphasis  is  placed  upon  the  amount  of  heat  recovered  than  in  the  combined 
heat  recovery  system.   The  reason  for  this  is  that  since  the  furnace  is 
inefficient,  because  of  the  combustion  of  fuel  any  lessening  of  the  fuel 
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requirement  means  that  the  quantity  of  heat  is  substituted  for  a  higher 
grade  source.   Thus  large  initial  increases  in  heat  recovery  have  more 
effect  on  the  combined  heat  recovery  system's  objective  function  resulting 
in  a  larger  initial  slope  for  the  normalized  exchanger-plus-furnace  curve. 
The  reason  that  the  global  optimum  of  the  normalized  exchanger  curve  occurs 
at  a  fewer  number  of  heat  exchangers  than  does  the  global  optimum  of  the 
normalized  exchanger-plus-furnace  can  be  somewhat  similarly  explained. 
First,  consider  what  factors  actually  cause  a  global  thermodynamic  optimum 
to  occur  with  respect  to  the  number  of  heat  exchangers.   It  is  found  that 
as  the  number  of  heat  exchangers  increase  the  energy  debit  associated  with 
the  equipment  (see  Chapter  7)  finally  becomes  a  large  enough  fraction  of 
the  energy  input  so  that  its  incremental  increase  is  greater  than  any 
incremental  decrease  in  energy  cost  due  to  a  decreased  pressure  drop  or 
increased  heat  recovery.   Now,  as  stated  above,  recall  that  the  heat  ex- 
changer addition  system,  associated  with  the  normalized  exchanger  curve, 
places  less  value  on  the  heat  recovered  than  does  the  combined  heat  recovery 
system,  associated  with  the  normalized  exchanger-plus-furnace  curve.   Thus 
the  equipment  energy  debit  is  a  larger  fraction  of  the  energy  input  in  the 
exchanger  case  than  in  the  exchanger-plus-furnace  case.   Therefore,  as 
the  number  of  heat  exchangers  increase  the  global  exchanger  optimum  occurs 
before  the  global  exchanger-plus- furnace  optimum.   It  may  be  noted  that 
the  equipment  energy  fraction  of  the  total  energy  input  is  0.045  for  the 
heat  exchanger  addition  system,  the  exchanger  case.   A  corresponding  value 
for  the  combined  heat  recovery  system,  however,  may  not  be  calculated  because 
as  seen  in  section  9.2  the  equipment,  electrical,  feed,  product  and  work 
equivalents  for  the  furnace  system  alone  are  not  known.   In  contrast  2n 
equipment  cost  fraction  of  the  global  net  present  value  optimum  is  0.48. 
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The  above  discussion  gives  global  optimums  for  the  economic  and  thermo- 
dynamic objective  functions  and  shows  the  functionality  of  these  objective 
functions  with  the  major  independent  design  variable,  the  number  of  heat 
exchangers.   However,  that  discussion  does  not  provide  for  a  direct  com- 
parison between  the  economic  and  thermodynamic  objective  functions.   It 
should  be  remembered  that  in  Figure  10  these  objective  functions  are  only 
compared  on  the  basis  of  the  sasie  number  of  heat  exchangers  and  as  it 
happens  the  same  tube  outside  diameter.   The  other  independent  design 
variables  are  not  fixed.  Thus  upon  further  consideration  a  new  comparison 
approach  can  be  found  that  yields  Figure  11.   In  this  figure  the  economic 
and  thermodynamic  objective  functions  are  directly  compared  for  the  combined 
heat  recovery  system.   Two  economic- thermodynamic  comparison  ratios  are 
plotted  versus  the  number  of  heat  exchangers.   One  ratio  is  referred  to  as 

NPV        /SPV,       and  the  other  is  referred  to  as  ri         /ri 

1l„__»     locopt  'NPV,       'loc  opt 

locopt        c  locopt      " 

To  determine  these  ratios  for  a  specific  number  of  heat  exchangers  first 
the  local  net  present  value  optimum,  :n>viOCODt'  and  t'le  l°cal  exchanger- 
plus- furnace  optimum,  ru      ,  are  found.   Next,  the  values  of  the  inde- 
pendent design  variables  corresponding  to  the  local  exchanger-plus-furnace 

optimum  are  used  to  calculate  a  new  net  present  value,  NPV        ,   Likewise 

loc  opt 
the  values  of  the  independent  design  variables  corresponding  to  the  local 

net  present  value  optimum  are  used  to  calculate  a  new  exchanger-plus-furnace 

value,  n„pv      -  Finally  the  new  values  are  divided  by  the  corresponding 

locopt 
local  optimums  to  arrive  at  the  two  economic- thermodynamic  comparison  ratios, 

NPV        /NPV,     „  and  n„_.„       /n,         It  should  be  remembered  that 
n,  locopt      NPV,        locopt. 

locopt  locopt 

local  optimums  are  used  here  implying  that  a  specific  number  of  heat  ex- 
changers are  used  for  each  ratio  calculation.  Therefore,  these  ratios  make 
possible  the  direct  comparison  between  the  net  present  value  objective 
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function  and  the  exchanger- plus- furnace  objective  function.   Returning  to 
Figure  11  it  is  first  noticed  that  for  one,  two,  and  three  heat  exchangers 
both  ratios  are  unity.   This  implies  that  both  the  economic  and  the  thermo- 
dynamic optimal  designs  are  identical  for  these  three  cases.   One  explana- 
tion for  this  occurrence  is  that  as  had  been  previously  discussed  the 
pressure  drop  is  the  dominant  factor  in  both  objective  functions  for  this 
range  of  heat  exchangers  thereby  producing  identical  designs  that  minimize 

its  affect.   Another  observation  is  that  tU—-,  /n,     „  gradually 

NPV,        loc  opt         ' 
loc  opt       r 

decreases  asymptotically  to  0.74  while  NPV        /NPV,       decreases 

n,  loc  opt 

loc  opt 

steeply  after  about   4   exchangers   and  becomes   negative  after   15   exchangers. 
It  may  be   instructive   to   examine    these  ratios   at   the  economic  and   thermo- 
dynamic global  optimums   of   8  heat   exchangers   and  approximately  15  heat 

exchangers   respectively.      For  8  heat  exchangers   n*rpv  ^locoot    *s 

loc  opt 
about  0.9  while  NPV  /NPV,  .     is   about  0.5.      These  values    imply 

n,        loc  opt 
loc opt 

that  the  use  of  global  economic  optimum  design  conditions  work  well  in 
approximating  the  actual  local  thermodynamic  optimum  but  that  the  reverse 
of  using  local  thermodynamic  optimum  design  conditions  does  not  satisfac- 
torily approximate  the  actual  global  economic  optimum.   The  other  case  being 
considered  that  of  15  heat  exchangers  gives  a  similar  conclusion  since 

n  /n,       is  0.76  and  NPV        /NPV,       is  0.   One  addi- 

NPV,       'loc  opt  n,     ..    loc  opt 

loc opt  loc opt 

tional  point  may  be  brought  out  concerning  the  global  economic  optimum  case. 

Since  the  normalized  exchanger-plus-furnace  curve  in  Figure  10  is  actually 

a  plot  of  n,      /n  ,  ,,   .  versus  the  number  of  heat  exchangers  its 

value  at  8  heat  exchangers,  0.S4,  may  be  multiplied  by  the  Ivroy       / 

loc  opt 

n.       value  of  0.91  to  give  about  0.77  for  the  new  ratio  nvrpV  / 

p  global  opt 

n  ,  ,  ,       This  represents  the  aoproach  of  the  global  economic  optimum 
global  opt.         "  rr  s 

design  to  the  global  thermodynamic  optimum  design. 
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CHAPTER  11 
CONCLUSION 

The  thermodynamic  optimization  of  a  flue  gas  heat  recovery  system  using 
the  principles  of  second  law  analysis  yielded  several  enlightening  conclu- 
sions.  The  thermodynamic  optimums  produced  heat  recovery  system  designs 

that  physically  seemed  impractical  because  they  required  approximately  30 

2 
large  heat  exchangers  of  about  20,000  ft  each.   However,  specific  designs 

for  only  13  heat  exchangers  were  determined  to  be  nearly  as  optimal  as  the 
30  exchanger  cases  since  those  designs  gave  values  of  the  thermodynamic 
objective  functions  within  95%  of  the  optimal  values.   For  comparison  a 
distinct  economic  optimum  occurred  at  8  heat  exchangers. 

In  addition,  the  design  conditions  of  an  economic  optimum  for  a  pre- 
determined number  of  heat  exchangers  was  deemed  to  approximate  within  at 
least  20%  a  corresponding  thermodynamic  optimum  for  the  same  number  of  heat 
exchangers  as  long  as  that  number  was  less  than  about  13.   The  reverse 
approximation  of  using  thermodynamic  optimal  conditions  to  calculate  an 
economic  optimum  for  a  preset  number  of  heat  exchangers  was  not  satisfac- 
tory.  This  points  out  that  the  cost  of  energy  is  now  becoming  high  enough 
so  that  it  is  economical  to  optimize  the  design  of  this  flue  gas  heat  re- 
covery system  within  80%  of  the  actual  thermodynamic  optimum.   However, 
as  this  percentage  slowly  approaches  100%  the  economics  rapidly  become  less 
favorable.   Thus  this  additional  20%  proves  to  still  be  economically  quite 
costly. 
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NOMENCLATURE 

a.,  bi    =  constants  of  a  two  parameter  heat  capacity  equation  for  com; 

ponent  i  used  in  Eq.  (4.12),  kcal/(kmol  K)  and  kcal/ (kmol  K  ), 
respectively 

Cp       =  constant  pressure  heat  capacity  for  component  a  (a  =  i  or  j) , 
a         kcal/ (kmol  K) 

C         =   total  direct  plant  cost  less  land,  $;  defined  in  Eq .  (7.1) 

C.        =   cost  component  of  the  total  direct  cost  (i  =  1,  2,  3,  4,  5), 
1  $;  defined  in  Eq .  (7.1) 

C,,  C2    =  constants;  defined  in  Eq . (9 . 3) 

D_        =  direct  sector  energy  coefficient,  kcal/$ 

E        •  energy  associated  with  direct  cost  i,  kcal;  defined  in  Eq.  (7.3) 

E         =  energy  credit  coefficient  for  scrap  steel,  kcal/$;  defined 
SC  in  Eq.  (7.5) 

E~        =   total  equipment  and  energy  construction  cost,  $;  defined 
in  Eq.  (7.4) 

f.  =  fractional  multiplying  factor;  defined  in  Eq.  (7.2) 

f         =  energy  credit  fraction  for  scrap  steel;  defined  in  Eq.  (7.6) 

g         =  combined  energy  credit  coefficient  for  direct  cost  i  determined 
i  from  T   's  for  appropriate  sectors,  kcal/$ 

g.  =  combined  energy  coefficient  for  direct  cost  i  determined  from 

T  's  for  appropriate  sectors,  kcal/$ 

AG        -  standard  free  energy  of  reaction,  kcal/kmol 

o 
AG-       =  standard  free  energy  of  formation  for  component  a,  kcal/kmol 

h         =  fraction  of  energy  saved  in  the  steel  sector  alone  by  using 
scrap  instead  of  ores 

H  =  enthalpy  of  system,  kcal 

H.  =  partial  molar  enthalpy  of  component  i,  kcal/kmol 

H  =  enthalpy  of  system's  reference  state,  kcal 

H  .  =  partial  molar  reference  enthalpy  of  component  i,  kcal/kmol 

AH  =  total  change  in  enthalpy  of  system,  kcal/hr 
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AH        =  standard  enthalpy  of  formation  for  component  a(a  =  i  or  j), 
kcal/kmol 

AH_  _     =  standard  enthalpy  change  in  reaction,  kcal/hr 

AH        =  enthalpy  change  as  products  in  their  standard  states  are 
taken  from  298K  and  1  atm  to  their  actual  temperatures, 
kcal/hr 

AH„       ■  analogous  to  AH  for  reactants  instead  of  products 

AH.       =  standard  heat  of  combustion  for  component  i,  kcal/kmol 

AH^qo     =  standard  heat  of  combustion  (HO,  J,  kcal/kmol 

"c  c 

AH        =  analogous  to  AH    except  with  kcal/lb  units 

c    c 

AH  ,  AH   =  standard  heat  of  combustion  for  products  and  reactants, 

respectively,  kcal/kmol 

m        =  number  of  product  streams 

n         =  number  of  reactant  streams 

n        =  stoichiometric  coefficient  for  oxygen  in  Eq.  (4.8) 

n.        =  molar  flowrate  of  species  i,  kmol/hr 

NPV-      =  local  net  present  value  optimum,  S 
loc  opt  r     * 

NPV       =  net  present  value  associated  with  independent  design  variable 
loc  opt   values  corresponding  to  the  local  exchanger-plus-furnace  opti- 
mum, $ 

p  =  stoichiometric  coefficient  for  hydrogen  in  Eq.  (4.8) 

p.  =  product  stream  j 

p  =  reference  pressure,  1  atm 

q  =  stoichiometric  coefficient  for  nitrogen  in  Eq.  (4.8) 

Q         =  net  heat  flow  across  system  boundaries  with  heat  flow  in  being 

a  positive  quantity  and  heat  flow  out  being  a  negative  quantity, 
kcal/hr 

Q.  ,Q     =  heat  flows  into  and  out  of  the  system,  respectivelv,  taken 

IS       out  .,.-,/, 

as   positive,    kcal/hr 
R  =     ideal  gas   law   constant,    kcal/(kmol  K) 

r.  =      reactant  stream  i 


S  =      entropy   of  system,    kcal/K 
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S^  =  partial  molar  entropy  of  component  1,  kcal/(kmol  K) 

S  =  entropy  of  system's  reference  state,  kcal/K 

S  .  =  partial  molar  reference  entropy  of  component  i,  kcal/(kmol  K) 

T  =  heat  transfer  temperature,  K 

T  ■  temperature  of  component  a  (a  ■  i  or  j),  K 

T  =  total  sector  energy  coefficient,  kcal/$ 

T         =   total  sector  energy  credit  coefficient,  kcal/$ 
cc 

T        »  reference  temperature,  298  K 
o 

V.        =  molar  volume  of  species  i,  £/kmol 

W        =  net  shaft  work  between  system  and  surroundings  with  the  con- 
vention that  work  done  on  the  system  is  positive  and  work 
done  by  the  system  is  negative,  kcal/hr 

wj  >W  .  "  shaft  work  done  on  the  system  and  by  the  system,  respectively, 
in  out      ,         ■-•11 /i  >r      " 

taken  as  positive,  kcal/hr 

i 

W.        =  shaft  work  into  the  furnace  unit  before  flue  gas  heat  recovery, 

kcal/hr;  see  Eq.  (9.2) 

t 

AW       =  shaft  work  into  the  flue  gas  heat  recovery  addition,  kcal/hr; 

111         see  Eq.  (9.2) 

W         "  actual  work  equivalent  used  by  a  process,  kcal/hr;  defined 
net  in  Eq.  (3.8) 

W        =  maximum  absolute  value  shaft  work  done  by  system  under  rever- 
sible conditions,  kcal/hr 


i 


=  mole  fraction  of  species  i 


a,g,Y     "  stoichiometric  coefficients  for  carbon,  hydrogen,  and  oxygen 
respectively 

Y.        =  quantity  used  in  estimating  the  heat  of  combustion  of  a 
variety  of  organic  material;  defined  in  Eq.  (4.9) 

6  =  stoichiometric  coefficient  for  nitrogen 

€  =      total  work  equivalent  of  a  system,  kcal 

-„u  .     =  work  equivalent  due  to  chemical  energy  at  T  and  p  ,  kcal/hr; 
chem        j  ^   1  .   «    r,    ™%  o     *o         ' 

defined  m  Eq.  (4.20) 
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equip 


work  equivalent  associated  with  equipment  and  construction 
energy  cost,  amoritized,  kcal/hr;  see  Eq .  (9.1) 


equip 


equip 


=  analogous  to 


auip 


for  the  furnace  unit  before  flue  gas  heat 


Ae. 


Ae. 


Ae, 


Ae. 


i 
ein 


Ae 


"phys 

1T 


loc  opt 


recovery,  kcal/r.r;  defined  in  Eq.  (9.2) 

=  analogous  to  e   .   for  the  flue  gas  heat  recovery  addition, 
kcal/hr;  definllPin  Eq.  (9.2) 

=  partial  molar  work  equivalent  of  species  i 

=  change  in  work  equivalent  from  pure  component  i  at  T  and 
p   to  other  conditions,  kcal/kmol;  defined  in  Eq .  (4.19) 

change  in  work  equivalent  due  to  a  change  in  pressure  from 
p  at  T  and  constant  x.  for  component  1,  kcal/kmol;  defined 
in  Eq.  ?4.17)  1 

=  change  in  work  equivalent  due  to  a  change  in  temperature 
from  T  at  constant  p  and  x  for  component  i,  kcal/kmol; 
defined  in  Eq.  (4.12)  for  a  two  parameter  heat  capacity  and 
in  Eq.  (4.13)  for  a  constant  heat  capacity 

=  change  in  work  equivalent  due  to  a  change  in  mole  fraction 
from  x.  =  1  at  T  and  p  ,  kcal/kmol;  defined  in  Eq.  (4.18) 

■  standard  state  work  equivalent  for  species  i,  kcal/kmol 

=  work  equivalent  input  due  to  process  streams,  kcal/hr; 
defined  in  Eq .  (3.7) 

■  work  equivalent  of  any  stream  n  (n  =  1,  2,  ...),  kcal/hr 

»  work  equivalent  of  any  stream  n  (n  =  1,  2,  ...)  associated 
with  the  furnace  unit  before  flue  gas  heat  recovery,  kcal/hr; 
defined  in  Eq .  (9.2) 

=  analogous  to  e   except  associated  with  the  flue  gas  heat 
recovery  system,  kcal/hr;  defined  in  Eq .  (10.2) 

=  work  equivalent  output  due  to  process  streams,  kcal/hr; 
defined  in  Eq .  (3.7) 

■  work  equivalent  due  to  changes  from  pure  components  at  T 
and  p  ,  kcal/hr;  defined  in  Eq.  (4.21) 

=  second  law  efficiency;  defined  in  Eq.  (3.7) 

=  incremental  efficiency;  defined  in  Eq.  (3.8) 

=      local  exchanger-plus-furnace  optimum,    kcal/hr 


'NPV 


=  exchanger-plus- furnace  value  associated  with  independent 
loc  opt   design  variable  values  corresponding  to  the  local  net  present 
value  oDtimum 
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overall  second  law  efficiency  of  the  combined  heat  recovery 
system;  defined  in  Eq .  (9.2) 

■  overall  second  law  efficiency  of  the  heat  exchanger  addition 
system;  defined  in  Eq.  (9.4) 

=  stoichiometric  coefficient  for  sulfur 
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APPENDIX  A 
This  appendix  refers  to  the  COED  process  analyzed  in  Part  II.  A 
detailed  block  flow  diagram  listing  all  the  stream  numbers  is  given  in 
Figure  12.   Corresponding  stream  descriptions  are  presented  in  Table  14. 
Following  that  table  a  listing  of  the  COED  thermodynamic  computer  program 
with  output  is  given.   Throughout  the  program-listing  comments  are  included 
that  describe  the  procedures  used  and  give  many  of  the  minor  assumptions. 
In  reference  to  the  computer  output  the  first  table  lists  the  composition 
and  physical  conditions  of  each  stream.  The  next  table  essentially  gives 
the  stream  information  entailed  In  Figure  12  plus  the  utility  useage  infor- 
mation.  After  that  the  material  balance  check  is  given  which  shows  a  slight 
discrepancy  in  the  balance  that  is  considered  insignificant.   Following 
that  table  the  molar  composition,  physical  conditions,  energy  values,  and 
work  equivalent  values  are  presented  for  each  stream.   It  may  be  helpful 
to  define  the  notation  that  describe  the  last  two  items  mentioned  above. 
ACH  and  ECH  refer  to  the  total  standard  heat  of  combustion  and  total  standard 
state  chemical  work  equivalent  respectively.   APH  refers  to  the  total  physi- 
cal enthalpy  change  from  standard  state  while  EPH  refers  to  the  total  physi- 
cal work  equivalent  change  from  standard  state.   The  E  and  A  represent  the 
total  thermal  energy  and  total  work  equivalent  respectively.   The  next  table 
lists  the  total  input  thermal  energy  and  work  equivalent  for  each  process 
unit.   Finally  the  last  table  gives  the  final  efficiencies  and  energy  balance 
closure. 
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Table  14.   COED  Coramer 

cial  Design  Process 

Stri 

earn  Desc 

riptions 

Stream 

Stream 

Number 

Description 

Number 

Description 

1 

Feed  Coal 

45 

Spent  Catalyst  Plus  Coke 

2 

Coal  Fines 

46 

Waste  Water 

3 

Fuel  Gas 

47 

Syncrude 

4 

Air 

48 

Purge  Gas 

5 

Stack  Gas 

49 

Stripper  Off-Gas 

6 

Crushed  Coal 

50 

Water  to  Scrubber 

7 

Waste  Water 

51 

Acid  Gas 

8 

Fuel  Gas 

52 

Clean  Purge  Gas 

9 

Air 

53 

Hydrogen 

10 

Stack  Gas 

54 

Hydrogen 

11 

Second  Stage  Transport  Gas 

55 

Hydrocarbons 

12 

First  Stage  Char 

56 

Water 

13 

First  Stage  Waste  Liquor 

57 

Air 

14 

Second  Stage  Transport  Gas 

58 

Steam 

15 

Oxygen  to  Stage  4 

59 

Ash 

16 

Steam  to  Stage  4 

60 

Waste  Water 

17 

Second  Stage  Off-Gas 

61 

Gasifier  Product  Gas 

18 

Char 

62 

Waste  Water 

19 

Second  Stage  Transport  Gas 

63 

Gasifier  Product  Gas 

20 

Unfiltered  Oil 

64 

Acid  Gas 

21 

Waste  Water 

65 

Acid  Gas 

22 

Pyrolysis  Product  Gas 

66 

Stack  Gas 

23 

Water  to  Scrubber 

67 

Product  Sulfur 

24 

Water  to  Stripper 

68 

Air 

25 

Recycled  Pyrolysis  Gas 

69 

Fuel  Gas 

26 

Stripper  Off-Gas 

70 

Clean  Gasifier  Product  Gas 

27 

Waste  Water 

71 

Clean  Gasifier  Product  Gas 

28 

Pyrolysis  Product  Gas 

72 

Plant  Fuel  Gas 

29 

Acid  Gas 

73 

Clean  Gasifier  Product  Gas 

30 

Clean  Pyrolysis  Product 

Gas 

74 

Plant  Product  Gas 

31 

Clean  Pyrolysis  Product 

Gas 

75 

Stack  Gas 

32 

Clean  Pyrolysis  Product 

Gas 

33 

Clean  Pyrolysis  Product 

Gas 

34 

Clean  Pyrolysis  Product 

Gas 

35 

Clean  Pyrolysis  Product 

Gas 

36 

Steam 

37 

Water 

38 

Waste  Water 

39 

Hydrogen 

40 

Purge  Gas 

41 

Filter  Cake 

42 

Filteraid 

43 

Hydrogen 

44 

Filtered  Oil 
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c 

C    SECOND    LAW    THEHMOOY NAM  1 C    ANALYSIS    UK     A    GcNlRAL    COMMERCIAL    COAL 
C    CONVERSION    PROOESS — APPlIEJ    S Pt CI F I CA LL Y    TO    FMC    CORP.'S     CHArf    OIL 
C    ENERGY    DEVELOPMENT     PROCESS 


C    WRiTTEN    BY    TERKY    L.     UNRUH    UNDER     Tht    ADVISEMENT    CF     DR.     6.     S.     KYLE, 
C    OEPT.    OH     CHEMICAL    ENGINEERING.     K.ANSAS    STATc     UNIV. 


C  *«»« 
C**** 
C 
C 

c 

c  **** 

c 

c 
c 


«*******?* *.*********( 


CATA    DESCRIPTION 


I. 


NO.  OF  PROCESS  STREAMS. 
OVEkALL  UNIT  IF>3 


NO  . 


CF    UNITS.     AND    NO.    Of     STREAMS     li 


STREAM 
S( If Jl 

S(  1,  I)  - 
S(L»2)  -- 

SI  1,3)  -- 
S( I r4J  -- 
S( I ,5)  -- 
SI  1,6)  -- 
SI  I  ,71  — 
S(  I ,8)  — 
St  1,9)  -- 
St  I ,101-- 
St  I,  III— 
St  1  , 121- 
St  1  ,13)  — 
St  1,14)-- 
St  I  ,15)-- 
S(  I,lo)-- 

sti  ,m-- 

S( I  ,  13)-- 
St  1,19)-- 
Sl  1,20)  — 
St  I  ,21)-- 
S(  1,22)  — 


COMPCNENT  ARKAYll  HR  OASIS,  ALL  AMOUNTS  IN  TONS)-- 
nHERE  I-NO.  OF  STREAM,  AND  J=CCMPCNENT  CF  PROPERTY 
--    COAL    OR    LHAk     IF     VALUE    NEGATIVE 


OIL 

NITROGEN,     N2 

OXYGEN,    02 

riATER,    H2C 

CARdON    DIOXIDE ,    C02 

C ARSON    MONOXIDE,     CO 

HYDROGEN, H2 

METHANE,     CH4 

ETHYLcNE,     C2H'. 

ETHANE,    C2Ht> 

STATE     (l.=GAS, 

TEMPERATURE,     DEGREES    F 

PRESSURE,     PSIA 

PROPYLENE,     C3H6 

PROPANE,     C3H3 

NORMAL    BUTANE,     C4H10 

HYDROGEN    SULFIDE,     H2 S 

AMMONIA,    NH3 

SULFUR,     S 

ASH 


LlW.    UK    SOLID) 


C    3. 


PROCESS    UNIT    ARRAfS     (1    HR    BASIS) 

UNAMEI!,J)     WHERfc    I  =  NO.    OF    UNIT,     J-NAME    OF     UNIT 

UNlTtl.U)     WHrRE     I=NL.    OF    UNIT 

UNIT! I, 1-8) —    NC.     CF    STREAM    ENTERING    (  *)     OR    LEAVING     1-)     UNIT     I 

STKEAM    NO. <1 .--UNjT IL UcO    STRtAM,     STREAM    NO. MOO. 

—  VALUE    ALSO    INCLUDED     IN    UTILITIES 
UN  IT  (  I  fS>  )     --    FUEL    GAS,    MM    BTG 
UNITtl.lD)--    SnAFT     *CRK,     1000    HP     II.N-*) 
UNIT  (  I  ,11)--    NET     STEAM,     MM    3TU     (Iii«*l 
UNIT (1,12)--    NET    CONDENSATE,     TONS     (OOT-H 


C 

C 

c 


uNITU,l3) —    NET    COOLING    to  AT  ER  ,     lojj    GPM     I  IN-*  J 


#**«:?tt$»9ft»*«#»9««**«.i 


c 

c 
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c 

c 

«       OfcCLARt    ANO    INITIALIZE    ARRAYS                                                                                * 

REAL    UNAME120.30I 

REAL    UNIT<20»50)  rPHYSUOl  »CHE.M(  tJO  J  ,EFF(2  0),£(dO),Sl3D,25) 

REAL     BttYSI 30)  , CHAM (30) , AF F 1 20 ) , A(  80 )  , Or F ( 20 ) iSNI T ( 2  0 , 1 i I 

REAL    X<80i25J.ANIT(20il31  , ETOT I  20 ) , 3T3T ( 20) 

REAL*d    RWrlEAOl  Z3  J  /  "CUAL-TUN  '  ,  '          OIL       ','          AIR       ', 

*'          N2          '  ,  '          02          •  ,  •          H20        '  ,  '          CJ2       •  ,  •          CO          •, 

*'          H2          ','          CH4       •,'        C2H<,        ','       C2H6       ','        C3no        ', 

*'        Cthti        '  , 

*'        CHIO    ','          H2S        '.'           NH3        ','             S           '  ,  -ASH-  TUNS'  . 

*'CHAk-TC:M'  ,  'TOTAL-KH'  ,  'PRES-Pil  ■  ,  '     TEMP-K    '/ 

OAT A    S, PHYS.CMEM.EFF , E, UN  IT ,UNaME/ J  860* 0. / 

DATA    OHYS iCHAP, AFF , A,3FF/23  0*0./ 

DATA    X/2000«l./ 

c 
c 

c 

DATA    ANIT.ONIT.ETUT  ,  BTOT/ 5o  0*0.  / 

*       READ     IN    CATA                                                                                                                               * 

READ(5,10J)     NSTRM.NUNIT.OVERAL 

Re  AO  (5, 200)     1  I  j 1  I , J  )  ,  J=  1 ,  22 ) ,  I=1,NSThM1 

REA0I5.300)     I  1'JNAMEIK,  I  )  ,1  =  1,71  ,  {  UNIT  U«M  J  ,H»l  ,  13  J  ,  K»  I,  NUN  1 1 J 

If  (OVE*AL.EU.O.)     00    TO    2 

K=NUN1T 

c 
c 
c 

READ!  5,40  0)     t  UMTl  K,M)  ,M«l*  ,*5] 

*       ECHO    CHECK    OATA                                                                                                                       * 

2    WRITEI6, I ) 

».ftlTE(o,50) NSTRM.NUNIT 

WRITEI6,  I6<t) 

WrITE(6, 123) 

UO     153    I=1,NSTRM 

rtKlTE(6, 150JI I ,IS( I,  J »,  4- 1, 22 11 

IT  l  I.NE.-»2  )     GO    TO     153 

*RITE(6,lob) 

WRITEI6, 123) 

153    CONTINUE 

WRITEI6,  167) 

WRlTt(6,165) 

*RITE(6,  180) 

WRITE  (6,  250)     (  (UNAMEU,  I)  ,  I  «l  ,7  )  ,  («Nl  T(  K,«l  ,M*l»l31  ,K-l  ,  hU.il  T) 

If luVEKAL.EU.O.)     CO    TC    70 

K=NUNIT 

HRITE(6,25U(UM  T(K.,M)  ,1=  !<,,<,!>) 

70    CONTINUE 

c 
c 
c 
c 
c 
c 
c 

r«RITEt6,185) 

*♦**»*******»****»********«=  ******#«*****»*«-*»#**-e#*»**#****** 

■"ATERIAL     BALANCE    CnECK 

*****»*********•**»**»»♦*«»*«***«♦**********«»***«***«  .  w  .  ,  -   . 

WRITE (6 , 174) 

«RlTE<fc  ,175) 
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*  MOVE    AIR    VALUE     SI    REPLACE     IT    rflTH    CHAR    VALUE     IF    NECESSARY     * 

00  73  l=l,NSTKrt 
S( I  ,24)=S(  I  ,Zl\ 
S( I ,221=0. 

1FISI  Ii  L)  .GE.O. )    GO    TO     73 
S(I , 22) =-S (  1,1) 
SUti  1  =  0. 
73    CONTINUE 

00    oO    K=l .NUNIT 

*  iERU  MATERIAL  BALANCE  IN  C  OUT  VARIABLES  * 

T  IN  =  0. 
TOUT  =  0. 

*  FINU  THE  TOTAL  MASS  ENTERING  £  LEAVING  THE  L.NI  T  * 

DU  65  N=l  ,  4  5 
M  =  N 

*  SKIP  ITEMS  9-13  * 

IF(M.Ct.9l  M  =  M*5 
U  M  T  =  U  N 1 T U,M) 
1FIUNT.EJ.0.)  GO  TO  65 

«   PUT  STREAM  NO.  IN  PRCPER  FGcM  « 

IFl ABSIUNIT  U,M)  ).LT. 1. )     LNT-UNT«103. 
IK  ASS  I  UNIT  (K..M)  )  .GE.10  0.  )     UNT=UNT/  U. 
IU=UNT 

*  CORRECT     FOR    COMPUTER    RUJNOCrh    ERROR    OF    utCIMALS  * 

AA=ASS(UNT  )  +  .4 
J  =  AA 

IK J.LE .NSTRMJ     GO    TU    67 
WRITE16.299)     J 
GC    10    65 
67    CONTINUE 

*  ADO    VALUES     IN    STREAM    TC     [N    08    OUT     VARIABLE  * 

DO    6d    LL=1 ,21 
L  =  LL 


SKIP    STREAM     ITEMS     12 


14 


IKL.GE.12)    L  =  L*3 
IF( IU.LT.O)     GO    TC    69 
TlN=TIN  +  AdSI S( J,L) I 
GO    TO    68 

69     TOUT=TGUT*ABS! S( J.L) ) 

6d    CONTINUE 

65    CONTINUE 
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c     *  sobtract  total  ogt  valle  from  total  in  valge  * 

c 

3ALCHK=TIN-TGUT 
rfRlTE<6,17&HUNAME(K,l)  tt*l  i7)  ,BALCHK 

60    CONTINUE 

c 

c 

C  UTILITY    CALCULATIONS    £.    UNIT    CONVERSIONS 

C 

c 

C  THE    AiNIT    ARRAY     IS    USED     IN     THE    ENERGY    dALANCE    CHECK.     * 

C  »       THE    BNIT    ARRAY     IS    USED     IN    FIkST    LAV,     IFFICIENCY    CALCULA-       * 

C  »       T10N.        THE    UNIT    ARRAY    BECOMES    THE    WORK    EQUIVALENT    FOR  » 

C  »       USE    IN    THE    SECUNO    LAW     EFFICIENCY    CALCULATIONS.        THc  * 

C  *       CORRESPONDING    ELEMENTS     IN     THt    FIRST    TWC    ARRAYS    ARE  * 

C  *       QEFINFO    THE    SAME    AS     IN    THE    UNIT    ARRAY.        THE    VALUES    ARE  * 

C  *       CONVERTED    TO    kCAL.  * 

C 

7    DO    3    J,M  =  1,NUNIT 

ANITt JM,9)=    UNITUM,9)*.2  52«10.»*6 

BNITI JM,9)=ANITI JM,9) 
C 

C  *       ASSUME    FUEL    GAS     PRIMARILY    CH4     THUS    STANOARD    STATc    WORK  * 

C  *       EQUIVALENT    ■     .93*(HEAr    OF    CuMBUSTIGN)  * 

C 

UNIT! JM,9)=UNIT< JM,9)». 93*. 252*1 0.«»c 

ANIT(JM.IO)    =    UNITl  JM,  101  *6 '»!  .4»LQ0Q. 
C 

C  *       USE    AN    ELECTRICITY    TG    WORK    CONVERSION    EFFICIENCY    OF    0  .  d       * 

C 

1F|UNIT(JM,10).GT.0.)     .)NlT(JM,lJ)-O.MIT(JM,LO)*t^l.i*lCLj./.o 

IF1UNITI JM,10) .GT.G. )     JNI  T(  JM  ,1  1)  =JNIT(  J«  ,  L0)*6<U  .•»  *1  0  )  3  ./ .  d 

IFlUNITUM.lOJ.LT.G.)    BNI TUM ,  10)=UNI T( JM, 101*641. 4*1  CJO.*. J 

I  Ft  UNIT  (  JM,  10)  .LT.O.  )     JN  I  T  (  JM  ,  10  )  =  UNI  Tl  JM  ,  1  0  )  *6  4  1  .-.*  100  0  .  *  .  d 

ANITI JM, U)=UNITI JM, 11] *.252*10.**6 

BNITI  JM,U>=A.NIT(J,M,  11) 
C 

C  *       USE    411.5K    L     50    PSiA    AS    THE    CONDITIONS    FUR     STEAM    CALCULA-* 

C  *       TIuNS 

C 


o.tl  1.5)  i  ♦ 


C 

C  *      CONVERT    FROM    DECREES    F     TO    K 

C 

S(IT,13)  =  (StIT,U)-32.)/l.d    '27  3. 15 
27    CONTINUE 
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*  CONVERT     THE    STREAM    COMPONENTS    FROM    TONS    TO    K.MQLS 
DO    23    L=1,NSTRM 

*  USE    A    GIVEN    OIL    MCLbCULAR    nclGHT    OF    300    K.G/KMUL 


S(L 
SIL 
S(L 
S(L 
SIL 
SIL 
SIL 
SIL 
SIL 
SIL 
SIL 
S(L 
SIL 
SIL 
SIL 
SIL 
SIL 
23    CCN 


II 


30 


32 


33't 


92 


33 

61 


,2I=S 
,3)  =  S 
,4)  =  S 
,5)  =  S 
,6)  =  S 
,7)  =  S 
,8)  =  S 
,9)  =  S 
,10)  = 
,11)  = 
,151- 
il6>* 
,  L7)  = 

,ld)  = 

,19)  = 
,201  = 
,241- 

TINUE 


IL.2) 
I  L  .31 
(L»41 

(L,5) 
(L,6) 
(L,7) 
I  L  ,8) 
(L,9) 
SIL.l 
SIL  ,1 
SlL,L 
SlL,l 
SIL  ,  L 
SlLrl 
SIL  ,1 
SIL, 2 
SIL, 2 


♦3.024 

♦32.385 

♦28.351 

♦50.357 

♦20.61* 

♦32.383 

♦450.013 

♦56.548 

0)^32.337 

11*30. 169 

5)*21.553 

6>*20.573 

7)*15.608 

8)^26.62 

9)^53.269 

0)^28.293 

4)^23.99 


TUTAL    KMGLS    AND    CALCULATE    HOLE    FRACTIONS 


THE  TUTAL  KMOLS  OF 
SISTR.  NO. ,23)  AND 
X    ARRAY. 


EACH    STREAM    AkE    PUT    INTO 

THE    MU.E    FRACTIONS    ARE    PUT     INTO    THE 


DO    32    LT=1 .NSTRM 

00    31    LA=2, 11 

S(lT,23)=S(LT,23)*A3S(S(LT,LA)) 

CONTINUE 

00    30    LA=15,20 

S(LT,23)=S(LT,23)*A3S(S1LT,LA)I 

CONTINUE 

5(LT,23)=S(LT,23)+S(LT,24) 

CONTINUE 

00    81    J=l, NSTRM 

IFISI J,23) .GT.O. )     GO    TO    334 

Xt  J,K)=l. 

GO    TO    8  1 

JO    82    x.  =  2  ,  LI 

XI  J,K)=     SI  J,t<l/S(  J, 23) 

IFIXI J,iO  .Lt.O. )     X  U,K)  =  1. 

CONTINUE 

OU    83    K-  15,20 

XU,K)»    SU,K.)/S  (J.23) 

IFIXI J.KI .LE.O. )     XI J,<) =1. 

CONTINUE 

CONTINUE 


****#**«•«***#* 


***»******•»****#***««  *****»*****s»*-?«i**#***©« 


CALCULATE    DATA    NEEOEC    TL    PeRFiiiiM     THE    ENERGY    BALANCE 
AND    TO    DETERMINE    THE     1ST    £,    2ND    LAW    EFFICIENCES 


c 
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c 

fc***************^*************************^**************** 

c 

c 

*                   AFTER    SOME    PRELIMINARY    CONSTANTS    Art E    ESTABLISHcu, 

* 

c 

*       VALUES    ARt    CALCULATtO    FOR    VARlASLtS    CORRESPONDING    TO 

* 

c 

*       EACH    COMPONENT     IN     THE     STREAM — FUUR    VARIABLES    PEk    CCM- 

* 

c 

*       PUNtNl:        c|     COMPONENT     Atl  BR  EV  I A  T  [UN     I        USED    FUR    THfc    CHEM- 

« 

c 

*       ICAL    wORK.    EQUIVALENT,     Pi     COMPONENT    ABBREVIATION     1       USED 

* 

c 

*      FOR    THE    PHYSICAL    hJRK    EQUIVALENT,     A|     COMPONENT    ABdREVIA- 

e 

c 

*       TIGN     I       USED    FOR    TllE    HEAT    OF    COMBUSTION    AT    REFtKtNCE    CON- 

-* 

c 

♦       OITION,     AND    B|     COMP.     A2b.     1       USED    FOR    THE    ENTHALPY    CHANGE* 

c 

*       FROM    REFERENCE    CONDITION. 

* 

c 

*       THE    CPI     AND    CPA    FUNCTION    SUBPROGRAMS     ARt    USED    TO    FIND 

* 

c 

*       VALUES    FOR     ThE    P    AND    B    VARIABLES    RESPECTIVELY.       CPI     ANO 

* 

c 

*       CPA    CALCULATE    THE    CHANGE     IN    "CRK    EQUIVALENT    AND    ENTHALPY 

* 

c 

*       RESPECTIVELY    AS    A    COMPCNtNl     IS     ADJUSTED    TO    THE    STANDARD 

t 

c 
c 

*       STATE     1298.151. 

» 

DO    125    [*ltNSTRM 

T=St  I  ,li) 

TI  =  293.  15 

RTI  =     1.9o7*TI 

c 

c 

HVCOAL    ■     I2O0O. 

»      USE    THE    FEAT    OF    CQMSUSTI GN-«0»rt    EJOIVALINT    APPROXIMATION 

« 

c 
c 

c 
c 

*       hITH    THE     12000     BTU/LQ    VALJt    GIVuN    FOR    COAL 

» 

EC0AL=HVCGAL*SII,1)*5Q4. 

*       TAKE     THE    AVERAGE    HEAT    CAPACITY    FUR    COAL,     CHAR,     AND    ASH 

* 

c 

•       TO    BE    G.224    bTU/LB/OEG    F    AS    GIVEN    BY    PERRY     IN    ThE    CtiEM- 

* 

c 

c 

*       ICAL    tNGINtEK'S    HANGBuC*. 

* 

PCUAL  =  S  I  I  ,  I  I'lU.M    *  !  T-T  I     -II        «ALQG(T/TI     ))*l.d 

ACOAL=    St  I  ,  1)*HVCGAL*504. 

c 
c 

BCCAL=    Silt  U*(T-Tl)  *U2.S*i.<i 

*       CHECK    FOR    TYPE     OF    OIL     AND    CALCULATE    E     £.    A    VALUES 

» 

c 
c 

♦       ACCORDINGLY. 

* 

HVOIL    =     15300 

IF(S( I ,2) .G6.0. )     GO    TO    14 

S(l  ,2)=    -SI  I  ,2) 

EUIL=     L9100.»S( I ,21*166.67 

AUIL=    19100. *S( I ,2) »l6o.o7 

c 

c 

GO    TO    19 

*       USE    THE    GIVEN    SYNCRUOE    HluHCR    HEATING    VALUE    OF     19100 

• 

c 

*       BTO/Lu    ALONG    k»ITH    AN    APPROXIMATE    VALOE    uF    14700    BTU/Ld 

« 

c 
c 

*       FOR    THE    RAb    OIL. 

* 

14    EOIL=S(  I  ,2)*HVCIL«16o.  j7 

c 
c 

AU1L=    St  1  ,2  )*Hv>CIL*U6.  7 

*       CHECR    OIL    FUR    LljUID    OR    GAS    STATE    AND    CALCULATE     P    £    6 

* 

c 

*       VALUES     ACCORDINGLY, 

* 

c 

*       ASSUME    OIL    PHYSICAL    PROPERTIES    TO    BE    CLOSE    TC    THOSE    OF 

* 

c 

*       ANTHRACENE:       CP=- 14. 09 ♦. 20 4* T     KtAL/KMDLi     HtAT    OF     VAPC.-U- 

* 
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c 

*       ZATICN*135    KCAL/MOLEi    A, ID    VAPOR    PRES  SURE  =  .  CC2    HM    Hj    AT          * 

c 

c 

»       293.15     K.                                                                                                                                        * 

19 

P0IL«135.*SII.2)*(T-    TI    -    Tl    *AL3G(T/    Tl    )) 

8UIL=    Sll  ,2)*135*(T-TI I 

1F(S(  I,  12)  .NE.l.  )     GO    TO    11 

P01L=S(  I  ,2)*(  1.937*ri*ALOG(  38  CO  00.)  *CPI  (-14.09  ,204.2,  T)  t-RTI* 

"ALOGIXI I ,2) ) 1 

130  [L  =    S(I,2)*ICPA(-14.09,204.  2,  T)  ♦137.*I66.7] 

11 

EN2=S(I  ,3i*lh3. 
PN2*SU,3)*ICPH6.83,.9,T)*JtTl*ALUG(XU,3n  J 

AN2=    0. 

8N2=    5(1,3) *CPAI6.83, .9,T ) 

E02=S(l  t4)»932. 

P02=S(I  ,4)*(CPlir. 16, l.,T>»-Rr  I* ALUG  1X11,4)1) 

A02  =    0. 

c 
c 

302  =    S( 1,4) *CPA17.16, 1. ,T ) 

*       CHECK    hATER    FCR    LIUU.10    OR    GAS    STATE     AKC    CALCULATE    P    E    3       « 

c 
c 

*       VALUES    ACCORDINGLY.                                                                                                            * 

IFISI  I, 12). 60.1.)     GO    TO    10 

EH2G=0. 

AH20=    0. 

ph2c=s(  i  ,5)«(  r-Ti-n*ALOcir/r  i)  i*id. 

BH20  =    S ( I , 5 ) *1 d . «( T-T I ) 

GO    TO    20 

n 

EH2G=0. 
AH2G=0. 
PH2U=S(  I, 51*12054. 5<-CPI(7. 30, 2.  4o,T)     *RTI*ALQG  (X(  It  Si  )  1 

OhZU=    S(  I  ,5)*l  10519*-    A3SICPAI  7.  j,2.46,T)  )  )M  T-TI  1 /AdSiT-TI  1 

iO 

EC02=Sl I ,61*4776. 

PC02=S( I ,6) *(CP1 ! 10.57,2. 

10 ,T)*RT I*ALOolXI I  ,o)  )  ) 

AC02  =    0. 

3C02  =    S( 1 ,6l*CPA(lC.57,2. 

10  ,T  ) 

ECO=S 11 ,71*65770. 

PCO  =  S (1  ,  /) »(CPI (6.79  ,  .9d, 

T)  t-RTi  *ALOGI  XI  I,  7)  )  ) 

ACU=     S( I ,7)*67o36. 

6CJ=    SI  I ,7)*CPAlo. 79, .98,  T) 

EH2  =  S ( I  ,8)*5o224. 

PH2  =  S(  [  , 31 *(C?I I o.  52i.T8»T)+RTI*ALOGtX(  1,8)  )) 

AM2=    SI  I ,8 ) *6d31 I. 

3H2=    S( I ,3)*CPAIo.52 ,.73,T) 

ECH4=S( 1,9) *193420. 

PCH4=S(  I  ,91*1 CP 113. 331,  13.0  4t,T  )  CRT  I*AL0GIX  1  I  ,  9)  )  ) 

ACH4=    S( I ,9)*2l2dC0. 

SCH4=    S  (  1,9  )*CPA!3.  331,  18.044,  T  ) 

EC2H4=S( I ,101*303672. 

PC2H4^SI  1,  10)*(CPl!2.d3,  23.6,  Tl  »-R  T  I  *ALUG  (  X  I  1,10)1) 

AC2H4=    Sll, 101*337150. 

3C2H4=    S(  [  ,  10)*CPA(2.83  ,28.6,  T) 

EC2H6=S(  I, 111 *35703o. 

PC2H6  =  S(  I  ,  11  )  *(  CPU  2.  247,  33  .2,T  )  t-rtT  I*AluGU(  I  ,  11)  1  ) 

AC2H6=    SI  I  ,  11  1  *372820. 

uC2ri6  =    Sll,  It)  *C  PA  12.  247,  33.  2,  T) 

tC3H6=S  1  I  ,  15)*44501a  . 

PC  3H6=i([,15)*(CPI13.  253,  <»5  .  12,  T  l+R  T  I  »A  LOG  I  X  11,15)11 

AC3H6=    SI  1 ,  15 1*491990. 
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dC3H6 
EC3H8 
PC3H8 
AC3H6 
SC3M8 
EC4hl 
PC4HI 
AC4H1 
BC4H1 
£H2S  = 
PH2S  = 
AH2S  = 
QH2S- 
ENH3  = 
P,\H3  = 
ANH3  = 
3NM3  = 
ESLFR 
PSLFR 
ASLF  R 
3SLFR 
EASH  = 
PASH= 
AASH= 
BASH" 


SI 
=S(  I 

=  S(I 

=  si 

=  S( 
0=S  I 
0=SI 
0=  s 
0=  s 
SI  I  , 
SI  I, 

S(  I 

51  I 
SI  I, 
SI  t  , 

0. 

SI  I 
=S(  I 
=  S(  I 

Sll 
S( 
0. 
S(  I  , 

0. 

S(  I 


I  ,  15>*CPA(  3.253,  45.12!  T) 

rl(j)*5l36l8. 

,  16)*(CPI(2.4l,5,7.2,  Tl*RTl*Al.aG(X(  I,  L6>)  I 

[,  16 1*530600. 

I  ,16)*CPA(2.4l  ,57. 2, T) 

I  ,  17)  *6698L4. 

I,17)*(CFl(3.d44,73.35,T)  t-RT  I  *AL0G  I  X  I  1,17))) 

{  I  ,171*687640. 

(  I  ,17)*CPA(3.d44,73.35,T) 

181*187992. 

ld)*l CP [17. 81, 2. 96, T)*RTI*ALUG 1X11,18)1) 

,18)  *L3<tt62. 

,IB)*CPAI  7.81,2.96  ,[  ) 

L9)*30432 . 

19)*(CPI(7.U,6.,T)i-Rri*ALGG(X(I,19))l 

,  19)«CPA(  7.U.6.  ,T  1 

,201*139660. 

, 20>*CP II 3. 5  J, 6. 24, T ) 

,20) «70960. 

I  ,20)*CPA(3.53,6.24,  T) 

21)*L12.9*IT    -H-l  1*ALGGIT/T  I)  1*1.3 

,21  1*112.  9*1.  3*(T-n  ) 


USE    T*U    DIFFERENT    CHAR    HEAT    GF    CQM30STIUN    VACJE S I dT J/ Ld) 
DEPENUING    JPCN    STREAM    NUMBER. 


1FI1. 
If  I  I  . 
ECriAR 
PCHAR 
ACHAR 
8CHAR 
EAIR  = 
PAIR  = 
AAIR  = 
oAIR  = 
PHYSI 

*  PC2H 
CHEM( 

•EC2M6 
BHYS1 

*8C2H6 
CHAMt 

♦AC2H6 


EG. 12)  HVCHAR 
Ew.L3)  HVCHAR 
=S( I ,221*504. 
=  S(  1  ,22)  *U2. 
=  S( I ,221*504 
=  S(  I  ,221*112 
0. 
S( 1 ,24)*CPI (6 

0. 

SI  I ,241 *CPAl 
11  =  PCQAL+PU1 
6<-P^3H6»PC3H8 
11=  ECGAC+EOI 
♦  tC3H6»-EC3Ha  + 
II-  3CGAL+BGI 
*8C3H6*8C3H8  + 
1)=  AC0AL*AG1 
*AC3H6<-AC3H3«- 


=12000. 

=11100. 

*HVCHAR 

9*( T-Tl-T  I*ALGulT/TI)  ) *l. 

.*HVCHAR 

.9*Ld*ir-T  I  1 

.9  ,.92,  T) 


6.9, .92 ,T ) 

l_*PN2  +  PC2+Ph2U«-PL02+PCOPH2  +  PC.44*PC 

+PC4H10*PH2S*PNH3*PSLFK*PASh*PCHAK* 

L«-EN2*EC2*-EH2C  +  EC02*ECC+EF2»ECH4*EC 

EC4HL0»-£H23<-ENH3«-ESLFR  +  EASH  +  ECHAK 

l  +  BN2*BC2*3H2U*dC02*3Ca*BH2«-flCH4fBC 

BC4H10*BH2S+BNH3+BSLrR*BASH+dCriAK+3 

i_*AN2-i-AC2-t-AH2C»-ACG2*ACU-i-AH2«-ACH4t-AC 

AC4H10*AH2S  «-ANH3  +  ASLFR*AASri<-ACHAR 


2H4  + 
PAIR 
2H4» 

2H4» 
AI  K 
2H4  + 


IFtSI  1,12) .Nfc.l.JGC    TG    74 

PHYSI  l)=PHYS([)*SII,23)*l.987*TI*ALUG(S(I,14)/14.7) 

*  CALCJLATt    THE    10TAL    WGRK    EQUIVALENTS.     Eli),     AND    ThERMAL        • 

*  ENEnGY,     At  I  )  ,     FUR     STREAM    I.  * 


74    El  I )=PHV3I  I )+CritM(  I ) 
Al  I  )  =  bhYSl  1  1  t-CriAMl  1  ) 
125    CONTINUE 
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|c  ' 

c 

****♦**#********$*****<=*******?**#<?#*******♦***************** 

c 

c 

CALCULATE    HIRST    #N0    SECOND    LA«     EFF  IC  I ENC  IES 

c 

c 
c 

AMU    ENERGY    BALANCE    CLOSURES 

It *#^****  ♦*****♦****»♦»♦**♦***  *»*****»*«*#♦*******♦****#***»*• 

c 

c 
c 
c 

OU     130    K=l ,NUN[T 

*       ZERC    THE    FIRST    AND    SECONC    LAW    EFFICIENCY    VARIABLES                     * 

EIN=0. 

EOUT=0. 

AIN=0. 

BIN=G. 

AOUT=0. 

60UT=0. 

00    151    M-1.45 

M=N 

c 

c 

c 

*       SKIP     ITEMS    9    -     13                                                                                                                 * 

IFIM.GE.9)     M=M+5 

UNT=ONI TU,M) 

ABSU=ABS(UNT) 

c 

c 
c 

*       RUT    STREAM    NU.     IN    PROPER    FORM                                                                               * 

IF(  ABSIUNITIK.M)  )  .LT.l.i     UN  T=  UN  T*  LO  D  . 

IF(A3S(UNIT(K,M) ) .GE.IOO. )     ONT=0. 

IU=UNT 

c 

c 
c 

*       CORRECT    FOR    COMPUTER    ROUND    OFF    OF     DECIMALS                                           ♦ 

AA=AdS(UNT  )  *.* 

J=AA 

Ir(  J.NE.T)     GO     TO    1 9 

EUN1T=0. 

AUNl T=0. 

3UNIT=0. 

GO    TO    17 

29    EUNIT  =  EU) 

BUNIT=A( J ) 

AUN1T=A( J) 

c 
c 

1F(  IO.LT.O. ]     GC    TO    33 

*       CALCULATE    THE    TOTAL    USEFUL    nORK    EQUIVALENT,     EIN,     THE                « 

c 

»       TOTAL    USEFUL    THERMAL    ENERGY,     BIN,     AND     THE    TOTAL    THERMAL       •> 

c 
c 

«       ENERGY,     A1N,     ENTERING    THE    UNIT    AS     STREAMS.                                            * 

L7     1F< ABSG.LT. L )     EUNIT=3. 

IFIABSU.LT.  L  )     3UNIT=0. 

fclN=ElN*EUNIT 

B  IN=B  INfBUN  IT 

ain=ain*-aun;t 

GO    TO    151 

t 
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C  *       CALCULATE    ANALUGOUS    QUANTITIES    LEAVING    THE    UNIT    A5  * 

C  *       STREAMS.  * 

C 

ii     IF1A113U.LT. 1.  I     EUNIT  =  0. 

IF1A8SU.lT. 1. )     BUNIT=0. 

EUUT=EGUT+EUNIT 

BGUT=BOUT+SUNIT 

ACUT  =  AOUT«-AUMT 
151    CONTINUE 

H=0. 

Y  =  0. 

1  =  0. 

DO    47    KM=9,13 
C 

C  *       CALCULATE    THE    TOTAL    WChK    EQUIVALENT,     «l ,     THERMAL    ENESGYi        * 

C  *       Y,    AND    THtRMAL    ENERGY     INCLUDING    THE    GtNtRATING    EFFICIENCY* 

C  *       I,    FOR    EACH    UNITS    UTILITIES.  * 

C 

IFIUNITU.KM)  .LT.O.  )     GO    TO    4o 

W=W+UNIT(K,KM ) 

IF1ANIT IK, KM) .LT.O. )     GO    TO    46 

Y=    Yt-ANI  D  K  ,KM) 

Z=Z«-BNIT(K,KM) 

GO    TO    47 
C 

C      *   CALCULATE  TOTAL  QUANTITIES  LEAVING  A  UNIT.  * 

C 

4o  EUUT  =  CCUT«-A8SiUNIT(K,KM)  J 

AOUT  =  AUUT*ABS(ANIT  ( K , KM  )  ) 

B0UT  =  B0UT«-A6S(dNIT(K,K1)) 
47  CUNTINUE 

3»EIN+«' 

IF10.NE.0. I  GO  TC  43 

EFFIK  1  =  0. 

afF(K)=0. 

GO  TO  4* 
C 

C  *       DETERMINE    THE    SECCNO    LArt    EFFICIENCY,     EFFIK),    THE    ENERGY       * 

C  *       BALANCE    CLOSURE,     AFFIK),     AND     THE    FIxST    LAW    EFFICIENCY  * 

C  *       FOR    UNIT    K.  * 


C 


EFFlK)=EOUT/(EINi-W  ) 
AFr(K)  =  AOUT/(AIN*-Y)«lUO, 
3(-F(K)  =  dUUT/l  BIN+Z) 
BT0T1K)  =  ElNt-Z 
EIN+W 


ETCT1  K) 

4'/ 

EOUT=0. 

ACUT=0. 

BOUT=0. 

E1N=0. 

AIN=0. 

3IN»0. 

n  =  0. 
Y  =  0. 
Z=0. 

13U 

CONT  INUE 

C 

c 

*       ADJUS 

ADJUST  ARRAY  ELEMENTS  FUR  CUTPJT. 


c 

106 
00    51    K»XfNST8M 

rs»$(K,i3) 

PS=S IK,  14) 

DO    52    M»12  ,20 

S(K,N»«S(K.N+.3) 

52 

CONTINUE 

SIK,21)=PS 

S(K,22)=TS 

00    87    M=l ,20 

N=24-M 

SIK,N>SU,N-1) 

87 

CONTINUE 
S(K,3)=S<K,24) 

51 

CONTINUE 

NS  =  0 

M=l 

MF    =     10 

&d 

WRITE16.310) 

WRITEI6.320)     (J,J=M,MF) 

wRITE(6,330) 

WRITE(6,350) (RWHEAQ(J)  ,  IS (I  ,  J ) , I "M , MF ) , J= 1, 2  3 ) 

WRITE (6,551 KCHAH(N) ,N=M,MF ) 

WRITE (6, 651) (BHYSIM , N=M,MF) 

WRITElc,  751  MAI J  J)  ,JJ=M,MF) 

WRITE (6,550) (CHEM(N) ,N=M,MF  ) 

«RITE(6,650) (PHYSIN) ,N=M,MF) 

WRITE (6, 750) I  El  J J)  , JJ  =  M,MF) 

NS    =    NSHO 

IF(NS.GE.NSTRM)     GC    TC    89 

M    =    M  +  LO 

MF    =    HF*iQ 

GO    TO    88 

89 

WRITE (6 ,423 ) 

WRITE  (6,450)  (  (UNAMblK,  I  )  ,  1  =  1  ,  7 )  ,  3T0T  1  K.)  ,  ETOT  (  K )  ,  K.=  1  .NUNITI 

WR1TE16.670) 

WRlTE(o,675) 

WRITE(6,680) 

c 
c 
c 
c 
c 
c 
c 

WRITE  (6,477)  I  [UNAMEU,  U  ,  1=  1,  7)  ,BFF  IK.)  ,  Ef  F(  K)  ,  AFF  IK)  ,  K"  I,  WUtfTT) 

*»*******»****»»»» ******»»**»*S »*****»*»»**»¥»»»*«***»* ***>«* 

FORMATS 

v**************************************************?*^******* 

L 

FORMAT  CI','                EChO    CHECK.    OF    DATA     '//) 

50 

F0RMAT15X, ' NO.     STREAMS=     ',12,'       NO.    UNITS-     ',12//) 

100 

F0RMAT(2I 10.F10.4) 

123 

FORMAT  (IX,  'NO.  '  ,3X,  'COAL'  ,2X,  'U  If  ,  5X,  'N2',  4X,  '02'  ,4X,'H  20'  , 
*4X,'CQ2'  ,5X,'C0'  ,4X,'H2,,3X,'CH4'  ,1X,  'C2H4'  ,2X,'C2H6',1X, 
*'G»',2X,,TEMP',lX,'PREjS',lX,'C3H6',lX,'C3rid',lX,'C4H10l,LX, 
*'H2S'  ,2x,'NH3'  ,3X, 'S' ,4X,  'ASH',4X,  'AIR' /  ) 

150 

FORMAT! lX,I2,lX,F7.1,2F7.2,iFb.2,2F7.2,2F6.2,2F5.2,F3.),Fo.O, 
«F7.1 ,6F5.2,F6.2,F7.2) 

164 

FUKMATI52X, 'STREAM    -    COMPONENT    CATA**1/) 

165 

FORMAT!  •  1  '  ) 

lo7 

FORMATl/     IX,  ■*    GAS-LIQUID    I  CENT  I  F  IC  AT  luN       ■'!•■=    GAS,     ''0"', 

/     IX, '♦♦    All 

107 
RATURE     '  , 

*  ■  ■    L  I  UJ  I  U  • 

COMPONENTS     IN     TGNS/HR,     TEMPE 

••IN    OEG    F, 

ANtl    PRESSURE    IN    PSIA'J 

174 

F0RMATI/////////55X, ' MATERIAL    BALANCE    ChECK'/) 

1/5 

FORMAT) 7oX 
*' (  IN-OUT) * 

,'OIFFERENCE 

IN    T0NS/HR'/53X, 'UNIT    NAME', 

ISX, 

176 

F0RMATI43X 

,7A4, 10X.F3.3/) 

ISO 

F0RMAT148X 

, 'PROCESS    UNITS    STREAMS    ANO   UTILITIES*'/ 

52X, 

♦  ■STREAMS    ENTERING'  ,  18X, 

'FEUL'  ,1  IX, 'NET' /     17X,'UMT 

NAME'  ,29X  , 

••AND    LEAVING' .20X, 'GAS' 

,4X,  'WORK'  ,3X,  'STEAM  '  ,4X,  'CO 

NO. ' , 4X, 

•»Ch"//l 

ws 

FOKMAT!/    aX,'*    ENTERING 

.    STREAM    NO.     "*•',    LEAVING    STREAM    NO', 

♦< .     • ■-• ■ , i 

/  10X, 'STREAM 

N0.<1.— UNUTILUEO    S  TR  E  AM,  '/ 10X, 

♦'STREAM    NO 

.MOO.  —  VALUE 

ALSO     INCLUOEO    IN    UTILITIES' 

//10X, 

♦'UTILITIES 

: '/18X,'FUEL 

GAS,    MM    BTU'  ,9X,  •(  IN  =  H  '/18X 

,  '  wORK  ■     '  , 

♦'HP'  , 17X,  • 

I IN=  +  ) '/ldX, ' 

STEAM,     MM    BTU'  ,  12X,'  (  IN=  +  )  '/ 1BX. 

♦'CONDENSATE  ,     TUNS', 9X,' 

(OLT=*J'/ldX, 'COOLING    WATER, 

100  0     ', 

♦'GPM       (IN= 

*■)'} 

200 

FORMAT! 11F7.2) 

250 

F0RMATI8X, 

7A4, 1X.4F7.2, 

5F5.  0,  4Fd.  1/1 

251 

FORMAT! IX, 

21F6.2.4X/1 

299 

FORMAT!  10X 

,  'STREAM    NOT 

INPUTTED — • ,15) 

300 

F0RMAT(7A4 

,4F4.0,5F3.0, 

F5.1  ,F5.0,  2F5.  1) 

310 

FORMAT!  '  1' 
♦//) 

///////       50A, 

'CCED    PROCESS    STREAMS     (     1    hR 

LAS  IS)  ' 

320 

FOkMAT! IX, 

■  STREAM    *»|8X 

.12,9!  IDA,  121/  ) 

330 

FORMAT! IX, 

'  COMPOS  I  TICN-K.MOL  '/  ) 

350 

FORMAT! ( IX 

, A8,3X,  10(F  10 

.2,2X1  1  1 

40J 

FCP.MATI20F' 

4.0) 

42  3 

FORMAT!  •  1  ', 

////////51X, ' 

TOTAL    ENTERING    J  CANT  I  T  I  ES  •  // 

a7X,  'TUTAL'j 

' 

►'     ENERGY     IN' ,4X, 'TUTAL 

WOnK    El,.     I  N  '  /  4  IX  ,  •  uN  I  T    NAME' 

,  17X, 

' 

►'KCAL     (1ST 

LAW) • ,7X,'KCAL     ( 2N0    LArtl'/l 

4  50 

FORMAT (31X 

, 7A4.9X, El  1.4 

,9X,£l  1.4/J 

4  77 

FORMAT! 15X 

,  7A4,UX,F9.5 

,1  1X.F9.5,  12X,  Fv.3/1 

550 

FORMAT! IX, 

'ELH/KCAL' ,3X 

,10E12.*/1 

551 

FORMAT (/  IX 

,  'ACH/nCAL',  iX,  IJE12.4/I 

650 

FORMAT! IX, 

'EPH/KCAL' ,3X 

•  1UE12  .4/1 

03l 

FORMAT! IX, 

•APri/KCAL' ,3X 

r ICE  12 .4/1 

67J 

FORMAT!  '  1', 

'//////29X, 'THERMODYNAMIC    EFFICIENCICS    ANO    ENERGY'. 

••     BALANCE     ( 

:LCSURE    FOR    THE    CCED    PROCESS'//) 

675 

F0RMAT124X 

,  'UNIT    NAME'  , 

22X,  'FIRST    LA*' ,  10X,  '  SECOND    1 

-Aw1   ,  10X  , 

♦' ENERGY    BALANCE' ) 

680 

FORMAT! 55X 

,  'EFFICIENCY' 

,     9X,'  EFFICIENCY'  ,12A,  'CLOSURE     IZK//I 

750 

FORMATQX,  ' 

'  E/KCAL'  ,5a,  10E12  .4/1 

751 

FORMAT (IX, ' 

' A/KCAL' ,5X,10E12 .4/1 

STOP 

EN  J 

108 
FUNCTION    CPU  A,b,T) 

THIS    FUNCTION    SUBPROGRAM    CALCULATES    THh    CHANGE 

IN    WORK    tUUIVALENT/*MCL    Ao    A    DCSPONENT    IS 

ADJUSTED    TC    THE    STANDARD    STATE     (293.  15KI 

C=B/l00O. 

TI=298. 15 

CPl=(A-TI*C)*iT-TI)+C/2*lT»*2.-TI«*2.)-A*TI *ALGG( I/TII 

RETURN 

END 


FUNCTIUN    CPA(  A.B.T ) 

THIS    FUNCTION    SUBPROGRAM    CALCULATES     ThE 

ENTHALPY    CHANGE    AS    A    COMPONENT    IS    ADJUSTED 

TO    ITS    STANDARD    STATE     (293.15K) 

C=d/10C0. 

TI=29d.l5 

CPA=A*( T-TI ) *C/2.*(T«»2-T I»«2 ) 

RETURN 

END 


109 
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MATERIAL    BALANCE 

ChECK 

DIFFERENCE     IN 

TONS/HR 

UNIT    name 

(1N-OUT) 

COAL    PREPARATION 

-0.001 

DRYING    C    STAGE     I    PYRCLYSIS 

-0.0  00 

STAGE    2.3,    mO    *    PYPOLYS  IS 

-O.G0L 

PRODUCT    RECOVERY 

0.001 

PYRGLYS1S    CAS    COMPRESSION 

-0.0 Oi 

PYRULYSIS    ACID    GAS    REMOVAL 

O.OOl 

PYRQLYSIS    GAS    EXPANSION 

0.0 

OIL    FI  LTRATIO.N 

-o.d  oo 

OIL    hYCRCTREAT ING 

0.001 

PURGE    ACID    GAS    REMOVAL 

0.0 

HYCR0CAR3CN    REMOVAL 

u.O 

AMMONIA    REMOVAL 

-0.000 

HYDROGEN    PLANT 

0.  00b 

CHAR    GASIFICATION 

0.0  00 

GASIFIER    GAS    COMPRESSION 

0.0  00 

GASIFI  ER    ACID    GAS    REMOVAL 

-0.  300 

GASIFIER    GAS    EXPANSION 

0.0 

SULFLR    RECOVERY 

-0.0  00 

OVERALL    COED    PROCESS 

0.043 
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APPENDIX  B 
This  appendix  refers  to  the  optimization  of  the  flue  gas  heat  recovery 
system  in  Part  III.   A  computer  program  listing  with  sample  output  corre- 
sponding to  8  heat  exchangers  and  a  1  in.  tube  outside  diameter  is  presented. 
Throughout  the  program  listing  descriptive  comments  are  included  which  give 
most  of  the  assumptions.   The  program  may  briefly  be  described  in  terms  of 
three  sections.   The  first  section  physically  designs  the  heat  exchanger, 
a  counterflow  one-pass  fixed- tube-sheet  shell  and  tube  carbon  steel  ex- 
changer, once  the  inside  and  outside  fluid  properties  have  been  entered  as 
data  and  once  the  independent  design  variables  have  been  given  values.   In 
this  section  the  shell-side  heat  transfer  and  pressure  drop  quantities  are 
calculated  from  empirical  relations  given  by  Perry  (16)  and  the  corresponding 
tube-side  quantities  are  calculated  from  relations  given  by  Peters  and 
Timmerhaus  (24).   Because  of  the  moderate  stream  temperatures,  less  than 
550K,  radiation  has  been  neglected  in  the  heat  transfer  coefficient  calcu- 
lations.  The  next  section  of  the  program  estimates  all  of  the  equipment 
costs  from  correlations  given  by  Guthrie  (27)  and  Peters  and  Timmerhaus. 
The  final  program  section  calculates  the  remaining  quantities  such  as  power, 
annual  costs,  and  work  equivalents  necessary  to  determine  values  for  the 
economic  and  thermodynamic  objective  functions.   The  program  is  structured 
to  print  out  results  for  any  or  all  of  the  objective  function  optimizations. 
The  results  consist  of  two  major  sections.   The  first  gives  a  search  table 
for  a  specified  tube  outside  diameter  and  number  of  heat  exchangers.   The 
temperature  difference  column  listed  in  the  table  shows  the  temperature  dif- 
ference that  exists  at  the  hot  end  of  the  heat  exchanger,  i.e.  inlet  flue 
gas  temperature  minus  exit  preheated  air  temperature.   The  next  section 
presents  the  heat  exchanger  design  criteria  and  objective  function  values 

in  addition  to  several  intermediate  results  for  one  of  the  objective  function 
optimizations. 
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c 

c 

****************************************************************** 

SE 

CUNU    LAW    AND    CUS1     Of  T 1  MI  l  AT  RJ  M    OF    A    ELUE     GAS    HEAI    RECOVERY 

C 

SY 

STEM   USING    A    COUNT  EkfLUW    ON  E    PASS    M XeO-I UdE- SHE E 7    SHELL    A  NO 

c 
c 
c 

TU6E    HEAT    BXtHANCfcK-ALL    CAfJBGN    STElL 

WRITTEN    UY    TEkkY    L.    UNRUH    UNDER    T  hE    ADV  IS  EM  EN  T    Of    BR.    &.    G.    KYLE  , 

c 

DEPT.    DP    CHEMICAL    ENGINEERING,    KANSAS    STATE    EM  V. 

£,***********#********•*:***********$****************$**************** 

g*** 

****************************  *************  4;  ************y**3****# 

c 

c 

c 

c 

c 

c 

DATA    DESCRIPTION 

»****»»**********«******♦♦*«««**♦***«»(:***♦»**««»««*«»;* 

1. 

INSIDE    FLUID    PROPERTIES.    ttl-KG/HRt  MASS    FLCW    RATE),    CPI- 

c 

X.LAL/KG/DEG    KIHEAT    CAPAC  IT  Y  )  ,  V  IS  I-  S/CM/S  (  V  ISCGSI  TY)  ,  IK  I- 

c 

c 
c 

CAL/S/CM/OEG    K  (THERMAL    CONDUCTIVITY) 

2. 

OUTSIDE    FLUID    PROPERTIES,     V.U-KC7  HR  MASS    FLEW    RATc),     CPU- 

c 

KCAL/KC/OEG    KIHEAT    CAPACITY),     V  I SC-G/CM/ S i VI SCC3 ITY )  , TXU- 

c 
c 
c 

CAL/S/CM/DEG    MTH6RMAL    CONDUCTIVITY) 

3. 

RANGES    OF    SEARCH,    NUMBER    OF    HEAT    EXCHANGERS,     TUBE    CUTS  IDE 

C 

DIAMETER,     EXIT     PREHEATED    AIR    TEMP.,     INSIDE    SHELL    DIAMETER, 

t 

ANO    OPTIMIZATION    UUTPUT    CHOICE,    NET    PRESENT    VALUL,     rtC^N 

c 

EQUIVALINTICQMdINEG    HEAT    RtCUVERY     SYSTEM),    OVERALL    SECOND    LAW 

c 
c 
c 

EFFICIENCY (HEAI    EXChANGER    ADDITION    SYSTcM)     (1=YE3) 

►*# 

*******************************************  ^^SS^**^******^ 

REAL    JK.LC  ,JC,M»TG,NC,NC'»J  ,LMTD,LS  , Kl , MAXLS , LT ,  JL,  J  3 

RtAL    NPVOLU.EFOLO 

REAL    EINF125,25),TDIF(25),NPV(25,25),EFF(25,25) 

INTEGER    KNPV»*»hfc  ,nEPF 

INTEGER .103(25) 

DATA    EINF,TDIF,NPV  ,EFF/  1900*U  ./ 

DATA     105/25*0/ 

READ(5,lpO)    HI.CPI.VISIfTKI 

REAO(5,10Ci     WCiCPOfVISU.TKU 

READ  (5,  125  J'  NEXCH1  ,NEXCH2,I  TQ01  ,1  T0D2  , I TMP4A , I TMP4D , I S I  01 , 

*ISID2>r)NPV  ,H'«E,nEFF 

WRITE16,150) 

WRITE<6, 170) 

WRITE16.200)     Ml ,CPI ,VI3I , TK I 

c 

c 
c 

c 
c 
c 

c 

WR  I TE  (  6  ,  2  20  )    HO.CPOf  ViSO.  TKQ 

******t?***?##**##J$>¥$¥#**t¥*##*****ft4****4fr*««?****¥***(:^* 

INITIATE    OPTIMIZATION    LGLPS 

*********************************************** *********£#*$3 

NXCH1    ■    NEXCHI 

NXCH2    =  .NEXCH2 

I TEMPA    ■    ITMP4A 

IIEMP8   =     ITMP4d; 

c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
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i sai  =  isidi 
iso2  =  1S102 
t»i  ■  nl 

UU    -    WU 
SKIP     =     0. 
IECK    =    0 

NPVCLC    =    -*.*13.**3 
E1NCLO    =    9.*10.*«<* 
JOPT    -    0 
IUPT     =    0 
95    CONTINUE 

*  VARY    THE    NUMBER    Of    HfcAT    EXCHANGERS.  * 

UU    50    NEXCH   =    NEXCHI.NEXCH2 
rtl     =    ThI/NEXCH 
M    -    TwO/NcXCH 

*  VARY    THE    AlLCAOlE    TOdE    CUTSUe    UlAMtltKi     TCU     IN    CMi  * 

*  FUR    TriKtE    PUSSIBLE     VALUES,  b/d,     J/4,    CR     I     IN.  « 

DO    40    JTUO    =    ITCJUITC02 
ITOU    ■    JTUO-2 

»       VAnY     THc    EXIT    PREHEATED    AIR    rEMPERATERi     ITEHP4,     IN  • 

*  UEGktcS    K.       THE    AliSOLUTt    MAXIMUM    IS    533.15    K,     IMG    INLET       * 

*  FLOE    GAS    TEMPERATURE.  * 


I  TMP4A,  ITMP4S,  20 


SI  J,      IN    Crt 


J    ■    0 

DO    10    ITEMP4 

J    =    J»l*JUCT 

IDIFU)     »     S33.15-ITEMP4 

ITN    =    HTMP4d-ITKP4A)/2GH 

*       VARY     THc    SHELL     INSIDE    CIAMETE* 

I    =    0 

UC    20    I  SI  J    =     ISIDI ,IS  U2.2D 

I    ■     H-l  +  IGPT 

SID  •    1S10 

IDS!  I  )     =    I  SID 


H    h    A    T 

EXCHANGER 

DESIGN 

SECTION 

**********«»*«***»»«#»»»*♦♦***»«»»»«*#» ««**»*»»s*«*««»**»««t* 


*  RELATE    THE    TUBE    GLTSIuc    DIAMETER    TO    PITCH,     F,     AND  * 

*  PARALLEL     TO    FLOW    PITCH,    PP.        THE     TOTAL    NUMBER    OF     Ex-  * 

*  CHANGER     UuES,    NT,    ARE    CALCUlaTcD    FrQM    LATA    CORRELATIONS     » 
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IT 

— 

*  ""ACCURATE    MITHIN    3%.                                                                                                    * 

30 

IK  ITCD)     31,32,33 

31 

TCU    =    l.ba75 

P    =    2.0637b 

PP    =     l.7bc)2 

NT    =     lla./5.l»'2l*.^(16'.«IS!J/c.i',l«V.lW/b 

GO    TO    35 

32 

TQO    =    1.905 
P    =    2.3(3123 
PP    =    2.0B76 
NT   =    U*«/3«)»*2)*.40»04*tSl0/2.54J**2.0Jl733 

GO    TO    3b 

33 

TQO    ■    2.b4 

P    =    3.17a 

PP    =     2.7483 

NT    =     .367b4*(;>IQ/2  .b4)«*2  .3544b 

3b 

CU.MTINUfc 

IF(  iECK.NE.O)     GO    TG     132 

WRITE  (6,945) 

WRITc(6,9bO) 

WRITE  (o,9oO)     NXCHl,NXCrl2 

WRl TE  16  ,-WU)     TGJ, TLU 

»RITE(6,9B0)     ISUI.IS02 

«»R[TE  16  ,990)     ITEMi»Af  ITEMPE 

c 
c 

132 

IEoK    "     IECK+1 

«       FINU    THE    AMOUNT    OF    HEAT     TR ANSFcRtU,     ., ,     IN    CAL/FR    FROM    Trtc  * 

c 
c 

*       INSIUE    l-LUIO.                                                                                                                            * 

TEMPI    ■    333.15 

TEMP3    «     31J.93 

TEMP4    =     ITEMP4 

U    »    *1«  IU00./2O.  J*  (o.9»  (TLMP4-T  fcHPi  )  *  I. 92*1  J.»  *(-3)  1/2.  * 

c 
c 

»ITEMP4«»2-IEMP3*»2)  1 

*       CALCULATE    THE    SHElL-SIUE    FLO  LO    EXIT    TEMPERATURE,     TdP2,        * 

c 

«       IN    DEU    K    USING    AN    ENERGY    BALANCE    ANU     A    T„C    PARAMETER    h£.\T» 

c 

c 

*       CAPACITY                                                                                                                                          « 

A   =     1.14/2.«l0.**l-3) 

B    =    7.53 

C    ■    -(B*TcMPl<-A*TEHPl«*2-&/«0/1000.*30.9) 

c 
c 

TEMP2    =     I-d*SURTlB**2-4*A*C)  1/1  2»A) 

*  *  c**    *  *  *  fr  r-  *  ?**  r.  *t?  5*  **  *  *c  **  *  *  t*  fe+  *  *  *  v  s  <-  r  *  **  *  ^  :-  -:■  e  c  *  »  .-  .-'  **  *  4  s*  *  *- * 

c 

c 
c 
c 

CALCULATE    SnELL-SlJt    1-tAT    TRANSFER    PARAMcTtRS 

*******e****9*4:*9********«*****e*t*#******«««4*«*»«#**»*9eS«« 

c 
c 

*       CALCULATE    CRGsSFLCW    AREA,     SM  ,     AT    Uk    NEAR    uENTERLINE    FC.-v       « 

c 

»       ONE    CRGSSFLOw    SECTION    ASSUMING    A    1/2     IN.     CLcARANCE                     » 

c 

*       BETnEEN    SIJ    AND    OUTER     TUBE    LIMIT,    AlSC    FIND    THE     3AFFLE          * 

c 
c 

*       SPACING,     LS.     IN    CM                                                                                                              ♦ 

NNN    =     >J 

Li    =    S1U 

127 

IFUS  .LT.5.03)     LS    =    5.3b 

MAXLS    =    (  /4.»2.54)*(TG0/2.54J  **.75 

IFUS.GT. MAXLS)     LS    «    MAXLS 

OQTL    »    SiQ- .5*2  .54*2. 

c 

c 

47 

SM    =     LS*l  SlD-OOTL+1  jOTl-TOJ)*  IP-TOlD/Pl 

*       CALLULATE    THE    FRACTION    OF    ruTAL     rjdci     IN    CRGSSFlOW,    f-C, 

« 

c 

c 

*       ASSUMING    A    BAf-FLfc    CUT,     LC,     OF    49  5    S I C I ONoBSTkUCTEG    rLCW) 

» 

PI    =    3.1416 

LC    =     .49*SI0 

FCC    =     (S 10-2. *LC 1/DQTL 

c 
c 
c 

FC    =    I. /Pi*  [PI*2.*fCC*SINURCOS  IFoC)  ) -2  .*ARCQS  IFCC)  ) 

*       CALCULATE    ShELL-SiGE    REYNGLJS    NUMBER  t$R£  VN 

* 

SREYN    =    TCJ*WG*10GO./VI SO/SM/3600. 

c 
c 

iPRNDI     =    ^PC*VIS0/TKG 

«       CALCULATE    SHELL-SIUE    H  EA  T- TKANil- ER    CGEFF  IC  I  E  NT  ,  HC  .     IN 

* 

c 

»       CAL/tiEC    SiCM    DEl.    C)     FCK    AN     IGEA    TUdE    dANK. 

* 

c 
c 

c 
c 

*       FIND    JK    FROM    A    TUBE-BANK.    L    REYNGLU    NC.    R  EL  AT  I  CM  R  E>300  .  ) 

t 

jK   -    0»2i3*i»SKtYN»*l-.i4V5) 

*       ASSUME    uuLK.    ANU    WALL    AVERAGE    VISCOSITIES    ARE    ECUAL    SO 

* 

c 

L 

c 
c 

♦       THAT     1 VI SB/VISn)**G. 14    =    1. 

♦ 

HO    =     JK*uPC*WC/SM/3.6UJ*lTKC/CPO/ViSC)*«! Z./5. ) 

*       CORRECT     hO    FCk    bAFFLc-CLNF  I  j  jr  AT  ION    EFFECTS,     JC  ,     LSINo 

* 

c 

c 

*       FC. 

« 

IF1FC.GE..9!    GO    TO     190 

JO    =    .ol25*l-C*G.61  73 

GO    TC     ISS 

c 
c 

190 

JC    =    -FC»2.05 

♦       FINu    THE    CORRECTION    FACILR,     JL  ,     FOR    d AFF LE- LEAKAGE 

* 

c 

*       EFFECTS     ASSUMING    A    1/36     INCH    JlAMcUAL    CLEARANCE     utTnEEN 

v 

c 

«       TUBc    ANl)    BAFFLE    AND    A    1/4     INCH    UIAMETRAL    ShELL-dAFFLE 

* 

c 

*       CLfcARANCE.       THE     TCdE-TC-  bAFF  LE    AND    jHElL-TC-    dAFFLt    LEA*> 

-* 

c 

*       AGE    AREAS    Fun    ONE     UAFFLL    ARL    STB    ANJ    SSd    RESPECTIVELY. 

* 

c 

*       THE    JL    CURVE-FIT    IS    VALlO    FOR    XS.Od     WITHIN    ii    ANC    FOR 

? 

c 
c 

*      X<.33    WITHIN    5J. 

* 

195 

STB    =     .024:>*TOO*NT*(l.*-FC)/2.  54 

SSd    =    ilO/2.54/a.*(PI-AKLGS 11 ,-2*LC/SIJ) ] 

X    =     ( SSB*STB)/SM*2.54 

I    =    SSB/ISSB+STB) 

1FU.LT..03)    GO    TO     196 

JL    =     I-  .4222*Z-.<.7<J52)  *  At  I- .0974*2  +  . 9041  7  J 

GO    TO    197 

!->& 

JL    ■    -2 .02»X*. 973 

c 

*       FINJ    ThL    CoRrECTHN    FACToR,     Jo,     FOR    JUNO  LE-o  Y- PASS!  N  J 

* 

c 

*       EFFcCTS    ASSUMING    NO    SEALING    STRIPS.       THE    FRACTION    OF 

* 

128 

r 

"•"'CRCSSFLOW    AREA    AVAILABLE    FCR    3Y-PA3S     FLOW    IS    F8P.                      * 

i 

147 

FSP    =     (SID-DOR  )*LS/SM 

Jb    =    EXP(-1.23y7*hBP-. 00234) 

\ 

c 

c 

*       CALCOLATE    PINAL    HO.                                                                                                           * 

HO    =    HO*JC*JL*J3 

**#*##«*#****$■((»&****».$##*****$*************«**************** 

c 

c 
c 
c 
c 
c 
c 
c 

c 
c 
c 

c 
c 

c 

.- 

CALCULATE    TUBE-SIDE    HEAT     TKANS^ER    PARAMETERS 

****?#***««**«**#««?**«**  tv****??*****)?*****^*****?***?*****? 

»       CALCULATE    TG3E-SICE    RtYNOLCS    NUMBER ,  TREY  N ,     ASSUMING    A            * 
*       BWG    OF     14.                                                                                                                                    * 

TIU    =    TOO-. 033*2.54*2. 

TREYN    =    TID*Wl*10OO./VISI/lPI*(TIO**2)/'4.)/Nr/2fcOC. 

*      CALCULATE    THE     INSIDE-FLUID    PKAMUTL    NCfBER,     FRNCTL.                     * 

PRNOTL    =    C  P I*V I SI/TKI 
MMM    =    0 

*  INITIALIZE    TUBE    LENGTH     IN    CASE    FLOW     IS    LAMINAR    OF     IN    THE    * 

*  TRANSITION    REGION    AND     ITERATION    IS    REQUIRED.                                     * 

c 
c 
c 
c 
c 
c 

39 

LT    ■     6. 
LTI    =    LT 

*  CALCOLATE    THE    INSIDE    nEA T- TRANSFER    COEFFICIENT,     HI,     USING* 

*  THE    SIEUER-TATE    CURR ELAT IONj .       ASSUMING    VISB/VISn    ■    I.          * 

*  FCR    EITHER    A    TUnBULENT,     TRANSITION,     OF    LAMINAR    REGION             * 

*  DEPENDING    ON     TREYN.                                                                                                            * 

66 
87 

IF(TREYN. LT. 10000. )     GC    TC    36 

HI    =     .023*TM/TID*ITREYN*».  3)  *(  PRNOTL**  U./3.  J  J 
GC    TO    83 

IFITREYN.LT.21J0.)     GC    TC    87 

HI    =     .116*TM/TI0*(  UREYN*»l2./3.  i  )  -  1  25  .  )  *(  PRND  !  L  **  C  1./3.)  1  * 
•I l.M  TID/LT)*l2./3.)  1 
GO    TO    88 
HI    =    1.18*TM/I  IU*(4.*WI*LP  1/  3.  It/TisI/LT  )♦*(  I. /J.  ) 

c 
c 
c 
c 

c 
c 
c 

c 
c 
c 
c 

***♦+♦♦♦**#♦****#***♦**#*** *##*+?**«#****»»*#*#********ff*v?^* 

CALCULATE     THE     TLBING    LENGTH    NElJED    TO 

TRANSFER     ThE    AMOUNT    uF    HEAT,     1,    DETERMINED 

FSCM    THE    TEMPERA  TORE    APPRJACH    CONDITION 

*,n»K;*^Hfn;:^M»;in->fljiJ^f!rf:»;i=ntr^:fi)(*).;Y*- 

*       CAlCULAIE    IriE    LOG    MEAN     TcMPERATURE    DIFFERENCE                                * 

3d 

LMTO    =     (TEMP1-TEMP4-TEMP2*:E1P3)/AL0GI ( rgfPi-TEKJHl / 
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"  *(TEMP2-TEMP3> ) 

t  *      CALCULATE    THE    OVERALL    HfcAT-TKA.NSFER    COEFFICIENT    3ASEU   CN    * 

C  »       THE    OUT  S  IDE    TUdL    AREA,     UC.        FUJLING     IS    ASSUKccJ    NEG-  * 

C     .  *       LIGAtJLE.  * 

C 

TKW   =    .1075 

UO   ■    l./(l./Ha*TOO/HI/Tll3*TOJ*(TOa-TIO»/TKH/HTlO»TOO»/2.ll 

c 

C  ♦       CALLULATE    THE    OUTSIDE    An EA ,    AG,     THE     TUBE    LENGTH,     LT,     AND    « 

L  »       THE    NUMBER    GF    UAFFLcS,    Nd.        IF    MtJ=0    THEN     ITERATE    OY  * 

C  *      REDUCING    THE    BAFFLE     SPACING,     LS.  * 

C 

AO    =    C/UC/LMTD/30GG. 

LT    =     AG/PI/TGU/NT 

IFITREYN.GE.IOGOO. )     GO    Tb     II 

IF(ABSlLTl-LT) .LT..3I  GC  TU  II 

MKM  =  MMM  +  1 

1FIMMN.GE.15)  GO  TO  77 

GO  TO  89 
77  NB  =  LT/LS-l 

NNN  =  NNN*1 

IFINNN.tQ.20)  GO  TG  44 

IF(NB.NE.O)  GG  TU  44 

LS  =  LT/2.-.1 

GC  TO  47 

44  IF(LT.GT.SID)  GU  TO  45 
NPVI  I  ,J)  =  9999999. 
EINFI  I, J)  =  9999999. 
EFFI 1  ,J)  =  .99999 

GO  TO  2  0 

45  IF1LT  .LT.610. )  GC  Tu  46 
NPV( 1  ,J)  ■  tdBdddo. 

einfi  i, j)  =  aoBaeta. 

cFF I  I  ,  J )     =     .odddd 
GO    TG    20 
C 

C 

C  CALCOLATE    SHELL-SIOE    PRESSURE    OKOP 

C 

C 

C  *       FIND     THt     lUEAL-TUBE-dANK    FKICTILN    FACTOR,     F*,     FKCM    A  • 

C  *       RELATION    VALID    FOR    aREYN>300.  « 

C 

46  FK    =    .91o7«SR£YN»» (-.17729) 
C 

C  *       ITTERATE    WITH    RESPECT    TU    T hE    AVERAGE    ShtLL     PRESSURE,  * 

C  «       SPAVG,     TO    CALCULATE    THE    SHtLL-FLUIU    CE.NSITY,    RHUS,     IN  ♦ 

C  »       G/CM»»J.        INITIATE     THE    LCQP    dlTri    A    VALUE    CF     5    PSI     FUR  * 

C  *       THE    SHELL-SIDE    PRESSURE    LRCP    GUESS,     UPSo.  • 

C 

NN    =     1. 

DPSU    =    5. 

MkTO    =    30.9 

AVUfC    »     lTEMPl<-TEMP2)/2. 
3    SPAVG    «     i 14. 7*1 14. 7*UPSGI 1/2. 
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RHOS    =    SPAVG/14.7/(.C<j205)/AVGTO/ LOGO. *M»TG 

C  *      CALCULATE    THE    NUMBER    Of    fu  B£    RJwS    CROSSED    IN   ONE    CRUSS- 

C  »       FLU»    SECTION,    NC. 


NC    =    SID* li.-2.*LL/SI0)/PP 


C  *       CALCULATE    THE    PRESSURE    JKUP,    uPC  ,     lli    PSI     FUK    AN    [DEAL  * 

C  *      CRGSSFLGW    SECTIGN    ASSUMING    VISB/VlSrf    =1.  * 

UPC    =    2.44.014*10. **l-4)*.33613*l  I  1./3.6)  **2)*FK*UG**2)*NC/ 
*RHUS/(SM**2) 


c 
c 

C  *       CALCULATE    ThE    NUMoER    CF    EFFECTIVE    LRGSSFLGW    ROWS    IN    EACH    * 

C  *       hINUQh,    NCW.  « 

C 

HCn  =   .b*lc/pp 
c 

C  »      CALCULATE    THE    ARfcA    FCR    F  LG »    THROLGH    EACH    wIr-OG„,     Srt,     dY       * 

C  *       THE    DIFFERENCE    BETWEEN     THE    GROSS    WINDOW    AREA    A  NO     THc     PART* 

C  *       CF    THAT     ARcA    OCCUPIED    BY    TUBES.  * 

C 

SHC    =    i.-2*LC/SI0 

S«    =     (SlU*»2)/'..*(ARCUSlSnC)-SwC*SjRT(l.-SrtC**2))-NT/8.»(l. - 
*FCJ*P 1*T0G 


C      *   CALCULATE  THE  PRESSURE  JKUP,  uP*  ,  FGR  AN  IDEAL  olNCOw 
C      *   SECTION  IN  PSI. 


UP«    =     .  5**.  01^*10.  *»[-4J*.  Gjbl_t*(  *J**Z)  *(2.+.b*NCW)  /S,*VW 
•RHOS* I (  l./3.b)**2  ) 
C 

C  *       FIND    THE    CORRECTION    FACTOR,    RL,    FGR    EFFECT     CF    EAFFuE  * 

0  *       LEAKAGE    ON    PRESSURE    CROP.       THE    LINEAR    cORVE-FIT     IS    VAlID    * 

C  *       FOR    X<.7    AND    RL>.l.        THE    FIT    APPROXIMATES    THE    CURVES  * 

C  «       WITHIN    54.  * 

C 

IFl X.LT.. 1)     GO     TO    laC 

RL    =     l-.o53l66*Z-.618o/4)  *X*l-.2'<  J3a*Z*.  734491 

GU    TO    IBS 
lao    RL    ■    -3.i*X  +  .<;25 
C 

C  *       FIND     THt    CORRECTION    FACTOR,    Ru,     FOR    BUNGLE    BY-PASS  * 

C  *      ASSUMING    NG    SEALING    STRIPS.       THt    EXPONENTIAL    FIT    ACCURACY* 

C  *       IS    WITHIN     U.  « 

C 

135    Rd    =    EXP(-3.754U<,*FdP-.ClCrC.l  I 
C 

C  *       CALCULATE     THE    PRESSURE    JRuP.    UPS,     IN    PSI     ACROSS     THE    SHLlL* 

C  *       SIDE    EXCLUUING    NUZZLES.  « 

C 

OPS    =    U.Nd-1.  )*UPC*Ra*NE*CP»)  *kt*2*0PC*RB»(  i.  ♦NCw/NCI 

C  *       CHECK,     IF     THE    CALCULATED    PrESoUkE    JauP     IS    nlTHIN     1.    PSI    JF* 

*       THE    INITIAL    GUfcSS.     IF    NOT    THEN    LET     THIS    CALCULATED    VALUE     * 

D  *       BE    THE    NEW    GUESS.  » 

t 

IFINN.GT.luO)    GG    TC     15 
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1 

c 
c 
c 
c 
c 
c 
c 
c 

c 
c 
c 
c 
c 

c 

c 
c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

PUlFF    =    AttSIOPS-UPSG) 
iF(PUIFF.LE.O.l)    GO    TU    15 
UPSG    =    OPS 
NN  *  NN*l« 

Gu    TO    5 

«|«tt»t»»tt*t**«»S**«***t*»t»t**»«**»»*»»*»**»»tt***»t*»»«*** 

CALCULATE    THE     T0i3E-SI0E    PRcSSLRE    DKCP 

»tt(t*»titt*>4tt»tf!»{»fft!ttH(;Hlt(*tmt*tttJ»t»tl»ttl,*!! 

*       FINu    ThE    FANNING    FRICTION    FACTuR.     F,     FROM    TFE    BLASluS             * 
»       FORMULA.                                                                                                                                          * 

15    F    =    .C79I /lTREYN**.25 ) 

*  CALCULATE    Sit    THE    CUKRECTILN    FACTOR     TO    ACCOUNT    FOR    SCDJEN» 

*  CONTRACTION,    SUDDEN    EXPANSION,     AND    REVERSAL    OF    FLOW                  * 

*  DIRECTION.                                                                                                                                     * 

M    =     (  1-(TI0»«2)*NT/SIG»*2)  »*2<-.J5 

FILMOT    =     UEMPl+TEMP2)/2.-l  TE  MP3+ T£*P<V)  /2. 

dl    =     l.*l .51*Kl«FILMCT)/lT£MP4-TEMP3)/lPRNJTL**l2./3. )) 

*       CALCULATE    THE    PRESSURE    UKGP,     TOP,     In     INChES    CF    H20                     * 

MWTI    =    26.8 

AVGTI    »     (T£MP3*TEMP4J/2. 

RHCI    =    l./.Cfl205/AVGTI/1000.»f«T I 

TOP    =    ul»2.»F*UI/  1000. /(  PI  »(  UU**2)M.  1  »»2  )  *L  T/RHO  I /TI  0/1  .  0  1/ 
*24<30.8*(  11  ./36.)**2) 

TUP    =    bI*2.*F«l  1*1  »lCOO./(?I/<f.»T  I  J*«2)/NT)  «»2  J*LT/RHCI/TI  J/ 
•1.02/  (3o0  0.»*2l*4.0  1^2*10.*«(  -4) 

c  c   s   r 

ESTIMATION 
SECTION 

c 
c 
c 
c 
c 
c 
c 
c 
c 

c 

c 

• 

tt,t*n»*t*»utt**¥rt)vtMMtt-»»tjtrntj;r»utj)i*JS)(?^fii 

CALCULATE    ThE    JAN.     17  7v    COSTS    OF     THE    AIR    BLUMtK  i 
AR3CST,    AND    THE    FLUE    GAS     3LC*ER.     FG2GST. 

******»*****♦#*#*******$.*»$********«***♦»***«$****$*#♦******* 

*  CALCULATE    AIR    AND    FLUE    GAS    VoLJMETRlC    FLU»    RATES    ARCFM          » 
«       ANU    FGCFM    RESPECTIVELY,     IN    CJdIC    FT./MlN.    ASSUMING   A.N             » 

*  INLeT    PRtSSURc    OF     1.    ATM.    ALSO    CONVERT    TOP    TC    PS  I .                     * 

AkCFM    =    Wl*1000.*.Od205*TEMP3/Mf,T  I  * .3 35 32/o 0 . « NEXCH 
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FGCf-M    =    w0*lGG0.*.Go2G5*TEMH/M^Ta*.G3532/oG.*NEXCH 

PST    =    .C36L»T0P 

PUS    =    UPS 

c 

c 

♦       USE    THE    BLUnER    CORRELATIONS    FUR    A    PRESSURE     INCREASE    OF 

* 

c 

*       BETWEEN    3    ANU    30    PSI.       IF    THE    PkESSURE     INCREASE     IS    LESS 

• 

c 

♦       THAN    3    PSI     THEN    SET     IT    ECUAL    TO    THAT     LCntR    LIMIT-       HOrt- 

* 

c 

*       EVER,     IF    THIS     INCREASE     IS    GREATER    THAN    _s  0    PSI     THEN    USE    AN* 

c 

*       AXIAL    COMPRESSOR    CORRELATION. 

* 

IFIPST.LI  .3.)    PST    =     3. 

1FI0PS.LT.3-)     PUS    =    3. 

IF(PST.Lb.3C. )    GCIO    60 

ARBCST    =    3.5»10.*»5 

Gu    TO    65 

6U    ARBCST    =     (  13.3847*PST**1.C93W)*ARCFM**  (.  57695*  PST*  * 

*l-. 0170*1)  ) 

65    IFIDPS.LE.30. )    GO     TO    70 

1 

FGBCST    =     6.*10.**5 

GO    TC    75 

70    FGBCST    =     '.  13.38*7*PuS**l.  C93*7)*FGCFM**(  ,57695*PDS** 

*  (-.0  17041)  ) 

c 
c 

75    CUNTI NUE 

*       CONVERT     TFESE    JAN.    1967    CuSTo     lu    A    JAN.     1979    BASIS    JY 

* 

c 

•       USING    A    COST    COEFFICIENT,    CSTCFl,     OfcTERMINEC    FROM    ME 

* 

c 
c 

*       CHEMICAL    ENGINEERING    PLmNT    COST     INDEX. 

* 

CSICFl    =    229. /IC9. 1 

ARBCST    =    ARuCST*CSTCFl 

c 

c 

FGBCST    =    FGBCST'CSTCFl 

*       APPLY    A    FIELU    INSTAlLATILN    FACT-iR    Tu    UBTAIN    A    TOTAL    OUST 

* 

c 
c 

c 
c 
c 
c 

*       OF    THE    BLOnERS. 

* 

6HRCST    =     (ARbCST+FGBi-Sr  1*1.  1 

***************************    **»«***********V*fc*$$«*$*******V*4; 

CALCULATE    THE    JAN.     1979    CCSI 

c 
c 
c 
t 
c 

OF    THE    FEAT     EXCHANGER ( S ) 

•  »»»K,t|l(t(t,,H**l»t>(«HM>tllllltttttt,(t«»>«,„»»|t»J 

*       CONVERT    THE    AREA    TO    SJ.    FT.,     t Au ,     AND    CALCULATE    THE    HEAT 

, 

c 

c 

*       EXCHANGER    BASE    HID     197C    COST,    BASCSF. 

# 

EAO    =    Au/luJUG.*IU.7o* 

c 

c 

BASCST    =     73.7*EAU**.701 

*       CALCULATE    The    DESIGN   FACTOR,     FJ,     FOR    A    FIXEC    TUBE    SHEET 

* 

c 
c 

c 

c 

*       hEAT     EXCHANGER. 

* 

FO    =     .B09*UASCST**l-.  007664) 

*       CALCULATE    ThE    JAN.     1979    PURCHASED    cjUIPMtNT    COST,     hxCljST, 

,* 

c 
c 

*       USING    A    SECONU    COST    COEFFICIENT,     CSICF2. 

* 
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CSTC?F»'    229./I22. 

hXCUST    =    CSTCF2*dA3CSr*F0*NEXCH 

c 

C  *  CALCULATE    THE    TOTAL    MATERIAL    COST,    HXMTL,     USING    A    MATER-    * 

C  *  IAL    FACTOR    UF    0.7.        INCLUDED    ARE    ITEMS    SUuH    AS    PIPING    AND* 

C  «  INSULATION. 
C 

HXMTL    =    . 7*HXC0ST 
C 

C  *  CALCULATE    ThE    TOTAL    HEAT     EXCHANGER    LOST     INCLUDING    LAtJUR,     * 

,C  *  THXCST.     USING    A    MODULAR    FACTOR    OF    3.17.                                                    » 

THXCST    =    3.17*HXCUST 
C 

C  *       U3TAIN    ThE    TOTAL    COST    OF    ThE    HEAT    RECOVERY    ADDITION,  * 

C  *       TOTCST.  * 

c 

TOTCST    ■    THXCST+dHRCST 
C 
C 
C  »»****»***<■»»♦****«****»*»**»*»****♦*****»*******»««********» 

C  t»«M«ii>»tlftltt!*l(»»>tit>ltst«»;»tntlii(ili,t!)tt>t;j» 

c 

C  t    C    C    N    U    M    I    C 

C  A    fl    0 

iC  h    N    E    R    G    Y 

C  OPTIMIZATION 

C  SECTION 

C  ♦»♦**»*»»*««»#**♦***»»«*«»«»**«♦«*«*<■  *♦««♦«««»*««*(.»«♦*  (.»»»** 

c 

c 

C 

C  CALCULATE  ThE  NET  PRESENT  VALJE  OF  THE  MEAT 

C  RECOVERY  ADDITION  FOR  A  12  YEAR  SERVICE  LIFE 

C 

C 

C  *       CALCULATE    TFE    ELECTRIC    POWER    ANNUAL    COST,     P*RCST,  * 

C  *      ASSUMING:       (II     ISOTHERMAL    COMPRESSION,     (2)     i<fO    UAYj/YR  * 

C  *       OF    OPERATION,     1 3  J     An    AVERAGE    ELtCTRIC    RATE     CF     S.Q32/K«H.     * 

C  *       3LG»ER    PChcRS    ARE    IN    CAL/HR,     I t )    dO{     ISOTHERMAL  * 

C  *       EFFICIENCY.  * 

L 

Po«FGS    =    wU*lC0O./MfcTU*l.'*a7*TLMPl*ALOG(  (LPS*l<i.7)/  lt.M/.d* 
•NEXCH 

PCWAIR    =    M  1*1 00  0.  /MHT  l*L98  7*  TEMP  3*  AlOGU  TO  P  +  AC6.S)/<.  Co.  d  )  /  .  d  * 
♦NEXCH 

PWRCST    =     (  PCWA  IR  +  POWFGS)/ JbOO.*<,.  186*13.  **(  -3)  »24.*340.  *.  332 
C 

C  *      CALCULATE    ThE    FUEL    UlL    ANNUAL     SAVINGS,    hLSVNC,     ASSUMING:     * 

L  *        11)     33     A.P.I.     12)     Oi    SuLFJk,     13)     AN    AVERAGE    FUEL    CCST   CF     « 

C  *       i.Bb/GAL    CR    io.l/MM    ETU.        THIS    LEADS     TC    A    hlOFER    HEATING    « 

C  *       VALUE    OF     19370.     dTU/Ld    OR     107oO.CAL/G.  * 

C 

FLSVNG    =    J/2b2. *0.1/(  10. **o  )*2-*.*3<»0. *NeXCH 
C 
C  *       CALCULATE    THE    ANNUAL    COST,     ANLCST,     ASSUMING    ThE    FOLLOWING* 


134 


C 
C 
C 

c 

c 
c 
c 

c 
c 
c 
c 
c 

c 
c 

c 
|c 
c 
c 

£ 

c 
c 

c 
c 
I. 
c 
c 

< 

c 

c 
c 
c 
c 
c 
c 


*  PERCENTAGES  UF  ThE  TOTAL  ADuiTIuN  COST:   MAlNTENANCe-4.SE,* 

*  PLANT  UVcRHEAD-2.6i,  AND  TAXlS  I    I  NSURANCE-  L  .  5  t .   ALSO    « 

*  ASSUME  NO  ADDITIONAL  LA30R  AND  STRAIGHT  LINE  DEPRECIATION* 

ANLCST    =    PnR0SI*(  I  4.*2.6+l.  5) /LU0.H./12.  )*TOTCST 

*  FINJ    THE    NET    ANNUAL    CASH    fCi-nr     AlF,     USING    A    i61    TAX    RATt.* 
ACF    =     IFLSVNG-ANLLST )*l l.-.4d ) tTLTCST/li. 

*  EVALUATE    ThE    NET    PRESENT    VALuE,    NPVi     ASSCMlNu    A    15*    RATE    * 

*  OF    KETUriN,    UNIFORM    CASh    FLOi<A,     NU    WURMNG    CAPITAL,    AND  * 
»       NO    SALVAGE    VALUE.  * 

NPVII  ,J)     =    ACF*l  (EXPl. 15*12  .1-1  )/ .15)*EXP(-.  15»12.)-TOTCST 

***************************************** 4* *««»***S ♦**«*»**** 

CALCULATE     THE    WORK    EQUIVALENT     INTO 

THE    FURNACE    UMT    CN    A    ONE    HOUR    SASIS 

i\cSULTING    FROM    THE    hEAl    RECOVERY    ADDITION 

*  DETERMINE    ThE    nGRK.    ECUIVALENT    OF    ThE     SAVED    FGcL    Old  * 

*  EOIL,     IN    K.CAL.       USE     THE    HI  GH-MOlECUlAR-k  t  IGFT    HYCROCARBON* 

*  ASSUMPTION    THAT     THE    WORK.    ECUIVALENT     IS    APPROXIMATELY  * 

*  ECUAL    TU    THE    HEAT    Or    COM  BUST  ION  ( L  I U.     h.2C  )  .       ALSO    ASSUME        * 

*  OIL    AT    hEFcRENCE    TtHPEKATUkE    AND    PRESSURE.  * 


EOI  L 


-C/  100Q.*NfcXCH 


*  CALCULATE    ThE    WORK    ECUIVALENT    UF     The    ELECTRICAL    WORK.,  * 

*  tELcC     IN    NLAL    ASSUMING    AN    ELECTRICAL    GENERATING  * 

*  EFFGCIENCY    OF    J  d  4  .  « 

EELEC    ■     ( PChAlR'PCwFGSl/lCuC./. 38 

*  CALCULATE    THE    ECUIFMENT    AN  C    CONSTRUCTION     ENERGY    FLCi  ,  * 

*  EOQ.N,     IN    KCAL    FRCM    ENfcRGY/i( 1963)     COEFFICIENTS    ASSUMING  * 

*  NO    SCRAP    VALUE.       USE    A    ThIRD    CUiT     COEFFICIENT,     C S T SF 3 ,  * 

*  TO    CONVERT    TO    THE    JAN.     1979    BASIS.  * 

CSTCF3    =    229./102.3 

EgCON    =    (2J20C.*riXCuST*-207DC.  *HXMTL*l  AK8CST  +  FGBCST  I  *1<W0U  . )/ 
♦  CSTCF3/12./J40./2',. 

*  CALCULA1E     THE    WORK    EqUIVALENT     INPUT     1 1 Tri    RESPECT     TO    T ME       ♦ 

*  ADDITIuN    ASSUMING     INPUT    AIR    AT    REFERENCE    CONDITIONS    FCK       » 

*  THE    COMBINED    HEAT    RECOVERY    SYSTEM.  * 

EAIR    =    TnI/2S.d»CPAl6.6  ,.92  ,TEMP4J 
ESGAS    =     T*0/3U.9*CPA(  7.  53  ,1.1',, TEMPI) 
ESGAS    =     ESGAS*. 013*TWU/J0. 9*13519 
EINF11.JJ     =    EUIL«-EELEC*-E0.CON 

*  CALCULATE    THE    OVERALL     SECOND    LAW    EFFICIENCY    OF     The    MEAT       " 

*  EXCHANGE*    ADDITION    SYSTEM.  t 
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[r 

c 
c 

EFFII.J)       =    EAIR/(ESGAS*EELEC>EuCGNl 

*       ROUTINE    TO    FIND    ANU    SAVE    VARIABLE    VALUkS    CORRESPONDING          » 

c 

*       TO    OPTIMAL    OBJECTIVE    FUNCTION    VALUES.                                                         * 

c 

IF ISKIP.GE.l.)     00    TC    90 

1FINPV1 I, J) .LT.NPVCLO)    GO    TC    60 
NPVOLU    =    NPVt I iJ) 
NJ1    =    J  TOO 

NIA    =     ITEMP4 

Nil    =    ISIG 

NE1    =    NEXCH 

101    =     I-l 

JOI    =    J-1 

dO 

IFIEINK I , J I.GT.fclNOLD)     GO    TO    35 

EINOLD    =    EINMlf  Jl 

EJ1    ■    JTCD 

tIA   =     ITEMP4 

£11    =    ISID 

fcfcl    =    NEXCH 

102    =     I-l 

J02    -    J-1 

d5 

IFIEFF1 I ,J J.LT.tFOLD)     GO    TO    20 

EFOLO    =    EFFl I, J) 

EFJ1    =    JTCD 

EFIA    =     !TcMP4 

EFI1    =    ISlU 

EFU    =    NEXCH 

103    =     1-1 

J03    =    J-1 

10 

CONTINUE 

c 
c 
c 
c 
c 

C 

c 

10 

CONT  INOE 

»*♦*»******»*»****»*«*»*»  *»***«*««*  **«****««*»»«»*J?»«»<i»*»*** 

PRINTCJT    ROUTINE 

»*»»******«*******«*♦♦»«»«♦»«»*»***«♦**»*»«  iitt>tllt*ltlt:>>> 

IKnNPV.uE.l)    GO    TO    120 

1 

-RITElo  ,d50) 

WRITE(6,S6G) 

NR1TE  16,870)    TOO 

«RlTE I6,a74)     NEXCH 

WRITEI6.876) 

WRITE  (6  ,d80) 

WRITEI6.900)     ( IUS( I ) ,1=1, III 

WKiTh(o,910  )     1TDIFCJ), I  NPVUiJI, 1  =  1,11)  ,  J  =  l  ,  ITM 

WRITE  (6  ,94CI 

WKITE(o,930) 

120 

IFIWWE.NE.1)     GC    TO     130 
"RITE (o,920) 
"RITE (6, d6G) 
wRlTE(S>,d70)     TUO 
wRITE(6,a/4J     NEXCH 
W  R  i  T  E  ( 6  , 8  7  d  ) 
WRITE (o ,o60) 

136 

WKlTE(6,900)     II0S(  I )  ,1  =  1,11  I 

WrtITc(6,910)     (TOIF(J)  ,  (EINFl  I  ,J  )  ,  1  =  1,11  )  ,J=1  ,  UN) 

WRITE(6,940) 

nKI TE (6,930) 

130 

IF(rtEFF.NE.l)    GO    TO    40 

WRITE  16  ,925) 

«RITE(6,a60) 

WRIT£(6,ti70)    TOO 

HRITE<6,874)     NEXCH 

wRITE(6,876) 

WRITE(6,880) 

-KITE  16,900)     ( 10SI I)  ,1  =  1,11  ) 

WRITE  (6,890)     (Tl)IF(J),IEFF(I,J)  ,1=1,11)  ,J=1,ITN) 

WRITE  (6  ,935) 

>.RITE(6,S43) 

40 

CONTINUE 

50 

CONTINUE 

90 

SMP    =    SKIP+1. 
IFISK1P.NE.1. )     GO    TO     105 

lFl»NPV.Nt.  1)     GC    TG    140 

IT001    =    NJ 1 

ITGG2    =    NJ1 

1TMP4A    =    N1A 

ITMP48    =    NIA 

1 SIOl    =    Nil 

1SI02    =    M  1 

NEXCH1    =    Nfcl 

NEXCH2    =    Ntl 

10PT    =     101 

JGPT    =    JGL 

wKiTt(6,2aoi 

c 
c 

GO    TO    95 

♦       CONVERT     TO    AMERICAN    ENGINEERING    UNITS    AND    CHECK    vALjES 

* 

L 

*       OF    X ,    RL,    ANO    SREYN. 

v 

u 

105 

IFIX.LE..7)    GC    TO     186 
WRITE16.992) 

136 

IFIRL.GE.. 1)     GC    TO    lo7 
WRITEI6.994) 

137 

1FISREYN.GE.3G0)    GO    TO    ld8 
wRlTE(0,996J 

Ida 

AG    =     AO/ 10000. 

Q    ■    >**NEXCrt 

LT    =     LT/100. 

SOPS    =    DPS/14. o96 

STOP    =    TuP*1.86d 

ESID    =    SIO/2.54 

ELT    =    LT/2. 54/12. *100. 

ELS    =    LS/2.54 

tTGO    =    TOO/2.54 

tHG    =    HO/. 0001355 

6HI    =    HI/. 0001355 

EUO    =    UU/. 3001355 

fcO   =    0/252. 

EIEMP1    =     (TEMPl-273.15)»l.d*32. 

ETEMP2    =     (TEMP2-273.15)*l.d*32. 

ETEMP3    =     I  TEMP3-273.  15)*1. 8*32. 

137 

ETEMP4    =     (TEMP4-273. l5J*A.8+32. 

ElMTJ   =    LMT0*1.8 

EAG    =    AG*10.7o't 

ARFLUW    =    ARCFrt/2113.S 

FGFLCW    =    FGCFM/2113.9 

PUxFGS    =    PGWFGS/6G./252.». 0175a 

POWAIR    =    PCWA1R/60./252.*. 0175a 

AtiPwR    =    PCWAIR*.7457 

FGP»k    =    PChFGS*.7457 

WRITE(6,290 ) 

WR1TE(6 ,300) 

WRITE(6,310)    TCD.ETQD 

WR1TE(6,320)    SID,ES10 

WRITE(6,330)     LT.ELT 

WRITE16.340)    NT 

WRIT£(6,350)    US, ELS 

wRITt(fa,36C)     NB 

WRITElfa,370)    NEXGH 

*RITE(6  ,380) 

WRITE! 6,400) 

WRiTEl6,410)    SREYN 

WR1TEI6.420)     SPRNOT 

*RITE(6,43G)    FO.EHG 

WRITElfa,440)     TREYN 

»RITEl6,450)     PRNOTL 

WRITE<6,460)     HI.EH1 

WRITE16.470)     UC,EUC 

WRITE<6,480)    i.EO 

WRITE(6,490)     TEMPI, ETEMP1 

WRITE(6,500I     TEMP2.ETEHP2 

»RITE<6,510)     TEMP3.ETEMP3 

WR1T£(6,520)    TEMP4 .ETEMF4 

hRITE(6,530)    LMTO.tLMTO 

WRITE(6,540)    AG, EAC 

WRITE16.550)     SOPS, DPS 

»RITE16,56  0J     STOP,  ID P 

WRITE (6,600) 

WRITE(6,610) 

WR1TE(6,620)    ARFLGW.ARCFM 

WRITE  16,630)     FGFL.G  W  ,  FGCFrt 

WRITE16.640)     PC*AIR,A6P»R 

WRlTE(6,obO)    PCHFGS.FGPWR 

■  Rl  TE  (o  ,660) 

WRITElo, 67C) 

WRITE(6,o80J    ARbcSI 

kRITE I6,o90)     FG3CST 

WRUE16.700)    riXCOST 

WRITE  (6,710)     HXMTL 

wRlTE(6,720J     THXCST 

WRITE(6,730)     TGTCST 

«RITE(6  ,740)     PWRCST 

wfilTE(o,75U)    FLSVNG 

WRITc(6,760)    ANLCST 

WRITE(6,770)     Ai.F 

WRITE(6,780J    NPV(  I ,J) 

WRITE  (6,790) 

»RlTt(6,d00) 

WRIT  t(o, 810)     ECJIL 

138 

WRITE(6,d20)     EELEC 

WR1TE(6,830)    EyCCN 

*RITE(6,8<tO)    EINFl I , JJ 

WRITE(6,845)     EFF1I.J) 

IW 

IFlSKlP.Eg.ll     SKIP    =    2 

IFISKIP.NE.2.)    GO    TG    16C 

IFtWWE.NE.U     GG    TO    160 

ITC01    =    tJl 

1T0U2    =     EJl 

ITMP4A    =    EIA 

ITMP43    =    ElA 

IS  101    ■     Ell 

1S1U2    •    Ell 

NEXCHl    =    EEl 

NEXCH2    =    EEl 

lUPT    =    1U2 

JOPT    =    JG2 

■  R1TE16  ,2d5) 

GO    TO    9  5 

loO 

IF1SK1P.E3.2)     SKIP    =     3 

1FISKIP.NE.3)     GG    TC    999 

IFlWEFF.NE.l)    GO    TO    S99 

ITUU1    =    EFJ1 

1 TUD2    =    EFJ1 

ITMP4A    ■     i.F  IA 

ITMP4E    =    EFIA 

1S101    =    EFll 

IS1D2    =     EFll 

NEXCHl    =    EFE1 

NEXCH2    =    tFfcl 

I  OPT    =103 

JUPT     *    JU3 

WRlTE(fa,287l 

c 
c 
c 
c 
c 
c 

GO    TC    95 

l>>(tl>>t>>l(l>tU>t>««li 

i **  **  «***  «*  ****  ****  ***  ***  ******  ****** 

FORMATS 

************************( 

it*********************************** 

c 

100 

FCRMAT18F10.3  J 

12b 

FU«MATU1I5J 

150 

FGRMATI//////////56X, "EO 

10    CHECK    OF    DATA'////) 

170 

FORMA7l<t3X,'MASS    Ft_GW    RA1 
»5X, 'THERM.    CCNU.  '/47X.  '  Kl 
»6X, 'CAL/S/CM/OfcG    K'//J 

t'  t2X,  'HEAT    CAPACITY', 4X i  WISCOSITY'i 

./HR'  ,7X,  'KCAL/KG/UEG    K  '  ,  5X  ,  •  G/CM/ S  •  , 

200 

FCRMATI2/X, ' INSIDE    FLUID 

i5X,Mtl2.4,3X  )//) 

220 

FGRMAT127X,  'GLTSIOE    FLUIC 

>',<iX  ,ME  12.4>3X)  //) 

2d0 

F0RMAT('1'///<.0X  ,'  **♦*»    5 
•'OPTIMUM    ««*»♦•////) 

EARCH    RESULTS    FOR    ThE    ECONOMIC     ', 

2di 

FORM ATI' l'///,3  7X, ****** 
»'K    EQUIVALENT    »**»*'//// 

SEARCH    RESULTS    FUR     THE    OPTIMUM    WGR '  , 

287 

FORMAT!  •  IV//30X,  •**»*«    . 

.EARCH    RESULTS    FCR     THE    CPTIflUM    CVE.-C, 

►'ALL     SECOND    LAW    EFFICIENC 

.r     »****•////) 

290 

FGnMATl///55X,'HcAT-EXCH/ 

iNGER    UtSIGN'  ,///  J 

300 

FGRMATi 31X,  'hfcAT-EXCHANGt 

.R    SPECIF  ICATlGN',2oX,'  VALUE'//) 

310 

FCRMAT(32X ,' TUBE    GUTSIOE 

CI  METER'  ,  2<»X,F5.3  ,IX,'CM'  .8X.F5.3,  IX 

139 

*• IN'/) 

320    POHMAI432X,  'SHELL    1NSI0E    Ci  AMET  eR  ■  ,ZZX,  F  7  .3  ,1X  ,  •  CM'  ,6X,  F6.2  .  2X 

*• IN'/  ) 

330    F0HMAT(32X,'TUBE    L  ENGTH'  ,  34X  ,  Fi  .3  ,  1  X  ,  >H  '  .  oX  ,  F5  .2,  2X  ,  •  FT  '  /  ) 

340    FURMAT132X, 'NUMBER    OF    TUB ES ' , 27X , I  4/ ) 

350    F URM AT  132 X,' BAFFLE    SP AC  IN G '  , 29X , F 7 . 3 ,  IX ,  'CM  ' , 7X ,F 5. 2 , 2X ,  •  IN ' / ] 

360    FURMAT132X, 'NUMBER    OF    B AF FL ES ' . 27X , 12/ > 

370    FORMAT  I  3.2Xi  'NUMBER    OF    HEAT    EXCHANGERS ',  19X,  I  2/ // ] 

3d0    FURMAT!  '  1'  //  ,55X,  'HEAT-EXCHANGER    UESIGN',///) 

400    FORMAT! 7X , 'HEAT     TRANSFER    ANC/CR    PRESSURE    CROP    *UAM  I  T  Y  '  , -UX  , 

•'VALUE'// I 

410    FCRMAT!  IOX  ,'SHELi.~SICE    RtYNCLC'S    NUM  EER  '  ,24X  ,  i  12  .5/ ) 

420    FORMATUOX,  'SHELL-SIDE    PRANOTL    NUMBER  •  ,  2  7X  ,  F  6.  4/) 

430    FCRMATUOX,'ShELL-SICE    HE  AT -T  RANS  FER    CO  EFFI  C  I  EN  T  •  ,  1  6X  ,f  9 .  7  ,  3  X  , 

»'CAL/S/Sg    CM/C  , 12 X , Fo. 3 , 7X  ,  • BTu/HR/SC    FT/F'/J 

440    FORMAT!  LOX,  'TUBE-S1UE    REYNOLO"S    NUMBER  '  ,  25X ,  E  12.  5/ i 

450    FCRMAT110X,'TUBE-SIOE    PRANUTL    NUM BcR • , 2BX , F6 .4/ ) 

460    FORMAT!  10X,  'TUBE-SUE    HEA  T- TRANSFER    COt  FF  IC  I  ENT  '  ,  17X  ,  F9  .  7  ,  3X  , 

♦'CAL/S/SS    CM/C ,12X, F6.3, 7X, 'BTU/HR/SO    FT/F'/J 

470    FORMAT!  10X, 'OVERALL    HEAT     TRANSFER    CCEFF I C IENT •  ,  19X, F9  .  7 , 3X, 

*'CAL/S/SO    CM/C  ,  12X.F6.3, 7X  ,  • B I U/HR/ SU    FT/F'/J 

480    FORMATUOX,'  TOTAL    AMOUNT    lF    HEAT     TRANSFER  EC,  2CX  ,  E  12.  5,  IX  , 

*'CAL/Hk' ,19X,E12.5,1X,'BTL/HR'/) 

490    FORMATUOX,  'INLET     FLUE    GAS     TEMP  ER  ATURE  '  ,  24X  ,F  6.  2  ,  EX  ,  '  JE  G    <■, 

*19X,F6.2,SX,'0EG    F'  /) 

500    FORMAT!  IOX,  'OUTLET    FLUE    GAS    T EM PE RATURE ', 23X  ,F 6. 2  ,  3X ,' 0  EG    K '  , 

*19X,F6.2,8X,'DEG    F '/ ) 

510    FORMAT!  IOX, 'INLET    CCMBUSTIUN    AIR    T EMPERATCRE •  ,  18X , Fo  .  2 , 3X , 

*'OEG    K' ,19X,F6.2,8X, 'OEG    F'/) 

520    FORMAT! IOX, 'CUTLET    COMBUSTION    AIM    T EMPERATURE • , 17X, F6 .2, 3X, 

*'DcG    K' ,  19X,F6.2,8X, 'UEG    F'/J 

530    FCRMAI UOX.'LOG    MEAN    TEMPERATURE    J  1  FF  ER  ENCE  '  ,  1  9x  ,  F  6  .  2  ,  BX  , 

»*DLG    C  ,  19X.F6.2  ,dX,'OEG    F'/l 

540    FORMAT!  10X,  'OLTSlLiE     TUBULAR    ARE  A/  EXCHANGE  k  '  ,  2  1  X  ,  E  11  .  4  ,2  X  , 

*'SQ    M'  ,21X  ,EU.4,2X, 'S,    FT'/J 

550    FCRMAT! IOX, 'SHELL-SICE    PRESSURE    JriOP  '  ,2  7X  ,  F  7  .4  ,6X  ,  '  AT  A'  ,  2  IX  , 

♦F7.3.7X,  'PSI'/J 

560    FORMATUOA,'TUBt-SIOE    PRESSURE    URUP  '  ,  2o  X  ,  F  7  .2  ,  dX  ,  'r.M    hG',l-yX, 

*F7.3,  7X,  '  IN    H20*/i 

600    FCRMAT  1////59X,'  cLCHER    DESIGNS'///] 

610    FORM ATI 35X , • SPEC  I F I CA T I CN '  , 3* X , ' V ALUE ' / / 1 

620    F0RMAT131X, 'AIR    FLG»    RA  TE  '  ,  23X  ,F6.  2  ,  1  X,  •  C  C    K/S ' 1 7%, El 1.4 ,IX, 

*'CFM'  /) 

630    F0RMAT131X, 'FLCE    GAS    FLCi.    R  AT  L  •  ,  1  dX  ,  FB  .  2  ,  IX  ,  •  CI    M/S  •  ,  7X  ,  E  1 1  .  4  , 

♦  lX.'CFM'/J 

640    F0RMAT!31X,'AIR    BlCiER    FOwE  R'  ,  1  /X  ,  F*  .  2  ,  IX  ,  '  KW  ,  7X  ,  F  9  .  2,  7X,  •  HP  ' 

*/) 
650    FCRMAT131X,  'FLUE    GAS    SLOWER    PCwER  '  ,  12X  ,  F9  .  c  ,  IX  ,  '  Ki, '  ,  7X  ,  F  9.  2  .  7X  , 

*  'HP '  ) 

06O    FORMAT!  '!•///, 5CX,  'F.EAT    RECOVERY    AODiTlCN    E  CC  NCf  I  C  S  '  ///  ) 

670    FGRMAT147X, 'CUANTITY' ,25X,' JAN.     1979    VALUE'//) 

6d0    FORMAT!  36  X, 'AIR    BLC#£R'  ,3oX  ,'  f  ,F9.J/  J 

690    F0RMAT(3oX, 'FLUE    GAS    BLOWER ', 32 X , F 9. J/J 

700    FURMAT!36X,'HtAT    EXChANGE  R'   ,  j  3X  ,  F  9  .  0/ ) 

710    FORMAT! 36X, ' TOTAL    HEAT    EXCHANGcR    MA T ER I AL '  , 1 dX  ,  F9  .0/ ) 

720    FORMAT 136X, "HEAT    EXChANCER    ML  DOLE ' , 26X , F9 . 0/ I 

730    FORMAIl  J6X,'HEAT    RECOVERY    A  CO  I  T  ION  ',  25a  ,  F9  .  0/  ) 

740    FCRMAT!  36/,  'A/ERAGE    ANNUAL    PO  hE  R  '  ,  2  7x  ,F  9.  0/  ) 
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"750_F0RMAT(3oXr' ANNUAL    FlEL    SAVINGS'  ,2dX,F9.C/) 
760    FCRMAT136X, 'ANNUAL    CCST ' , 36 X, F9 .0/ ) 
7/0    FORMAT! 36X, 'NET     ANNUAL    CASH    F  LOrt  •  ,  2  7X  ,F9  .  0/ ) 
730    FGRMAT136X, 'NET    PRESENT    V  AL  Ufc  '  ,  28X,  Fll.  0/  ) 
790    FGfiMAT(////49X,'iiCRK    EQUIVALENTS     (ONE    HOUR    BASIS)'///] 
800    FGRMAT148X,  'QUANTI  TY  '  ,23X,'  VALUE     IN    KXAL'//) 
dlO    FCRMATUOX, 'OIL    S  AVEC  '  i  3UX ,  E 12  .  5/  ) 
820    F0RMATl40X,'eLECTRlClTY'  ,2dX,£12.  5/  ) 
d30    FCKMAT140X , 'EQUIPMENT     i    C  CNST  RUCT  I  UN  '  ,  1  5X  ,  E  L2  .  5  /  J 
840    FORMAT  1 40X, '  TLTAL '  ,34X , El  2 .  5/  ) 

d45    FCRMATl////45Xi 'OVERALL    SECCNC    LA*    EFF I C I ENCY '  ,3X ,F8 .5 ) 
850    FCRMATl  ■  1'///  ,44X,  'NET    PRESENT    VALUES    FOR    A    RANGE    OF    TEMPER', 

•'ATURES' /) 
860    FCRMAT146X, 'AND    SHELL     INSIDE    DIAMETERS    xITH    A    CCNSTAM'/J 
870    FCRMAI150X, « TUBE    OUTSIDE    DIAMETER    GF',Fo.3,   '    CM1/) 
d74    FCRMATl 55X, 'AND' ,13. '    HEAT    EXCHANGERS'/) 
876    FGRMATI63X , '{ JAN.     1979     i)'///) 
87a    F0RMAT161X ,' IKCAL/hR) '/// ) 

880    FORMAT!  7X,  'TtrtP.  '  ,  <,5X  ,  '  SHELL     INSIDE    JIAMETEk     (  Cf )  >  I  7X  , '  D  IFF  .  ■  j 
890    FLRMATlGX,F6.2,3X, 11F1G.5/) 
900    FORMAT!  /X.'DEG    C •  , 7X , 1 11  I  3 , 7X ) / ) 
910    F0RMATl6X,F6.2,3X,UF10.0/) 
920    FCRMATl' 1' /// ,45X, 'WORK    EQUIVALENTS    FUR     A    KANCt    CF     TEMPERA!', 

*'URES  '/) 
925     FCRMATl ' 1 '///37X, 'OVERALL     SECOND    LAW    EFFICIENCIES    FCR    A    RAN', 

»'GE    OF    TEMPERATURES'/) 
930    FORMAT120X,  'ERROR    CONDITION    9S99,)99.:        LENGTH    CF    TUEES    LESS', 

*•     THAN    SHELL     INSIDE    DIAMETER') 
935    F0RMAT(////2CX,'ERRCR    CCNSITICN    .djJbd:        LENGTH    CF    TUBES     ', 

*• EXCEED    6.1    M    I  20    FT )  •//) 
940    FGRMATI////20X, 'ERROR    CCNU1TICN    aaddtjfld.:        LENGTH    OF     TUQeS     ', 

♦•EXCEED     6.1    M    (20    FT)  '//) 
943     FCRMAT120X,  'ERROR    CONDITION     .99999:        LENGTH    CF     TUocS    LESS     ', 

*'THAN    SHELL    INSIDE    DIAMETER') 
945    FORMAT! 5dX,  'RANGE    CF     SEARCH'///) 
950    FuRMAT(35X, ' INDEX    VAR  1  ABLE'  ,2  5X  ,  '  IN  1  T  I  AL    V  ALUE  '  ,  8X  ,  '  F  In  AL     ,"  , 

»'ALUE '//) 
960    F0RMATI26X,  'NUMBER    OF    HEAT     EXCHAN oERS • , 2 5X , I  3 , 1 6X  ,  I  3 / ) 
9/0    FCRMAT126X  ,' TUdE    OUSIDE    D I AMETcR' , 1 IX , F5 . 3 , IX , ■ CM ' , 12X, F 5. i , IX, 

«'CM  •/  ) 
9d0    FCRMAT126X,  'SHELL     INSIDE    CI  AM  ET  EK ',  28  X  ,  i  3  ,  5  X  ,' C."'  ,  1  CX  ,1  3  ,  5X  , 

»'CM'/) 
990    FURMAT126X,  'CUTLET    CCMBLSIICN    AIR     TEMPE RA TJkE ',  loX , I  J , 5 X ,' K •  , 

*11X,  13, 5X,  •«.'///) 
992    FORMAT!'     X    GREATER    THAN    ./') 
994    FGRMAT!*     RL    GREATER    ThAN     .1') 
996    FCRMATl'     SREYN    LESS    ThAN    800') 
999    STOP 
END 
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FUNCTION    CPA(A,S,T) 
C 

C 

C  THIS    FUNCTION    SUBPROGRAM    CALCULATES    THE    CHANGE 

C  IN    WORK.    EQUIVALENT/nMOL    AS    A    COMPONENT    IS 

C  ADJUSTED  TO  THE  STANDARD  STATE  (298.15K) 

C 

£  ******  ************$*******$**«$**$***************♦*********** 

C 

c  =  a/iooo. 

TI    =    2S13.15 

CPA    =     (  A-T  I*C)«(T-Tl)+C/2»{  T**2.-TI**2.  )-A«TI*ALCOl  I/II  I 

RETURN 

END 
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HEAT    RECOVERY    AOCITICN    ECCNCMICS 

QUANTITY 

JAN.     1979    VALUE 

AIR 

3LCWER 

$          8  10  05. 

FLUE 

GAS    BLCWER 

115593. 

HEAT 

EXCHANGER 

654146. 

TOTAL    HEAT    EXCHANGER    MAT 

ERIAL 

4579  02. 

HEAT 

EXCHANGER    MODULE 

2073642. 

hEAT 

RECOVERY    ADDITION 

22  398  59. 

AVERAGE    ANNUAL    POWER 

2  616  99. 

ANNUAL    FUEL    SAVINGS 

257C8E4. 

ANNUAL    COST 

6380C6. 

NET 

ANNUAL    CASH    FLOW 

1195921. 

NET 

PRESENT    VALUE 

43650  C6. 

WORK    EQUIVALENTS 

ICNE    HOUR 

BASIS) 

QUANTITY 

VAL 

,UE    IN    KCAL 

OIL    SAVED 

-0. 

13U16E    C3 

ELECTRIC! TY 

0. 

22682E    C7 

EQUIPMENT    S    CONSTRUCTION 

0. 

.  12566E    C6 

TOTAL 

-0. 

10622E    08 

OVERALL    SECC.ND 

LAW    c 

FFICIENCY 

0.32135 
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HEAT    RECOV 

ERY    AODI TICN    ECONOMICS 

QUANTITY 

JAN.     1979    VALUE 

AIR 

BLOWER 

$          810C5. 

FLUE 

GAS    BLOWER 

115593. 

HEAT 

EXCHANGER 

1156419. 

TOTAL    HEAT    EXCHANGER    MAT 

ERIAL 

3094  93. 

HEAT 

EXCHANGER    MODULE 

3665847. 

1-EAT 

RECOVERY    A  DO  IT  ION 

38821 C4. 

AVERAGE    ANNUAL    POWER 

182100. 

ANNUAL    FUEL    SAVINGS 

257C364. 

ANNUAL    COST 

82C059. 

NET 

ANNUAL    CASH    FLOW 

1233937. 

NET 

PRESENT    VALUE 

29d4347. 

WORK    EQUIVALENTS 

(ONE    HOUR 

BASIS) 

QUANTITY 

VALUE    IN    i<.CAL 

CIL    SAVED 

-0. 

13016E    C8 

ELECTRICI  TY 

J. 

15733E    C7 

EQUIPMENT    &    CCNSTRUCTIQN 

0. 

21203E     Co 

TOTAL 

-3. 

11225E    08 

OVERALL    SECOND 

LA'*    EFFICIENCY 

0.34313 
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HEaT  recovery  audition  eccnchics 

cuantity  jan.   19  79  value 

air  slower  i       310c5. 

flue  gas  blcwer  u5593. 

heat  exchanger  1156419. 

total  heat  exchanger  material  809493. 

heat  exchanger  module  3665347 . 

heat  recovery  addition  38321c4. 

average  annual  power  182100. 

annual  fuel  savings  257c3t4. 

annual  cost  820059. 

net  annual  cash  flow  1233937. 

net  present  value  ^934347. 


WORK    EQUIVALENTS  ICNE    HUUR   fiASIS) 

yUANT ITY  VALUE    IN    *CAL 

OIL    SAVED  -0.13016E    C8 

ELECTRICITY  0.15783E    C7 

EQUIPMENT    &    CONSTRUCTION  0.21203E     C6 

TOTAL  -0.11225E    08 


OVERALL    SECOND    LA«    EFFICIENCY  0.34313 
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ABSTRACT 

The  basis  for  second  law  thermodynamic  analysis  of  processes  involving 
chemical  transformations  is  reviewed  and  the  COED  process  and  its  attendant 
chemistry  is  described.   The  energetics  of  this  process  and  of  each  step 
therein  are  calculated  from  a  commercial  plant  design  study  and  are  examined 
from  a  second  law  perspective.   Two  efficiencies,  an  incremental  and  an 
absolute,  are  found  to  be  useful  in  this  analysis.   The  incremental  effi- 
ciency is  most  useful  with  physical  processes  while  the  absolute  or  second 
law  efficiency  is  most  useful  with  processes  involving  chemical  transfor- 
mations.  The  second  law  (absolute)  efficiency  is  shown  to  be  0.75.   In 
contrast,  the  thermal  or  first  law  efficiency  is  0.80.   When  the  ineffi- 
ciencies of  oxygen  production  and  electricity  generation  are  charged  to 
the  COED  process  the  second  law  efficiency  decreases  to  0.68.   An  energy 
cost  assigned  to  the  equipment  and  construction  is  evaluated  using  energy/ 
dollar  data  obtained  from  an  economic  Input-Output  analysis  for  the  U.S. 
economy.  When  this  energy  cost .plus  an  energy  credit  for  salvage  is  in- 
cluded in  the  analysis  a  final  overall  second  law  efficiency  of  0.67  is 
found. 

In  addition,  the  optimization  of  a  flue  gas  heat  recovery  system  from 
a  second  law  perspective  is  performed.   This  system  is  described  and  two 
thermodynamic  and  one  economic  objective  function  are  determined.   These 
objective  functions  are  optimized  from  the  standpoint  of  four  independent 
design  variables,  the  number  of  heat  exchangers,  the  tube  outside  diameter, 
the  shell  inside  diameter,  and  the  exit  preheated  air  temperature.   Results 
are  compared  on  the  basis  of  the  number  of  heat  exchangers.   The  thermo- 
dynamic optimums  are  found  at  approximately  30  heat  exchangers  but  are  taken 
at  13  heat  exchangers  since  the  objective  function  values  of  these  two 
amounts  are  within  5%  of  each  other.   In  contrast,  a  distinct  economic 


optimum  is  found  at  8  heat  exchangers.   The  design  conditions  of  the  eco- 
nomic optimum  are  determined  to  produce  a  thermodynamic  objective  function 
value  within  20%  of  the  optimal  value.   The  converse,  however,  is  found 
to  not  be  true. 


